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All manufacturing processes are shape realization processes. How the shape is 
realized is as important as how it is designed. The present thesis addresses the problem 
of modeling the shape realization processes mathematically. Unified models for the 
geometry of a generic cutter as well as for the kinematic structure of a generic machine 
tool have been proposed using the concepts of envelope theory and conjugate geometry. 
The methodology adopted here is called conjugate geometry formulation. These models 
of the generic machine tool and the generic cutter have been implemented as symbolic 
algorithms using the symbolic manipulation package MACSYMA for modeling a variety 
of shape generating processes. It has been shown that several types of machine tool 
kinematic structures can be realized using the proposed unified model of the machine tool. 
Also, it has been shown that a variety of cutter geometries can be instantiated from the 
unified model of the cutter. Furthermore, one can automatically obtain the equations of 
the geometry of the machined surface for any combination of a cutting tool geometry and 
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a machine tool structure. The proposed methodology enables a manufacturing engineer 
to evaluate the likely outcome of a particular cutting tool/ machine tool assembly of a 
manufacturing process. 

Any shape generation process has three constituents, viz., the geometiy of the 
cutter, the geometry of the surface it produces on the blank and the relative motion 
between the cutter and the blank. The conjugate geometry formulation is based on the 
fact that the surface of the cutter and the machined surface are mutually enveloping 
surfaces. A coordinate frame is attached to each of these two surfaces and the positions 
and orientations of these two moving frames are given by twelve parameters with respect 
to a fixed global frame of reference. These twelve motion parameters are algebraic 
expressions as functions of time and the cutter surface and the blank surface will be 
expressed as algebraic expressions in the bi-parametric form. Homogeneous 
transformation matrices are used to map the instantaneous positions of the cutter surface 
and the blank surface from one frame to the other. Out of the three constituents of a 
machining process, if any two are known, die third constituent can be determined using 
the mapping functions and the condition of cutting. The condition of cutting is the 
orthogonolity condition between the relative velocity vector and the common normal of 
the cutter and blank surfaces at the point of cutting. The proposed model of a generic 
machine tool will be a useful tool for NC cutter path simulation, NC cutter path 
generation and selection of an appropriate combination of a machine tool and cutting tool 
to produce a given surface. In the traditional methods for cutter path simulation and 
cutter path generation, the blank is considered to be fixed and all the relative motions are 
imposed on the cutter by superimposing the reverse motion of the blank to both the cutter 
and the blank. In contrast to this, the degrees of freedom of both the cutter and the blank 
are directly used in the proposed conjugate geometry approach and hence when it is used 
for cutter path generation, the output is in terms of the joint coordinates of the machine 
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tool thus alleviating the need for post-processing and when it is used for cutter path 
simulation, the surface obtained is more realistic and free from post-processing 
inaccuracies. Added to this, the output of the symbolic programs will be algebraic 
expressions which can be reduced to numeric form of any desired accuracy. 

The models of the generic cutter uses nine parameters to define a two-dimensional 
profile. This profile consists of a set of three straight line and a circular arc segments. 
By appropriately sweeping this planar profile, the bi-parametric definition of the cutter 
surface is obtained. 

The usefulness of the two models of the generic machine tool and the generic 
cutter have been demonstrated by modeling a number of manufacturing processes. 
Modeling of gear manufacture using a gear planing process with a rack cutter and using 
a gear bobbing process has been illustrated. The proposed methodology will be helpful 
for a manufacturing engineer to evaluate alternative machining strategies in order to 
choose the optimum one. By simply changing the values of the nine parameters 
describing the cutter surface and the twelve motion parameters, an engineer can obtain 
the machined surface and can visualize it using a graphic display. Three alternative 
strategies to produce a helicoidal surface have been modelled to illustrate this. Two 
different ways of producing a hyperbolic paraboloidal surfaces are also given. Finally, 
a cylindrical cam and follower mechanism has been modelled. For a given cylindrical 
translating follower, the surface of the slot on the cam has been obtained using the 
proposed models. It is interesting to note that the same set of input and output will 
describe the manufacture of the cam surface using an end milling cutter. 

An extensive case study has been carried out for the design and manufacture of 
extruder screw surfaces. This case study establishes that the methodology proposed in 
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this thesis for modeling the manufacturing processes can also be useful for the 
geometrical design of extruder screws. First, the procedure for the design of one of the 
surfaces of the extruder screw pair is explained by considering it to be a set of contiguous 
arc segments. Since the relative motion between the screws is known and both the screw 
surfaces are mutually enveloping, the other screw surface can be determined using the 
generic model of a machine tool with twelve degrees of freedom. Examples for the 
design of extruder screws of all three types, viz., symmetrical screw pairs, unsymmetrical 
screw pairs of equal diameters and unsymmetrical screws of unequal diameters have been 
given. 


It is found that only the root and flank arc segments of an extruder screw influence 
the cutter design for the manufacture of the extruder screws since the tip segment is 
realized in the preparation of the cylindrical blank itself. A procedure to determine the 
characteristic cutter profile for any screw surface geometry has been developed. This 
characteristic profile is the intersection curve between the helical slot of the screw and a 
plane that is normal to the helix. The axial position of this plane is chosen such that the 
full slot is captured in it. The extruder screw which is basically a helicoidal surface can 
be machined using a turning process with a form tool. It can also be produced by a 
milling process with an end milling cutter, a side-and-face milling cutter or by a side 
milling cutter. The procedure to obtain the surfaces of the turning tool, end milling cutter 
and disc type of cutter have been explained with the help of examples. The design of the 
side milling cutter has been studied in detail and it was found that except for very simple 
helicoidal surfaces, such a cutter does not exist. However, a procedure has been given 
using which a side milling cutter can be designed if such a cutter exists and this 
procedure has been explained with an example. 
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Chapter 1 


INTRODUCTION 


1,1 Geometric Modeling in Design & Manufacturing 

Until recently the primary focus in geometric modeling has been on modeling of 
the basic geometry' and topology of a component. As design engineers continue to evolve 
powerful definitions of free form surfaces using these geometric modeling techniques, the 
problem of manufacturing such surfaces becomes a challenging task. Hence emphasis is 
now shifting towards the development of computer-based models of product shapes 
which take into account aspects of physical realizability in addition to its basic geometry 
and topology for representing a component. Computer-based models which take into 
account the shape generating capabilities of machining processes are essential to evaluate 
the manufacturability of a design. What is now referred to as "Design for Manufacture", 
"Concurrent Design" or "Simultaneous Engineering" is increasingly becoming relevant to 
the total life cycle value of a component. 

Physical realizability of a designed shape is a primary concern of any 
manufacturing activity. If the surface is to be realized using a milling-like process then 
it is necessary to select the most appropriate geometry of the cutting tool as well as the 
most suitable kinematic structure of the machine tool or the machining center. In some 
cases, a manufacturing engineer may have alternative options of cutting tool geometries 
and machine tool structures available for realizing a given designed surface. In such a 
case, it is necessary to find, before hand, the most appropriate combination of a machine 
tool type and a cutting tool geometry. Furthermore, a manufacturing engineer is usually 
called upon to define the law of relative motion between the cutting tool and the blank 
for a given type of kinematic structure of the machine tool. This is necessary since many 
a machines are now being controlled using CNC controllers. All such problems are now 
being considered as an important class of problems in the area of design for manufacture 
(Grossman, 1986; Dixon, 1989a, 1989b). Design for manufacture involves many aspects 
such as ease of manufacture, process planning, ease of inspection, design of tolerances. 
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design for reliability, testability and maintainability. Research interests in this area are 
growing. Many researchers are using the feature based representations to carry out these 
tasks of designing. 

Design, analysis and manufacturing constitute three basic components of any cycle 
of product development. The process of design in general, and engineering design in 
particular, requires simultaneous considerations of engineering analysis as well as 
considerations of realization or manufacturing. As the technology of Computer Aided 
Design started becoming available at a reasonable cost to a designer, it became obvious 
that the developments would be directed to make the CAD tools more effective and 
relevant to the needs of a designer. After an initial phase of development where 
considerations of geometric modeling played a significant role, it appears that this initial 
development has triggered research in several different areas of design and manufacturing. 

Behavioral scientists, computer science experts and engineers are trying to develop 
an appropriate design theory and methodology (Finger and Dixon, 1989a, 1989b). 
Descriptive, prescriptive and cognitive models are being used for developing such design 
theories. Engineers are trying to model all stages of the design process, viz., conceptual 
design, configurational design and detail design. There is also a considerable effort in 
developing new languages for design and manufacture. In the area of engineering analysis, 
considerable work is imderway to develop interfaces between analysis and design 
processes. Many designers are incorporating different methods of optimization for 
improving a design. A quick review of the literature (Wozny et al., 1990) shows that 
design implications of manufacturing have started to emerge as an iinportant consideration 
in the process of engineering design (Figure 1-1). 

In order to be able to visualize and evaluate the geometry of machined surfaces 
generated by employing alternative machining strategies, unified symbolic models and 
algorithms for conventional milling-like processes have been developed in the present 
work. In these models, the cutter is considered to be a biparametric axisymmetric surface. 
The machined surface is assumed to be a conjugate surface generated by the cutting tool 
surface due to the relative motion between the blank and the cutting tool. In order to 
unify different fy^pes of kinematic structures of machine tools, it has been proposed to 
model the relative motion using a twelve degrees of freedom model of a machine tool. 
The symbolic algorithm of machining enables the manufacturing engineer to obtain the 
biparametric equation of the machined surface (Symbolics, Inc, 1988; Wallis et al., 1989). 
One can visualize and evaluate computationally the accuracy of this surface before 
finalizing the manufacturing process. 

In the present work, unified models for the geometry of a generic cutter as well 
as for the kinematic structure of a generic machine tool have been proposed using the 
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Figure 1-1 . Model of Shape Synthesis Tools Based on Conjugate and Intrinsic Geometry 
Formulation. 


concepts of envelope theory and conjugate geometry. The three basic determinants of a 
manufacturing process as identified in this model are geometry of the cutter, geometry of 
the machined surface and description of the relative motion between these two surfaces. 
The twelve degrees of freedom model of a manufacturing process .proposed here can be 
used to find any one of the three determinants of a manufacturing process if the other two 
are known. Such a unified model will be useful for cutter path simulation, cutter path 
generation and selection of an appropriate combination of machine tool and cutting tool 
to produce a given surface. It can also be used in the design of cams, gears, extruder 
screws etc. which involve a pair of enveloping surfaces. 


1.2 Literature Review 


Recent developments in the area of geometric modeling in general, and solid 
modeling in particular, have provided design and manufacturing engineers with an elegant 
and precise definition of a component or a product that needs to be designed (Mortenson, 
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1985; Mantyla, 1988; Voelcker and Requicha, 1977; Faux and Pratt, 1983). Based on 
these definitions of a component, several research workers have developed methodologies 
to generate process plans for manufacturing the designed components (Chang and Wysk, 
1985; Kanumury and Chang, 1991; Joshi and Chang, 1988; Pemg et al., 1990; Shah, 
1991; Srinivasan and Ferreira, 1990). The process plans are generally based on the 
geometric features of a component on one hand and the capability of a manufacturing 
process to produce those features on the other hand. Such process plans do not generally 
take into account a micro-level geometric model of the tool and the blank along with the 
kinematic structure of the machine tool used for the manufacturing process. Researchers 
in the area of manufacturing processes have studied the physics and the mechanics of 
several processes (Acherkan et al., 1973; Ghosh and Mallik, 1981; Ber et al., 1988; Sen 
and Bhattacharyya, 1969; Pande and Shah, 1980). However, it has been found that except 
in the area of gear manufacturing, there is very little in the literature to model either 
symbolically or computationally the geometry of several manufacturing processes. 

Simulation of machining processes is one of the major areas of research in 
manufacturing. In order to model a machining process, one needs to take into accoimt 
several factors such as the geometry of the cutting tool, the kinematic structure of the 
machine tool, the cutting speed and the feed rate, the material properties of the cutting 
tool and the material being cut, the cutting forces of the machining process, the energy 
required for machining, etc. Srinivasan and Ferreira (1990) have proposed two levels for 
modeling a machining process. The lower level takes into account the kinematic chain of 
the machine tool and the upper level takes into account each component of the cutting and 
feed motions. Such a classification is helpful for feature definition and process selection 
activity of Computer Aided Process Planning (CAPP). 

The problem of defining the geometry of a machined surface is essentially a 
problem of defining the boundary of the swept volume of the cutter. Wang and Wang 
(1986a, 1986b) have defined the swept volume using the envelope theory approach. They 
however do not model explicitly the kinematic structure of the machine tool. Yap (1988) 
has described some algorithms for computing three dimensional tool paths for planar and 
non-planar machined surfaces. Dhande et al. (1990, 1992) have shown how the geometry 
of a machined surface can be derived for a non-traditional manufacturing process such as 
wire EDM. The discrete simulation of machining sculptured surfaces has been studied by 
Drysdale and Jerard (1989). An approach to define the swept volume using the sweep 
differential equations has been proposed by Blackmore and Leu (1990). All these works 
do not take into accoxmt explicitly the kinematic structure of the machine tool. 

The problem of conjugate geometry between a pair of moving as well as 
contacting surfaces is termed as a higher pair contact problem in kinematics (Boltyanskii, 
1964; Dudley andDarle, 1962; Chakraborty and Dhande, 1977). Such problems have been 
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analyzed by designers for design of cams and gears. Dhande and Chakraborty (1977) have 
described the methodology of defining the kinematics of relative motion as 'well as the 
geometry of the contacting surfaces for three dimensional cams. Several researchers have 
studied the problem of conjugate geometry for planar and spatial gear mechanisms 
(Dudley, 1962; Dhande and Chakraborty, 1977). One can see that there is some amount 
of similarity between the conjugate geometry problems encountered in machining 
processes on one hand and for gear and cam design on the other. 

One of the major computational difficulties in formulating and analyzing the 
problems of conjugate geometry is that for every given situation of a pair of contacting 
surfaces and their relative motion, one has to derive the biparametric equations of the 
surfaces every time. This can be improved by using the symbolic manipulation approach 
for the problems of conjugate geometry especially for modeling different manufacturing 
processes. Reinholtz and Dhande et al.(1990) have described an approach of how 
symbolic manipulation can be used for defining the machined surfaces. They have 
considered only a three degrees of freedom model to show how helically swept surfaces 
can be manufactured using alternative machining strategies. A model with more degrees 
of freedom such as the twelve degrees of freedom model proposed in this thesis, it is 
believed, can be used for instantiating a variety of machine tool structures. The proposed 
symbolic model of the generic cutter along with this model of the machine tool assembly 
provides a unified approach for modeling a variety of manufacturing situations. 


1.3 Envelope Theory 

A brief summary of the important results of the envelope theory are described in 
the present section. A family of curves can be specified in an implicit form by the 
equation 

fiXy y, a) = 0 ^ a ^ (1-1) 

where x andy are the Cartesian coordinates of a generic point and a is a parameter (Figure 
1-2). By giving a specific numerical value of a, say a„ one can obtain the equation of a 
member of the family as 

Six, y, a) = 0 (1-2) 

It is assumed that every member of the family is a C’-continuous curve. 
Furthermore, it is assumed that successive members of the family intersect at a finite 
point. Each of these points of intersections, it is assumed, has a limiting position when 
the successive member curves are considered infinitesimally close and approaching a limit 
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Y 



Figure 1-2. Envelope of a Family of Curves. 


(Boltyanskii, 1964). The point of intersection of a pair of infinitesimally close member 
curves can be obtained by solving the equations 

fix, y, ttj) = 0 

fix, y, a,+Ao) = 0 

Aa being an infinitesimal increase in a. The limiting position of this point of intersection 
can be obtained by eliminating a from the following equations; 

fix, T, «) = 0 
M , 0 

da 


All such limiting positions, in turn, define the envelope to the family of curves. 
In parametric form, this equation of a family of curves is given by 
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/(«, a) = 0 

and a generic point P of this family can be expressed as 

p(u, a) = a) y(u, a)] (1-6) 

where <u < and < a < In parametric form, the limiting positions of the 
points of intersections can be obtained by considering the collinearity of the vectors 

du 


and ^ . In other words, the equation of the envelope in the parametric form can be 

da 

between the following two equations 


(1-7) 


interpretation of Equation (1-7). If the 
will be the velocity and normal vectors at 
any point. When a curve is being swept, every point on the curve is moving along the 
direction of its velocity. If this velocity vector is orthogonal to the normal vector at any 
point on the curve, then that point (of the curve) is the point defining the envelope curve 
at that instant. All other points of the curve will lie within the swept area of the curve. 

In many cases, it is required to find the envelope of a curve or surface in another 
moving frame of reference rather than a fixed frame of reference. Furthermore, the 
parametric form allows us to describe the envelope of either a curve or a surface in a 
symbolic manner. Taking these two aspects into account, a general model of two mutually 
enveloping surfaces has been developed which is considered as a conjugate geometry 
model (Chapter 2). 


obtained by eliminating the parameter u 
(Chakraborty and Dhande, 1977); 

/(«, a) 
n . V 


where 


V = 




and n . 


dUj 


= -1 


It will be useful to give a physical 
parameter a represents time then Y ^nd n 
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1.4 Symbolic Computation 

Computers have been used traditionally to solve problems which could be 
expressed numerically, i.e., all symbols used in the problem are bound to numerical values 
before they are referred. On the other hand, some problems can be expressed best in 
symbolic terms, or perhaps can be expressed only that way. Solving such problems where 
the variables are manipulated even when they do not have numerical bindings is called 
symbolic computation. It is also known by many other names such as symbolic algebra, 
symbolic mathematics, symbolic manipulation, analytic calculation, computer algebra, 
algebraic manipulation, formula manipulation etc. 

The idea that computers, or calculating engines as they were known, could perform 
algebraic rather than numerical calculations was first suggested nearly 1 50 years ago by 
Lady Ada Lovelace'. However, only after 1960, the hardware and software capabilities 
became reasonably adequate to permit the development of useful computer algebra 
systems. Since this time, the systems developed for symbolic computation can be divided 
into three notional groups, viz.. Group I, Group II and Group III (Table 1-1). Initially the 
packages were written to solve problems of specific areas such as general relativity, 
quantum theory, celestial mechanics and high energy physics. Since the usefulness of 
these packages for other areas also was realized, general purpose symbolic algebra 
systems were developed. They were very general and hence they required very huge 
memory and they also were very sluggish. The early symbolic algebra systems were not 
interactive. To-day, the speed of computers have gone up by virtue of the latest 
technologies like RISC and parallel processing. The symbolic packages such as MathCAD 
and MACSYMA are available for a variety of platforms including personal computers. 
They are highly user-friendly as they make use of sophisticated user-interfaces like 
Windows. In view of this, ^mbolic computing is competing with even numerical 
computing. 

The symbolic manipulation systems can perform a variety of operations. They can 
perform matrix algebra and trigonometric calculations, solve a system of equations and 
find roots of an expression. They can also solve problems involving differentiation and 
integration. Some of the packages like MACSYMA can also find the numerical solutions 
wherever required. It can evaluate a function as one of its arguments approaches a 
specified limit. The results of the calculations are output in more readable multi-line 
format. These packages also have graphic capabilities with automatic scaling etc. for 
depicting the results graphically. Most of the symbolic packages can generate output in 


1 


Lady Ada Lovelace was patron of Charles Babbage who is usually credited with the development 
of World's first computer. 





Table 1-1. Classification of Symbolic Computation Systems 
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the format acceptable to other systems like FORTRAN and C compilers. 


The advantages of symbolic computing can be listed as given below. 


a) The purpose of computing is insight and not numbers. If we write 


X 







= 1 , it immediately occurs to a mathematician that it is an 


ellipse. This is a symbolic expression. On the other hand, if the coordinates of 
a set of points lying on this ellipse are given, one cannot visualize it as an 
ellipse unless he plots these points and performs a curve-fitting. Thus 
computing symbolically helps fulfill the purpose of providing analytical rather 
than numerical answers. In other words, an algebraic solution is usually more 
succinct and meaningful than a numerical solution. 


b) When a computation is performed symbolically, the algorithm maintains the 
generality of computation. Hence the symbolic ^stem is inherently well-suited 
for developing generic algorithms. Dhande (1976) has derived the expressions 
for the conjugate surfaces of more than 20 different cases of planar and spacial 
cam mechanisms. All those computations can be solved by a single symbolic 
algorithm. In fact, this was die motivation for the selection of the symbolic 
approach for implementing the models of the generic machine tool and the 
generic cutter. 


c) Numerical systems operate with floating point approximations whose precision 
is limited by the hardware whereas the algebraic solutions are always exact. 
Rounding-off can be delayed till the last step to minimize its influence on the 
final result when symbolic computation is used. The symbolic algorithms can 
deal with zero and infinity more elegantly. 

d) A symbolic algorithm can solve an analytical problem in a few minutes which 
could be taxing a capable mathematician for several months. One wonders how 
many tractable problems remain unsolved or have been forgotten about simply 
because they were making excessive demands on a researcher's time and 
sanity. Moreover, the solution will be free from errors caused due to fatigue. 

e) Symbolic algebra minimizes the need for tables of solutions to specific 
problems. Indeed, one of the triumphs of these systems is that they have 
highlighted several errors in tables of integrals used by engineers and 
scientists. Hence, one can concentrate on the problem formulation rather than 
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remembering various numeric techniques for solving the problem. 

f) For some problems where the analytical solution may not exist, it is preferable 
to simplify the expressions using the symbolic system to produce more 
efficient and accurate results. 

Mani (1989) has demonstrated the use of a hybrid system consisting of the 
symbolic manipulation system MACSYMA and the numerical compiler FORTRAN for the 
analysis of mechanical systems such as a robot arm. As the speed of computers increases, 
it is expected that more and more applications will make use of the symbolic approach 
which were hitherto solved numerically. 


1.5 Organization of the Work 

The objective of the present work is to develop the mathematical formulation and 
the symbolic algorithms required for the geometric modeling of manufacturing processes 
in a unified way. The methodology proposed here can be used in developing a user- 
friendly menu-driven software that a manufacturing engineer can use effectively. This 
report consists of eight chapters. In Chapter 1, brief descriptions of envelope theory and 
symbolic computation are presented following an extensive literature review. 

Chapters 2, 3 and 4 deal with the conjugate geometry formulation. In Chapter 2, 
mathematical interpretation of Acherkan's four methods of obtaining the generatrix and 
directrix of a manufacturing process is given (Section 2.2) followed by several examples 
of manufacturing processes identifying the generatrix, directrix and the degrees of 
freedom involved in each case (Section 2.3). After this, the conjugate geometry 
formulation and its symbolic algorithm are described in Section 2.4. The instantiation of 
the generic model of the machine tool for a variety of machine tools is presented in 
Chapter 3. Also given is the illustration of the model for machining an extruder screw 
surface. A generic model of the cutter and its instances are described in Chapter 4. 

In Chapter 5, the usefulness of the conjugate geometry model of a machine tool 
and the model of the generic cutter have been demonstrated by means of examples such 
as gear bobbing and shaping. Alternative machining strategies of a helicoidal and 
hyperbolic paraboloidal surfaces have also been modeled using the proposed models. It 
is concluded with an example on design and manufacture of a cylindrical cam and 
follower mechanism. 

The use of the proposed methodology for the design and manufacture of extruder 
screws has been investigated in Chapters 6 and 7. In Chapter 6, the procedure for the 
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design of a pair of extruder screws of different types is given. Chapter 7 is for the 
manufacture of extruder screws. Procedures for the design of cutter to machine a given 
screw surface and for the determination of the screw surface for a given cutter profile are 
given in this chapter. 

Finally, the conclusions and future extensions of this work are presented in 
Chapter 8. 


Chapter 2 


CONJUGATE GEOMETRY 
MODEL OF A MACHINING PROCESS 


2.1 Introduction 

The proposed methodology of a generic machine tool model is based on the 
premise that all manufacturing processes are essentially shape realization processes. 
Furthermore, all such shape realization processes can be classified into two groups, viz., 
(i) shape generation processes and (ii) shape formation or shape transformation processes. 
The shape generation processes are those processes which involve realization of a surface 
by removal of material. Generally, such material removal requires movement of a cutter 
through the volume of the blank. Processes such as milling, bobbing and shaping can be 
classified as shape generation processes (Figiue 2- la). Another class of shape realization 
processes involves transformation of a given volume from one state to another state. Such 
transformations generally involve plastic deformation of the material. Processes such as 
deep drawing, forging etc. can be cited as examples of this category (Figure 2-lb). 

To begin with, a shape can be considered as the mathematical description of the 
extrinsic and intrinsic geometry of the boundary of a closed regular set. Let iS be a closed 
regular set of points in either a 2-D or a 3-D space. Let b{S) be the set of points 
belonging to a closed regular set and defining the boundary of that set. The extrinsic 
properties are the properties such as the Cartesian coordinates of the points on the 
boimdary of a set while the intrinsic properties are the properties such as the arc length 
and the curvature at a point on the boundary of that set. 

It should be kept in mind that a mathematical description of the shape can be 
specified depending on die purpose of its use. For instance, if a shape is to be displayed 
on a computer-controlled display device then a parametric definition is appropriate 
(Figures 2-2b & 2-3b). If a shape is to be specified for some analysis procedure then an 
explicit or an implicit closed-form expression in Cartesian coordinates is useful (Figures 
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2-2a & 2-3a). However, if a shape is to be realized by a generation process then neither 
of these definitions are useful. 




(a) Example of a Shape Generation 
Process 


(b) Example of a Shape Transformation 
Process 


Figure 2-1. Shape Generation and Shape Transformation Processes for Shape Realization. 


Consider a case of generating a circle from a 2-D lamina. In order to realize the 
circular shape, one way could be to employ a series of straight line cuts (for example, 
using a pair of scissors). These straight line cuts constitute a family of curves (Figure 2- 
2c). An envelope to this family is the shape generated. In this case, it happens to be a 
circle. Similarly, if one needs to display a straight line dien it can be defined in a 
parametric form (Figure 2-3b). However, if one needs to realize a straight line shape then 
it can be generated as an envelope of a family of circles. The family is generated by 
moving the center of the circle on a linear path (Figure 2-3 c). Interestingly, it can be 
seen that in the examples shown in Figure 2-2 and Figure 2-3 a family of straight lines 
is generating a circular shape wdiile a family of circles is generating a linear shape. This 
observation is important for a manufacturing engineer. It shows that a shape specified by 
a designer is realized by a generation scheme. This generation scheme employs a shape 
that is not necessarily the same but a different one, called as the conjugate shape, for the 
realization of the designed shape. It is also important to note that the manner in which the 
family of curves or shapes is defined is as important as the definition of the conjugate 
shape for designing a shape generating process. Generally, the mariner in which the family 
of curves or surfaces of the cutting tool is defined is determined by the kinematic 
structure of the machine tool. 


Traditionally, the envelope theory formulation finds the definition of an enveloping 
curve or a surface in a fixed frame of reference (Boltyanskii, 1964; Wang and Wang, 
1986). This is not as useful for modeling several manufacturing processes. In many 
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manufacturing processes, both the cutter as well as the blank may be moved in order to 
realize a surface. In such cases, it is necessary to find the envelope of the tool geometry 
(curve or surface) in another moving frame of reference attached to the blank body. The 
kinematics of the machine tool determines the relative motion between the tool and the 
blank. Such an extension of the envelope theory formulation can be termed as the 
conjugate geometry formulation. The conjugate geometry formulation defines the envelope 
in one moving frame of reference due to a family generated by another shape (curve or 
surface) and due to the relative motion between these two shapes, i.e., the cutter and the 
blank. This conjugate geometry model which is an extension of the envelope theory is the 
mathematical basis for modeling the geometry of several shape generating processes. 


2.2 Shaping of Geometrical and Real Surfaces 

Acherkan et al. (1973) have classified surfaces as geometrical and real surfaces. 
A geometrical surface is the one conceived by a designer. It can be defined as the trace 
or the sweep obtained due to the motion of one geometrical generating line, the 
generatrix, with respect to another geometrical generating line, the directrix. Alternatively, 
a geometrical surface can also be thought of as the boundary of a set resulting from the 
Minkowski addition of a set of points belonging to the generatrix and the set of points 
belonging to the directrix. If A and B are two closed regular sets of geometric points then 
the Minkowski addition of A and B is defined as a set-theoretic union of all the mapping 
of the set A when it is copied at every element of B (Ghosh and Mudur, 1984; Kaul, 
1992). In the present case, the set A is the generatrix curve and the set B is the directrix 
curve. This definition of a geometrical surface is inadequate for defining surfaces realized 
in several manufacturing processes. For this purpose, it is necessary to consider the 
definition of real surfaces. 

Real surfaces are obtained on a solid body by any one of the several methods of 
metal working. All these methods produce a real surface by shaping or removing material 
with the aid of auxiliary material elements which could be either points, lines (i.e. straight 
or curved) or real surfaces. The designed geometrical surfaces are produced due to the 
relative motions of these real auxiliary material elements. These relative motions are 
called formative motions. In practice the shaping of a material surface involves realization 
of the geometrical generating line through the auxiliary material element. The designed 
surface is produced due to the formative motions of the geometrical generating lines. 

There are four methods of realizing geometrical generating lines from the auxiliary 
material elements which may be in the form of a line or a point. These are (i) forming 
method, (ii) generating method, (iii) tracing method and (iv) tangent method. 



Chapter 2. Conjugate Geometry Model of a Machining Process 


18 



Auxiliory moteriol line 



Produced line 


Auxiliary materia! point 
r Produced line 


rriuuuv.cu nil 


•Auxiiiory materiol point 
\ 'cAuxiliory geometrical circles 


/ Auxiliary 

Produced line geometricol 

straight line 


Figure 2-4. Methods of Producing Geometrical Surfaces 

(a) Forming Method 

(b) Generating Method 

(c) Tracing Method and 

(d) Tangent Method 





Chapter 2. Conjugate Geometry Model of a Machining Process 


19 


2.2.1 Forming Method (Figure 2-4a) 

In this method, the geometrical generating line is produced without any formative 
motion since it is a copy of the auxiliary material line. The geometrical generating line 
coincides with the shape and extent of the auxiliary material line. The geometrical 
generating line can be defined as, 

fix, y) = 0 (2-1) 

where / (x, y) = 0 is the geometry of the auxiliary line. 


2.2.2 Generating Method (Figure 2-4b) 

In this method, the geometrical generating line is obtained as an envelope of the 
auxiliary material line when the latter rolls along the line being produced. As such, only 
one formative motion is involved in the process. If /, is the parameter of the rolling 
motion then the geometrical generating line is an envelope of a family of curves due to 
the motion The equation of the enveloping geometrical generating line is obtained by 
eliminating the parameter r'l from the following equations, 


fix, y, ij) = 0 



( 2 - 2 ) 


2.2.3 Tracing Method (Figure 2-4c) 

The auxiliary element in this method is in the form of a material point. The 
geometrical generating line is found as the trace of the material point in its motion. Since 
a line or a curve can be described as a uni-variate geometrical entity, it can be said that 
the generating line is produced with a single formative motion. In practice, the material 
point will have a short length of the tool. However, it is assumed to be sufficiently small 
so as to idealize it as a point. If is the formative motion then their geometrical 
generating line can be described as. 


X = xii^ 

y = ^(y 


(2-3) 
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Formative 

Type of 

Auxiliary 

Element 

Method of Obtaining 
Geometrical Generating Lines 

Number of 
Formative 
Motions 

Processing Method 

Tool 

Generatrix 

Directrix 

Casting 

Foundry 

moulds 

Surface 

Forming 

Forming 

0 

Die forging 

Dies 

Surface 

Forming 

Forming 

0 

Blanking 

Dies 

Line 

Forming 

Tracing 

1 

Knurling and rolling 

Knurls or 
rolls 

Line 

Forming 

Generating 

1 

Ball burnishing 

Ball 

Surface 

Generating 

Generating 

2 

Metal cutting 
(machining) 

Cutting 

Edge 

Line or 
point 

Tangent 

Generating 
tracing or 
tangent 

1 to 3 

Grinding 

Grinding 

wheel 

Line or 
point 

T angent 

Genera ring 
tracing or 
tangent 

2 or 3 

Honing 

Abrasive 

sticks 

Point 

Tracing 

Tracing 

2 

Abrasive vapour 
blasting 

Abrasive 

slurry 

Surface 

Forming 

Tracing 

1 

Buffing 

Buffing 

wheel 

Point 

Tangent 

Tangent 

3 

Lappi ng 

Lap~ 

Point 

Tracing 

Tracing 

2 

Super- finishing 

Abrasive 

stones 

Point 

Tangent 

Tangent or 
tracing 

3 

Ball hole- burnishing 

Ball 

Line 

Forming 

Tracing 

1 

Electric discharge 
machining 

Electrode 

Line 

Forming 

Tracing 

t 

■plastic moulding 

Moulds 

Surface 

Forming 

forming 

-a 

Metal 1 izing 

Work piece 

Surface 

Forming 

Forming 

0 

Chromium and nickel 
electroplating to size 

Work piece 

Surface 

Forming 

Forming 

0 

Ultrasonic machining 

Sonic 

contour 

Line 

Forming 

Tracing 

1 

Wire-drawing 

Drawing 

die 

Line 

Forming 

Tracing 

1 


Table 2-1 Interrelation between the Processing Method and the Method of Geometrically 


Shaping a Surface. 
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2.2.4 Tangent Method (Figure 2-4d) 

In this method, the geometrical generating line is tangent to a series of 
supplementary auxiliary lines produced by the material point either by the tracing method 
(Figure 2-4c) or by the tangent method (Figure 2-4d). In either case, this method requires 
two or more formative motions. 

In Figure 2-4d, the auxiliary material point produces the auxiliary geometrical 
circle by a rotational formative motion described by the parameter u by the tracing 
method. The family of the auxiliary geometrical circles produces the auxiliary' geometrical 
straight lines by an oscillatory formative motion described by the parameter which is 
the generating method. The family of the auxiliary geometrical straight lines generates a 
circle by yet another rotational formative motion described by the parameter by the 
generating method. Thus, as shown in Figure 2-4d, the generating line is produced by one 
tracing method and two generating methods. If the equation of the material point at any 
instant is represented as. 


fix, y, u, iy = 0 (2-4) 

then the produced line (shown as dotted line in Figure 2-4d) is obtained by eliminating 
the parameters u, and ^4 from the following equations: 

fix, y, u, * 3 , 14 ) = 0 

3 m 3(3 3I4 

2.3 Shaping of Real Surfaces 

The foregoing discussion shows that a real surface can be obtained by defining a 
geometrical line for the generatrix as well as another one for the directrix of a surface. 
Since a geometrical line can be defined by any one of these four methods, one can 
conceive sixteen different conceptual schemes for shaping a real surface. These have been 
described by Acherkan et al. (1973). Schematic diagram of this classification is given in 
Table 2-1. 

As shown in Figure 2-5, manufacturing processes such as casting and die-forging 
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belong to the category where the generatrix as well as the directrix are realized by 
forming methods. It should be noted, however, that one auxiliary motion is required for 
setting the material curve or surface in its final position. Furthermore, there is no 
manufacturing process where the generatrix is realized by either the generating, tracing 
or tangent method and the directrix is realized by the forming method. Hence, in the 
present work these three cases will not be considered further. 

Manufacturing processes such as thread rolling and knurling belong to the category 
where the generatrix is realized by the forming method and the directrix is realized by the 
generating method. From Figure 2-6, one can see that the thread profile is the generatrix 
and the rolling and axial motion of the roll is the directrix. It should be noted that this 
manufacturing process has only one formative motion or the degree of freedom since the 
roll and the blank have a relative motion of pure rolling. 

Figure 2-7 shows the manufacturing process of ball burnishing. In this method, 
both the generatrix and the directrix are realized by the generating method. The ball has 
two formative motions and the resulting enveloping surface is the burnished contour 
surface. 

The manufacturing of spur gears using a circular shaping cutter is shown in Figure 
2-8. In this class of processes, the generatrix is realized by the tracing method and the 
directrix is realized by the generating method. The linear motion of the cutter is the 
formative motion for realizing the generatrix and the planetary motion of the cutter is the 
formative motion for realizing the directrix. 

If a disk type grinding wheel is used for grinding the involute shaped cylindrical 
surface of a spur gear then the process can be classified in a class where the generatrix 
is realized by the tangent method and the directrix is realized by the generating method. 
As shown in Figure 2-9, this manufacturing process has three formative motions. A 
variety of grinding processes belong to the category where the generatrix is realized by 
the tangent method and the directrix is realized either by the generating or the tracing or 
the tangent method. 

The case of cutting a thread on a cylindrical blank with a form tool using the 
screw cutting lathe belongs to a category where the generatix is realized by the forming 
method and the directrix is realized by the tracing method. As shown in Figure 2-10, such 
a method requires only one formative motion. A large number of conventional and 
advanced machining processes belong to this category. Processes such as ultrasonic 
machining, electric discharge machining, wire drawing, blanking and ball hole-bumishing 
belong to this category. In almost all these cases, one formative motion is required. 
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A relatively less frequently occurring category is the situation where the generatrix 
is realized by the generating process and the directrix is realized by the tracing process. 
The familiar process of gear bobbing belongs to this category (Figure 2-11). In this case, 
the rotations of the blank and the hob produce the generatrix. Here, the generatrix is the 
involute curve and the directrix is realized due to the linear motion of the hob along the 
axis of the blank. In this class of processes, the total number of formative motions 
required is two. It should also be noted that in this process the auxiliary material line 
responsible for generating the involute profile is a straight line cutting edge of the hob. 

Manufacturing processes such as honing and surface relieving belong to the 
category where the generatrix as well as directrix are realized by the tracing method. In 
case of generating the relief surfaces on the helical teeth of a plain milling cutter, a 
pointed tool traces a curve in the space due to two formative motions (Figure 2-12). The 
first motion generates an Archimedean spiral and the second motion is responsible for 
generating helix. In this class of manufacturing processes, the auxiliary material element 
can be idealized as a point or a short line representing the cutting edge of the tool. 

Manufacturing processes such as super-finishing using abrasive stones, some 
categories of grinding and some processes of metal cutting belong to the category where 
the generatrix is realized by the tangent method and the directrix is realized by the tracing 
method. Figure 2-13 shows the case of grinding a generalized cylindrical surface using 
a disk type grinding wheel. In this case, three formative motions are required. 

When a thread is milled with a single thread cutter, the generatrix is realized using 
the forming method and the directrix is realized by the tangent method (Figure 2-14). In 
this case, the directrix realization requires two formative motions. The realization of the 
generatrix requires one formative motion. However, the formative motion for realization 
of generatrix coincides with one of the formative motions required for the realization of 
the directrix. Since, these two motions are coinciding, the net number of the formative 
motions required for this process is two. 

Another class of manufacturing processes which is rarely encountered in practice 
is the category where the generatrix is realized by the generating method and the directrix 
is realized by the tangent method. The case of grinding the flank surfaces of a spur gear 
using a helically profiled worm-type grinding wheel belongs to this category. For this 
process, one would normally require three formative motions, one for the grinding wheel 
rotation, one for the longitudinal feed and one for the generation of involute. However, 
the grinding wheel rotation cannot be independent of the motion required for involute 
generation (Figure 2-15). Hence, the resulting number of formative motions required in 
this case is two. 
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Figure 2-5 A Manufacturing Process where both the Generatrix and the Directrix are 
Produced by the Forming Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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Figure 2-6 A Manufacturing Process where the Generatrix is Produced by the Forming 
Method and the Directrix is Produced by the Generating Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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Figure 2-7 A Manufacturing Process where the Generatrix is Produced by the Generating 
Method and the Directrix is Produced by the Generating Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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Figure 2-8 A Manufacturing Process where the Generatrix is Produced by the Tradng 
Method and the Directrix is Produced by the Generating Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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Figure 2-9 A Manufacturing Process where the Generatrix is Produced by the Tangent 
Method and the Drectrix is Produced by the Generating Method. 

(a) Schematic .Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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Figure 2-10 A Manufacturing Process where the Generatrix is Produced by the Forming 
Method and the Directrix is Produced by the Tracing Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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Figure 2-12 A Manufacturing Process where the Generatrix is Produced by the Tracing 
Method and the Directrix is Produced by the Tracing Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 



























Chapter 2. Conjugate Geometry Model of a Machi ning Process 


Grinding wheel 




Figure 2-13 A Manufacturing Process where the Generatrix is Produced by the Tangent 
Method and the Directrix is Produced by the Tracing Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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Figure 2-14 A Manufacturing Process where the Generatrix is Produced by the Forming 
Method and the Directrix is Produced by the Tangent Method. 

(a) Schematic Diagram of the Process 

(р) Coordinate Systems for the Conjugate Geometry Model and 

(с) Table of the Degrees of Freedom. 
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Figure 2-15 A Manufacturing Process where the Generatrix is Produced by the Generating 
Method and the Directrix is Produced by the Tangent Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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Figure 2-16 A Manufacturing Process where the Generatrix is Produced by the Tangent 
Method and the Directrix is Produced by the Tangent Method. 

(a) Schematic Diagram of the Process 

(b) Coordinate Systems for the Conjugate Geometry Model and 

(c) Table of the Degrees of Freedom. 
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As indicated by Acherkan et al. (1973), the category of shaping of surfaces where 
the generatrix is realized by the tracing method and the directrix is realized by the tangent 
method has found no application as yet in its true form, hence, this category will not be 
dealt with any further in our present analysis. 

The last class of manufacturing processes belongs to the category where the 
generatrix as well as the directrix are realized by the tangent method. Some of the 
machining processes belong to this category which is also referred to as the double 
tangent method. The case of milling a doubly curved surface with an end mill using a 
tracer controlled duplicating machine belongs to this category (Figure 2-16). In general, 
the double tangent method requires four formative motions, two for realizing the 
generatrix family of curves and two for realizing the directrix family of curves. How'ever, 
the cutting motion is common in both the cases and at any given time, only one curve is 
generated. Hence, the formative motions required are two at any given time. The 
manufacturing process of buffing also belongs to this category of double tangent method. 

The foregoing review of the manufacturing processes on one side and their 
characterization as one of the sixteen shape realization processes on the other side shows 
that different schemes of kinematic structures and tool geometries are employed to realize 
simple and complex surfaces. 


2.4 Definition of the Conjugate Geometry Model of a Generic Machine Tool 

There are basically three elements involved in any machining process. The first 
one is the surface of the cutting tool denoted as Z, and the second is the surface that is 
to be produced denoted as (Figure 2-17). The third element is the relative motion 
between these two surfaces which will be constrained by the t 3 q)e of kinematic chain of 
the machine tool. To produce the surface Z 2 correctly, it is essential that the two surfaces 
X] and Z 2 are conjugate to each other. In other words, these two surfaces will maintain 
a higher pair contact throughout die cutting motion. Using this principle of conjugate 
geometry and envelope theory, a unified model of the machine tool with twelve degrees 
of freedom has been proposed in the present thesis. Among the three elements of a 
machining process, viz., the cutter surface L,, the surface to be machined In and their 
relative motion, if any two are known, the third one can be determined using this model 
of the generic machine tool. Such a model can be used for the determination of the NC 
cutter path wherein the cutter surface X, and the surface to be machined In known and 
their relative motion is required to be determined. This model can also be used for NC 
path simulation which is required for off-line NC tape proving, error analysis of the 
produced surface, NC code optimization etc. The former will fall in the category of the 
reverse kinematic problem and the latter is a forward kinematic problem. Also one can 
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use the proposed conjugate geometry model for choosing the optima! combination of the 
machine tool and the cutting tool by an iterative process. 

A generic model of a machine tool is essentially a representation of the kinematic 
structure of the machine tool (Acherkan, 1973). A kinematic structure of a machine tool 
describes the kinematic chain necessary for providing the cutting motion to the cutter and 
the feed motions to the cutter and/br the blank. In this model, the cutter and the blank are 
considered to be two free bodies in space and hence their position and orientation can be 
described by six degrees of freedom each. These twelve degrees of freedom are expressed 
with respect to a fixed frame of reference {0^ - Xq, Yq, Zq). A coordinate system S. 
(O, - X 2 , ¥ 2 , Z,) is considered to be attached to the blank and another coordinate system 
(Oj - X^, y,, Z,) to the cutter (Figure 2.18). The initial positions of the coordinate 
systems and 5. are denoted as (<9f - X^, 7,-, Zf) and S 2 (O^ - AV, Tj, Z5) (Figure 
2-19). Let (<3,, bj, c,) and { 02 , ^2, Cj) denote the translation of the origins O, and Oj of the 
moving coordinate frames and let (6j, (|)„ \|/,) and (6^, (j). and Xfr,) denote the rotations of 
the coordinate frames and .S', respectively. As these twelve parameters of motion, viz., 
a,, bi, Cj, 0j, (j);, \|rj, Oj, ^>2, C 2 , 62, (J); and \|/2 as shown in Figure 2-18 represent the 
kinematic arrangement of the machine tool, they can be expressed as functions of time 
t. Since we assume that the cutter is rigid and the joints and links of the machine tool are 
perfect (i.e., no play and no deflection), the tool surface Z, and the surface produced on 
the blank Z2 will be time-invariants. The known surface among the two surfaces can be 
expressed as a biparametric surface in u and v. Since the two surfaces Zj and Zj are in 
higher pair contact, at the point of contact P, 

V[!^] = ^p[°2M] (2-6) 

and 


where. 



Point P with reference to the cutter coordinate system 5,. 

Point P with reference to the blank coordinate system S 2 . 

Unit normal vector of the cutter surface Z, at point P with reference to 
coordinate system S^. 
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Unit normal vector of the blank surface Zj point P with reference to 
coordinate system S 2 . 

Homogeneous transformation matrix for transforming a vector from the 
coordinate system S', to . 

Homogeneous transformation matrix for transforming a vector from the 
coordinate system Sj to Sq. 

It is to be noted that all vectors such as position vectors, normal vectors and velocity 
vectors are expressed as row vectors in the homogeneous coordinate s>'stem. Hence the 
fourth element of any position vector is unity and it is zero for the normal and velocity 
vectors. 

The surface Z, is the boundary of the swept volume of the cutter surface Z, which 
is nothing but the envelope to the various instances of the cutter surface. In order to 
identify the point on the cutter surface Z, which is cutting the blank and producing (or 
generating) the surface Zj at any given time, it is necessary to establish some additional 
constraints. This constraint is called as the condition of cutting or the condition of contact 
or the condition of envelope. It can be said that in order to ensure a smooth contact or 
cutting at any given time, the relative velocity vector between the cutter and the blank 
must be orthogonal to the common normal at the point of cutting (Chakraborty and 
Dhande, 1977). Mathematically, this condition can be expressed as, 

^ (2-8) 

and 

'""Xm = 0 (2-9) 



where, 

PI - Unit normal vector of the cutter surface Z, at point P with reference to 
coordinate system S^. 

. Unit normal vector of the cutter surface Zj at point P with reference to 
coordinate system S^. 

- Relative velocity vector between the tool and blank with reference to 
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coordinate system Sq. 

The condition of contact is a scalar equation relating parameters u, v and t. Thus, 
at any given instant of time t, the condition of contact reduces to a relation between u and 
V. For any value of u within its feasible range, one can find the corresponding value of 
the parameter v. Thus, one can locate a set of points or a curve on the surface Xj which 
defines the surface curve for Xj as well. The corresponding points of the surface X2 can 
be obtained using the equation, 

[Jm] (2-10) 

A symbolic algorithm, that can be programmed in any of the algebraic 
manipulation systems such as MACSYMA, MathCAD or REDUCE can now be stated as 
follows: 

1 . Define the position vector V of a generic point P on the cutter surface X, as a 
biparametric surface in u and v such that 

V («, v) = [ Xiiu, v) y^iu, v) Zi(u, v) 1] (2-11) 

2 . Define the twelve parameters of motion a,, i,, c,, 0 ,, 4 >,, x^i, 02, Cj, 62, (j)2 and 

in terms of time t. 

3 . Calculate the matrices j and [1^] as follows: 


where. 
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cos4)j 0 0 

0 10 0 

sin <l)j 0 cos <j)i 0 

0 0 0 1 


cos fj 

sinilTj 

0 

0 ■ 


■ 1 

0 

0 

0 ■ 

-sin tjfj 

COSljTj 

0 

0 

fjrl - 

0 

1 

0 

0 

0 

0 

1 

0 

J 

0 

0 

1 

0 

0 

0 

0 

1 . 


. 



1 


Matrix j is the time-differential of matrix jAf j . 


■ 1 

0 

0 

0 


0 

cos 6, 

sine, 

0 


0 

1 

-sinOj 

cosGj 

0 

0 

0 

0 

1 . 



4. 


Calculate the matrices j , [^ ] and j as follows; 


(2-13) 


where. 



■ 1 

0 

0 

0 


C06 4>2 

0 

-sin<j)2 

0 ■ 

0 

COS0, 

sine. 

0 

[^] = 

0 

1 

0 

0 

0 

-sin 02 

cos 02 

0 

sinij). 

0 

COS<})2 

0 

0 

0 

0 

1 . 


0 

0 

0 

1 . 



cosilTj smi|r2 0 0 
-sinijr2 cosiir2 0 0 
0 0 10 
0 0 0 1 



1 0 0 0 ' 
0 10 0 
0 0 10 
fl2 ^2 ^2 1 
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(2-14) 


Matrix ^ j is the time-differential of matrix ^ J 


5. Derive the position vector P . 


V = [ ^ 3^2 ^ 1 ] 




(2-15) 


6. Derive the normal vector p._i 






du 


dv 


(2-16) 


where, 

a a V 

and are the partial derivatives of !p with respect to u and v 

du dv “ ■ 

respectively. Therefore, 


av 



^ 0 

av 



^ 0 

du 

du 

du 

aa 

av 

av 

av 

dv 


d^p d^p 


When the vectors 


and 


du dv 

normal vector /»] is defined as 


are expressed in Cartesian coordinates, the 
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i 

j 

k 


dyi 

dz, 

du 

du 

du 


?Il 


dv 

dv 

dv 



d^p 


The same can be redefined when the vectors and 

du 


d^p 


— are expressed in 
dv 


homogeneous coordinates as 


7. 



Calculate the unit normal 


i 

j 

k 

djc, 


dz, 

Tu 

du 

du 



dz, 

dv 

dv 

dv 

0 

0 

0 

vector 

^p,i ■ 





1 

1 


(2-17) 


where I p,i I is the magnitude of 


8. Calculate the unit normal vector “w p,, in the global coordinate system. 


- l»f 


;m] 


(2-18) 
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9. 


Calculate the relative velocity vector ^,2 



(2-19) 


where 



10. Reduce symbolically the condition of contact 

~ ^ ( 2 - 20 ) 

2 

11 . P obtained in step 5 will be a function of u, v and t. Use the condition of 
contact obtained in step 10 to eliminate either of U and v. This will yield the 
surface Sj- 


It has been stated earlier that based on the method of obtaining the generatrix and 
directrix, the manufacturing processes can be classified into sixteen groups out of which 
only twelve are encountered in practice. An example for each of these twelve 
manufacturing processes has been given in Figures 2-5 through 2-16. Part (a) of these 
figures gives the schematic diagram of the process. The fixed global coordinate frame Sq, 
the moving coordinate frames Sj and Sj and their initial positions Sf and Sj are depicted 
in part (b). Part (c) contains a table in which the expressions for the twelve motion 
parameters are given as functions of time t. Out of these twelve motion parameters, it may 
be observed that at the most only six will be independent of each other and the rest of the 
parameters can be expressed as constants or as functions of these independent parameters. 
This is because of the fact that the cutter surface is always in contact with the blank 
siurface. Furthermore, for a milling process, the number of independent motion parameters 
are limited to five since the sixth independent parameter, viz., the spindle rotation, is a 
passive degree of freedom as far as the present conjugate geometry model is concerned 
(In a way, it can be said that the spindle rotation is related to the sweep parameter v in 
the definition of the cutter surface for a milling process - see Chapter 4). This is the 
reason why most of the CNC machines have only upto five axes. The generic model with 
twelve degrees of freedom assiunes the characteristics of a given machine tool when the 
twelve motion parameters pertaining to this machine are defined. 

In the next chapter, the conjugate geometry formulation has been explained with 
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the help of an example. Also the instantiation of the generic model of a machine tool for 
five different types of machine tools is illustrated. 



Chapter 3 


INSTANCES OF 

THE GENERIC MACHINE TOOL 


3.1 Introduction 

The model of a generic machine tool with twelve degrees of freedom was 
proposed in the previous chapter. Even though the proposed model is designed to handle 
twelve degrees of freedom, such a machine with twelve independent degrees of motion 
cannot exist in reality. The fact that the surfaces X, and are in continuous contact 
implies that there could be only six relative degrees of freedom, in other words, 
independent degrees of freedom, at the most And for a milling-like process, one cannot 
think of a machine tool which has more than five independent motion parameters since 
the sixth motion which is the spindle rotation is a passive degree of freedom. It vsnll be 
obvious from the following chapter that the spindle rotation may be used to describe the 
generatrix profile in the biparametric definition of the cutter surface. Even though all 
twelve degrees of freedom will not be used by a single process, by incorporating the 
highest number of degrees of freedom in this model, the greatest amount of versatility is 
achieved in order that any kind of machine tool can be instantiated using this model. The 
need for a unified model of this kind can be justified as follows: 

Other classical methods such as the APT algorithm for cutter path generation 
assume a machine tool which has a fixed table and five degrees of freedom for the cutter. 
The cutter path generated from these systems is a series of tool tip positions along with 
the corresponding direction cosines of the tool axis. The task of converting this cutter path 
called CL data into the joint coordinates (kinematic parameters) of a particular machine 
tool has to be performed by another program called post-processor. The design of the 
post-processor becomes more complicated when it is required to customize the CL data 
for a machine tool which has few degrees of freedom for the table also since it has to 





Chapter 3. Instances of the Generic Machine Tool 


48 


perform interpolation in order to adhere to the LINTOV limits. In the conjugate geometry 
model, the kinematic information of the machine tool is inherently taken care of thus 
alleviating the need for post-processing^. 

The output in the APT paradigm is discrete and numerical whereas the proposed 
conjugate geometry model generates the output surface in terms of the joint coordinates 
in the form of algebraic equations which can be reduced to any desired accuracy. 

Today, most of the NC programming packages come with built-in path simulation 
routines. There are several packages such as VeriCUT which are exclusively for NC 
program verification. But all these simulation programs take the CL file as input. Hence 
even if the CL file appears satisfactory during the simulation, one cannot be sure that the 
same surface would be produced during the actual machining since errors such as wrong 
start-up and retraction and LINTOL interpolation errors could creep in at the post- 
processing stage. The simulation based on the conjugate geometry formulation will be free 
from these problems. 


3.2 Instances of the Generic Machine Tool 

Figure 3-1 to Figure 3-5 illustrate the instantiation of the conjugate geometry 
model of the generic machine tool for five different kinds milling machines. The table 
given in each of these figures is the kinematic data of the corresponding machine tool 
input to the symbolic algorithm. 


3.2.1 Instantiation of a 4 Axis Vertical Milling Machine 

Figure 3-1 shows a four axis CNC vertical milling machine. In this machine, the 
translation along the Xj and Zj directions are effected by the table movement and the 7, 
movement is due to the movement of the spindle assembly. It also has a rotary table 
about the Zj axis. Thus it has three degrees of freedom for the blank and one degree of 
freedom for the cutter. Hence there are four entries viz., b^, 02 , V; in the table 


' LINTOL is a post-processor command wMch is an acronym for LINear TOLerance. The post- 
processor generates additional points between adjacent pairs of blocks of CL file in order to ensure 
that the path followed by the tool tip does not deviate beyond this tolerance band. 

^ The post-processing consists of two step. The first step is the mapping of CL data into the joint 
coordinate space and the second one is to generate output in the required format. The second step 
will still be required to be performed. 
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which are time-dependant and the rest are constant values. It may be noted that Gj which 
is related to the spindle rotation is entered as zero for all milling machines since it is a 
passive degree of freedom. 


3.2.2 Instantiation of a 5 Axis Horizontal Milling Machine with Two Rotary Tables 

A 5 axis CNC machine with the table capable of moving in the X 2 and Y 2 
directions and the spindle head contributing to the Z, movement is’ shown in Figure 3-2. 
This machine has a rotary table about the axis in addition to the conventional rotary 
table about the axis. Thus it has four degrees of freedom for the blank and one degree 
of freedom for the cutter. Hence c,, a^, b^, 62 and y, are the five time-dependant values 
and the rest are constants. 


3.2.3 Instantiation of a 5 Axis Milling Machine with a Rotary Tables and a Tilting 
Tool Head 

Another 5 axis CNC machine shown in Figure 3-3 has the X^ and movements 
for the table and the Z, movement for the spindle head. Apart from the conventional 
rotary table about the Zj axis, its spindle head also can be tilted about the Z, axis. Thus 
there are three degrees of freedom for the blank and two degrees of freedom for the 
cutter. Hence Cj, Gj, a^, and yj the five time-dependant entries in the table while 
the other seven are constants. 


3.2.4 Instantiation of a 5 Axis Milling Machine with Spindle Head Swivel 

Capability about Two Principal Directions 

The 5 axis CNC milling machine illustrated in Figure 3-4 is very similar to the 
machine shown in Figure 3-3 except that it has a tilting mechanism for the spindle head 
about the 7, axis instead of the rotary table about the Z 2 axis. Thus there are two degrees 
of freedom for the blank and three degrees of freedom for the cutter. Hence c,, G,, (}>,, o- 
and ^2 the five time-dependant entries in the input table while the other seven are 
constants. This machine is capable of producing two identical parts simultaneously since 
it has got two spindle heads. However, in this modeling, we are concerned with the 
spindle head on the right side only. 





Instantiation of the Generic Machine Tool for a 4 Axis Milling Machine with a Rotary Table. 

















Cutter Blank 
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Swivel Capability about Two Principal Directions. 
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3.2.5 Instantiation of a 5 Axis Milling Machine with Spindle Head Swivel 
Capability about Two Arbitrary Directions 


The four machine tools illustrated so far have their rotary motions about the 
principal directions. Since the ultimate aim of the two rotary movements in any machine 
tool is to orient the cutter in any desired direction, one can also think of machine tools 
with these two rotary motions about any two arbitrary axes. One such machine is shown 
in Figure 3-5. In this machine, the table is grouted to the floor and hence it has no degree 
of freedom. The cutter has all the five degrees of freedom. One of the rotary movements 
is about the A axis which is parallel to the Zq axis. The other rotary movement is about 
an axis called B axis which makes an angle 5 with the A axis. For this machine, a,, Z),, 
C], 6, and <)>, only will be the time-dependant variables. Yi will be a function of 6, and 
<{), and the rest six parameters will be constants. It may be noted that the purpose of B 
axis is to cause a throw or a tilt of the tool axis from its initial vertical position. Hence 
the effect of rotation about A axis is not felt when the tool axis is vertical since it will be 
equivalent to a change in the spindle rotation at this condition. If a is the angle of rotation 
about the A axis and j3 is the angle of rotation about the B axis then, a and P are related 
to the angles 6,, ({>, and \|r, as follows; 


01 = 0 

<Pj = cos”^[l - sin^6(l - cos P)] 

i|r - 5 sin g (1 - cos p) - cos g sin p 

^ [cos 6cosa(l-cosP)+smasmP 


Figure 3-6 describes the configuration of these two rotations. Figure 3-7 shows the 
cutter orientation for various values of p for a = 0. When P = 0, the rotation about the 
A axis has no influence in the orientation of the cutter. When the B axis is rotated, the 
cutter describes a tilted cone of fixed base diameter, say cone B, causing the cutter to 
deviate from its vertical position. When the A axis also rotates, the cone B itself rotates 
about the A axis. The enveloping cone may be called cone A whose apex angle is a 
function of p. Figure 3-8 shows the orientation of the cutter along with the corresponding 
cone A and cone B for various combinations of a and p. For these illustrations, 5 is 
assumed to be 7t/8. 


For a given cutter and the required surface Zj, the path generated by the conjugate 
geometry model will be directly in terms of the joint coordinates x, y, z, a and 6. In the 
case of the conventional APT method, the CL file will be generated which will contain 
the path information in terms of x, y, z, /, j and k where i, j and k are the directional 
cosines of the cutter axis. The post processor will have to map these into the joint space 
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using the following relations obtained from Equation (3-1): 

+ Icosflsind 

^ L^asmb 

^«w = 


2 tan"^ 


. b 

sm — 

2 




• 2 ^ • 2 ^ 
sm^ 0 - sin — 

2J 


71 . 

= — + 0 + tan *' 

new 2 


(cos 5 (1 - cos 


sinp„ 


(3-2a) 


(3-2b) 


where 


a = tan‘* ^ 
i 

b = cos'^ k 

and L is the distance between the cutter tip and the fulcrum point and the subscripts new 
and prev respectively refer to the new and the previous joint coordinate positions. The 
APT processor generates CL data such that between any two consecutive positions, when 
X, y, z, cos ' j, cos'^y and cos''ic are linearly varied, the resulting surface will represent the 
desired surface within the TOLER, INTOL and OUTOL^ limits. However, that the desired 
surfaces will be produced within these limits cannot be guaranteed when the joint 
coordinates x, y, z, a and P of the corresponding pair of CL data are linearly varied. In 
order to get the correct surface, the post-processor is required to calculate additional 
points between the two end positions and inserts them in the NC program. 


3 


INTOL, OUTOL and TOLER are APT commands for specifying the tolerance of the generated 
surface with respect to the desired surface. These stand for INside TOLerance, OUTside 
TOLerance and TOLERance. 



Kinematic 

Symbolic 
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Swivel Capability about Two Arbitrary Directions. 








Figure 3-6 SjiNndle Configuration of the 5 Axis Profiler Shown in Figure 3-5. 
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3^ Illustrative Example 

Let us consider an example of machining a helicoidal surface using a horizontal 
milling machine with an end milling cutter of the required form. The biparametnc 
definition of one such cutter in u and v is as follows: 

= tt cosv 

(3-3) 

3'j = « smv ' 

= au^ + bu + c 

where 

0 i ^ 5.33903^+01 
0 s V ^ 2 t: 

a = 1.68979^-02 
b = 1.02956^-06 and 
c = 4.91379e+01 

Figure 3-9a shows the kinematic scheme of this process. The cutter is stationary 
and the blank rotates about Xq axis at the rate of © radian/second while it translates along 
-Xq axis at the rate of p(£)f2ii meter/second where p is the pitch of the helicoidal surface. 
The twelve parameters of motion for this case will be: 


Cutter 

Blank 

a, =0 

pm 
a = -£ 

^ 1% 

II 

o 

b, = 0 

c, =0 

c, =0 

JX> 

II 

1 

0j = (£>t 
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o 

li 

o 

li 

< 

II 

o 

II 

o 


Table 3-1. Parameters Describing the Relative Motion in a Milling Process 

Having known the surface of the cutter L, and the twelve parameters of motion for the 
cutter and the blank, one can predict the machined surface 2^ using the algorithm 
described in the previous chapter. 

Step 1: The biparametric definition of the cutter surface is given below; 

^P=[^i 1 ] 


^/> = [« cosv u sinv au^ + bu + c l] "(3-4) 


Step 2: The twelve parameters of motion are described in Table 3-1. 


Step 3: 


The transformation matrix [ ] for transforming the position vector 

and the normal vector ” j>,ifrom the cutter coordinate system Sj to the 
global coordinate system Sq is 



10 0 0 

0 0-10 

0 10 0 

0 0 0 1 


(3-5) 


Matrix [ ] is the time-differential of the matrix jM 

the position vector V in the cutter coordinate system 
vector i>,i in the global coordinate system Sq. 


which converts 
Ij into velocity 
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0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 


(3-6) 


Step 4: 


The transformation matrices | ^ j and ^ j for transforming the 
position vector ^ from the blank coordinate system S; to the global 
coordinate system Sq and vice-versa are respectively: 



1 0 0 
0 cos o)f sin oir 

0 -sin (at cos cof 

0 0 

2 TC 


0 

0 

0 

1 


(3-7) 



10 0 0 
0 cos (Of -sin (Of 0 

0 sin (Of cos (Of 0 


0 0 1 

2 % 


(3-8) 


Matrix J is the time-differentialof the matrix | ^ 1 which converts 
the position vector P in the blank coordinate system Sj into velocity 
vector X global coordinate system Sq. 
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Step 5: 


[«M] = 


0 

0 

0 

po) 
' 2t. 


0 

of 

-cocos (Of 


0 

cocos (Of 
-CO sin (Of 


0 

0 

0 


(3-9) 


The position vector P is 

3-2 ^2 1] 




= 


P (0 f 

u cosv + ^ 

2 % 


{au^ + bu + c)cos(of-«smvsm(of 
-(o + bu + c) sin (0 f - M sin V cos co f 
1 


(3-10) 


The normal vector of the cutter surface is obtained as follows; 


d^p 



= [ COS V 

du ^ 

sin V (2 au + b) 

■0] 




= r -« sin V 

dv '■ 

« cos V 0 

0] 


Step 6: 
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Step 7; 


Step 8: 


d'p d'p 




p.i 


du dv 






-(lau + b) u cos V 
a u + i) M sin v 


IT 


U 

0 


(3-11) 


The unit normal vector ^ , of the cutter surface is obtained as follows: 


I I = ^ \ + Q.au + bf u 


p,i 




f,i 





Q.au + b) 

/l + (2 au + bf 
(^au + b) 

^ \ + Q,au + bf 


cos V 


sin V 


1 


f 1 + (lau + bf 
0 


(3-12) 


The normal vector '’w p_i of the cutter surface in the global coordinate 
system is; 
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Step 9: 


Step 10: 


“",.1 [>] 


”p,i 


(2 AM + b) 


/ 1 + Q,au + bf 

1 

\/ 1 + Q. a u + bf 
{2au + b) 


cos V 


yj 1 + {2 au + bf 


0 


sm V 


(3-13) 


The relative velocity vector f .,12 in the global coordinate system can 
be now calculated as: 


“F,, -■?[>}] 

=[0000] 


(3-14) 


oy 


P ,2 


owsinv a (ft + b u + c) 0 
2k 


(3-15) 


= Oy - Op 
p,n P.2 p,i 


P.12 


COMsinv 0 


2 Jl 


(3-16) 


The condition of contact p j * p j 2 - 0 yields the following 
relation between u and v: 
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tan V 


u + ^au + b) {a + bu + c) 


(3-17) 


Step 11: Surface can be now represented as a biparametric equation in u and t 

as given below which can be graphically represented as shown in Figure 
3-9b. 

pot 

a:, = « cosv + - — 

2 n 

y^= (au^ + bu + c) cos ot - u sin V sm ot 
^2 = -(fl «^ + & M + c) sin £0 f - « sin V cos 0) r 

where 


tan V 


u + Q,au + b) {a + bu + c) 


Incidentally, this happens to be the screw surface of a twin-screw extruder used in the 
polymer industry. This is a bi-lobed screw with zero minor the cross-sectional profile of 
which is shown in Section A-A of Figure 3-9c. 




Chapter 4 


SYMBOLIC DEFINITION OF A GENERIC CUTTER 


4.1 Introduction 

The geometry of a cutter can be modelled as an axisymmetric, biparametric 
surface. The axis of symmetry depends on whether the cutter is an end mill type cutter, 
a side mill type cutter, a face mill type cutter or any combination of some of these basic 
types. In the present work, a model of a generic cutter is developed in which the 
generatrix curve is a planar curve consisting of a set of contiguous straight line and 
circular fillet segments. The directrix motion is the axisymmetric sweep of this curve 
around an axis in the plane of the generatrix curve. In numerical control literature, a 
generic cutter is defined using seven parameters and the generatrix curve consists of two 
straight line segments forming the tip flank as well as the side flank and a circular fillet 
between these straight line segments (Chang, 1989). In this thesis, an extension of this 
concept is used. The generatrix curve is assumed to be a piecewise continuous curve 
consisting of three straight line segments and a circular fillet between the tip flank and 
the side flank lines (see Figure 4-1). This generatrix curve can be defined by means of 
eight parameters, viz., d, r, e,f, a, b, h and A, as shown in Figure 4-1. The cutter surface 
can be obtained by rotating the generatrix curve either along the Z, axis or along the X, 
axis as shown in Figure 4-2. When the sweep is a rotation about the Z, axis, one can 
generate the geometry of various types of end milling cutters used in 3 to 5 axis 
machining and also the cutters used in the non-traditional machining processes such as 
laser beam, water jet and wire of an EDM process. If the rotation is about the axis, 
then a disc type of cutting tool geometry is generated which will represent a side-and- 
face milling cutter or a grinding wheel. If the generatrix curve is swept in a helical 
manner, then the geometry of a gear hob or a tapping tool can be obtained. Thus by 
choosing appropriate values of the parameters describing the planar curve and the sweep 
motion, one can obtain a variety of conventional and non-traditional cutting tools. 





Zc, 



Figure 


4-1. Definition of a Generic Cutter 
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0 Sweep Matrix to Generate an End Mill Type Cutters (Rotation about Z, axis) : 



’Im] 


cos V sin V 0 
-sin V cos V 0 


0 0 


(ii) Sweep Matrix to Gnerate a Disc Type Cutters (Rotation about an Axis Parallel to X, 
axis and at a distance of q from X, Axis) : 



’Im' = 


cos V sin V 0 
-sin V cos V 0 


0 0 


(iii) Sweep Matrix to Generate a Gear Hob or a Tap Type Cutters (Helical Sweep along X^ 
axis with Lead / and Radius of Helix q) : 


Dll 


10 0 0 

0 cos V sin V 0 

= 0 -sin V cos V 0 

Iv . , 

q sm V -q cos v 1 


Figure 4-2. Examples of the Sweep Matrix 
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4.2 Description of the Generic Cutter 

The planar curve of the generic cutter surface consists of three straight line 
segments and a circular arc segment (Figure 4-1). For the case of an end mill type cutter, 
the first and second straight line segments will generate the tip and the side flanks of the 
cutter while the third straight line segment will generally represent the shank of the cutter. 
The circular arc segment will generate the comer radius of the end mill. A coordinate 
system (O^ - X„ ZJ is considered to be attached to this cutter profile. This planar 
curve of the cutter profile can be represented by means of a parameter u which is the 
profile length of any point on the profile with respect to the origin O^. The cutter surface 
Xj is obtained by sweeping this generatrix profile along a directrix curve which can be 
circular, helical, etc. This sweep motion is represented by means of a homogeneous 
transformation matrix whose elements will be functions of a parameter v. The four 

segments of the generatrix curve can be algebraically expressed as 


0 / 

PQ 

QS 

ST 


Zg = tan a 

(x^ - ef + (z, -f) = 


= z^tnn b - — tajiatanZ> + — 
^ ^ 2 2 

= AtanZ>--tanatanZ> + — 

‘'2 2 


(4-1) 


where 

- First straight line segment, 

PQ - The circular segment with center at A, 

QS - Second straight line segment, 

ST - Third straight line segment, 

R - Point of intersection of the first segment and the third segment 
d - Diameter of the cutter measured at R, 

r - Radius of the circular segment. This is positive for a convex arc and 
negative for a concave arc, 
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a - Angle of segment measured from axis. -Jc/2 < a < Till, 
b - Angle of segment QS measured from axis, -n/2 < b < 7t/2, 
e ~ X coordinate of point A in the coordinate system S^, 

f - Z coordinate of point A in the coordinate system 

h - Nonparallel height of the cutter, 


Aj - Parallel height of the cutter, 

is Jillet radius - Flag which is to be set to TRUE when the circular segment PQ is 
tangent to both the straight line segments OJ’ and QS and FALSE 
otherwise. 


u - Parameter which is the profile length of any generic point P on the 
profile measured from O^, 


V - Parameter which describes the sweep of the two dimensional profile. 



Homogeneous transformation matrix which describes the sweep of the 
two dimensional profile. 


s, - Length of the first straight line segment OJP, 

Sj - Length of the circular arc segment PQ, 

- Length of the second straight line segment QS, 


- Length of the third straight line segment ST, 

‘p _ Position vector of a generic point P on the two dimensional profile of 
the cutter in S^ coordinate system in terms of u and 

- Position vector of a generic point P on the surface of the cutter in 5, 
coordinate system in terms of u and v. 

The parametric equation of the generatrix curve P W can be defined in a generic 

manner with respect to a coordinate system S^ {O^ - X„ ZJ as- follows ; 
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"p(k)=[^c(«) 1 ] 


where, 

for the interval Q ^ u < s^, 
x^{u) = M cos a 
yc(«) = 0 
(«) = « sin a , 


for the interval ^ u < s^+s. 


2 > 


x^(u) = e + afc(r) cos^'^l 

yAu) =0 


5j sin 0 - f^ 


Zc («)=/+ abs(r) sin^ 


,-ii 


cos a - e 


sin a - 


Si cos a - e 


u - s. 


u - s. 


for the interval S1+S2 ^ u < 5 ^+ 53 + 53 , 
xA^) = [« - sin Z> + ;«/ 

JcC") =0 

Zc(**) = [« “ cos + z/ 

and for the interval ^ w ^ Si+s^+s^+s^, 

XcO*) = x^ 1 

3 ’c(«) ^ 0 

Zc(«) = [« " (-Si+Sj+^a)] + z! . 


Having known the planar profile ‘P (“) and the sweep, the 
definition of the cutter surface Xj can be obtamed from the equation 


(4-2) 


(4-3) 


(4-4) 


(4-5) 


(4-6) 


biparametric 
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‘p(w, V) =*>(«) . 


(4-7) 


where [ J is the matrix describing the sweep motion. The algorithm of this generic 
cutter described in Section 4.3 has been implemented in the symbolic manipulation system 
MACSYMA and the results obtained for a 27 cutters are shown in Figure 4-3 to Figure 4- 
29. 


4.3 Algorithm to Derive the Biparametric Equation of the Generic Cutter 

The symbolic algorithm to derive the biparametric equation of the surface Z, from 
the nine parameters of the cutter profile and the sweep matrix has been implemented in 
MACSYMA. It may be noted that all the nine parameters of the generic cutter may not be 
required in all cases. For instance, when the flag is Jillet radius is TRUE and b is zero, 
only d,r,a,b and /?] are required and when is _fillet_radius is TRUE and b is nonzero, then 
only d,r,a,b,h^ and h are needed. The algorithm can now be stated as follows: 


1. Calculate e and /if is Jillet radius = TRUE. 

g = ^ - - sin 

2 \ cos (a+b) ) 

f = ^ + g sin fl 
cos a 


(4-8) 


2. Calculate point P. 

Xc = {e + ftm a) cosP a 

+/tan a) cos^of - cosPo 

p p . 

Zc =x, tm a 


(4-9) 
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Calculate point Q. 


^ = \f 


e + — (tan a tan t - 1) 


tan b\ cos^b 




- 1/ + 


e + (tan tan - 1) 


tan Z>| cos?Z> 

T2 1 




- r" 


cos^i> 


tan b - - (tan atmb - 
2 


4. Calculate h if is^Uet radius = TRUE and b = 0. 

h = z ® 


Calculate point S. 


= Atant-— tanatani> + - 
^ 2 2 

z! = h 


Calculate s,, 5 ,, S 3 and s^. 

= i/W * W 


s, =r< 


tan" 


V-/ 




- tan'* 


4-f 


p 

Xc - e j 


s, - - x,^^ ^ {z! 

^4 = *1 


(4-10) 


(4-11) 


(4-12) 


(4-13) 
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7. Define the parametric equation 'P («) of the two dimensional profile of the 
cutter in terms of the parameter u. 

'>(»)= yjM) z,iu) 1] (4-14) 

where, 

for the interval Q ^ u < 
x^iu) = u cos a 
ycOi) = 0 
z^iu) = M an a , 


for the interval ^ u < 5^+52, 


x^(u) = e + abs(r) cos 
yc(«) =0 


i tan'^ 


'' sin a - 
cos a - e ^ 



Zc(u) =f + ak 5 (r) sin 


tan"^ 


Sj sin a - 
Sj cos 0 - 


+ 



for the interval 3^+82 ^ u < 5^+52+ Sy 
= [« - (Si+S2+s0] sin + jc/ 
VcOO = 0 

Zciu) = [u - (S1+52+S3)] cos + z/ 


and for the interval ^ ^ Sj+S2+-S3+'S4> 

x^ (tt) = Xc 

VcOi) = 0 

Zc(«) = [» " (51+^2+ ^ 3 )] + z/ . 

8. Define the sweep matrix [ 

(See examples in Figure ^ 2 ) 
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9. Define V(«, v) . 

V(tt, V) = 'p(«) .[^m] (4-15) 


4.4 Instances of the Generic Cutter 



Figure 4.3 Instantiation of a Gear Hob Using the Model of the Generic Cutter 

Robbing is gear manufacturing process for mass production. A typical gear 
bobbing process has been illustrated in Figure 2-11. In this section, it will be shown how 
the parametric equation of one such gear hob can be obtained using the model of the 
generic cutter. It is to be noted that the surface that will be considered in this modeling 
is a continuous surface as shown in Figure 4-3 and does not represent the flutes. The 
tooth profile of the hob and the coordinate frames and S, are also depicted in Figure 
4-3. The tooth profile can be obtained by substituting the values given in the table 4-1 in 
the definition of the generic cutter. 
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Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

3.000 

0.000 

0.000 

20.000 

0.000 

(1.500) 

(0.000) 

6.500 


Table 4-1. Input Data to Instantiate a Gear Hob.. 

Thus the parametric equation of the cross-sectional profile of the hob is: 

"Piu)=[xc{u) y^iu) z,(«) 1] (4-16) 


where, 

for the interval 0 s « < 1.500, 
x^iu) = H 

yM=o 

Zc («) = 0 


and for the interval 1.500 s ^ 8.417, 
x^iu) = 0.342 (M - 1.500) + 1.500‘ 
Tc(«) = ^ 

Z,(«) = 0.939 («- 1.500) . 


This profile is given a helical sweep to obtain the surface of the hob. The transformation 
matrix describing this sweep is 



1 

0 

0 

0 

cos V 

sin V 

0 

-sin V 

cos V 

Iv 

^sin V 

-q cos V 

In 


0 

0 

0 

1 


(4-17) 


where I is the lead of the hob and q is the distance of the origin from the axis as 
shown in Figure 4-2c. Hence, the biparametric surface obtained from Equation (4-15) is 
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V («, v) = [ ;tj (m, v) y, (m, v) (m, v) 1 ] (4-18) 

where 

x^(u, v) =x^(m) + ^ 

2lt 

j, (w, v) = (u) cos V - (w) sin v + ^ sin v 

Zj (u, v) = (?/) sin V + z^ (i/) cos V - ^ cos v 

This definition of the cutter surface X, can be used in the conjugate geometry model of 
the generic machine tool described in the previous chapter in order to determine the gear 
profile produced on the blank. 

In addition to the example of a gear hob, There are 26 examples shown in Figures 
4-4 through 4-29. These are several types of drills, end mills, T-slot cutters, side and face 
milling cutters, flange cutters, wire of the EDM process, laser beam, etc. Each of these 
figures consists of three parts. The first part of is the graphic representation of the two 
dimensional profile of the cutting tool. The second part gives the table of the parameters 
defining the cutter. This will be the input data for the symbolic algorithm. It may be noted 
that the values within the parentheses need not be provided as the algorithm will calculate 
them from the tangency conditions. The third part in these figures is the output from the 
symbolic algorithm which is the biparametric equation of the cutter surface Pm terms 
of u and v. In all the examples, the lengths are given in millimeter and angles are given 
in degrees. 
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(a) Profile of the Cutter : 




(b) Input Data to the Symbolic Algorithm ; 


Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

25.000 

4.000 

0.000 

0.000 

66.000 

(8.500) 

(4.000) 

(4.000) 


(c) Output from the Symbolic Model : 


Range of u 

Xc 

B 

Zc 

0<u<8.500 

u 

0 

0 

8.500<u<l 4.783 


0 


- 

- 

- 

- 

14.783<u<80.783 

12.500 

0 

u - 10.783 


Figure 4.4 Instantiation of a 25 dia. 4R Dome End Mill 



















Figure 4.5 Instantiation of a 25 dia. OR End Mill 






























Chapter 4. Symbolic Definition of a Generic Cutter 


83 


(a) Profile of the Cutter : 



' (b) Input Data to the Symbolic Algorithm : 


Is arc 
tangential 

d 

r 

a 

b 


e 

f 

h 

yes 

25.000 

12.500 

0.000 

0.000 

57.500 

(0.000) 

(12.500; 

o 

o 

<N 


(c) Output from the Symbolic Model : 


Range of u 

Xc 

Yc 

2c 

- 

- 

- 

- 

0.000<u<] 9.635 

12.500 cos [0.080 u-1.571] 

0 

12.500 sin [0.080 u-1.571] 

- 

- 

- 

- 

19.635<u^7.135 

12.500 

0 

n- 7.135 


Figure 4.6 Instantiation of a 25 dia. Ball Nose End Mill 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm : 


Is arc 
tangentia' 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

25.000 

4.000 

45.000 

0.000 

55.843 

(8.500) 

(14.157 

(5.561) 


(c) Output from the Symboiic Model : 


Range of u 

Xc 

y< 

Zc 

0<ii<16.021 

0.707 u 

0 

0.707 u 

16.021<u<19.162 


0 

■■lygiiiimivraiiiimiaiHi 

- 

- 

- 

- 

19.162<u<75.006 

12.500 

0 

u- 5.006 


Figure 4.7 Instantiation of a 25 dia. 90 deg. Tool 
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Figure 4.9 Instantiation of a 25dia. 4R Side Angle Cutter 
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(a) Profile of the Cutter : 


(b) Input Data to the Symbolic Algorithm : 


Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

25.000 

0.000 

0.000 

5.000 

50.000 

(12.500; 

(0.000) 

20.000 


(c) Output from the Symbolic Model : 


Range of u 

Xc 

fl 

Zc 

0u<l 2.500 

u 

0 

0 

- 

- 

- 

- 

12.500<u<32.576 

0.087 (u-32.576) + 14.250 

0 

0.996 (u-32.576) -i- 20.000 

32.576<u<82.576 

14.250 

0 

u - 12.576 


Figure 4.10 Instantiation of a 25dia. OR Side Angle Cutter 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm ; 


Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

ll 

yes 

25.000 

4.000 

10.000 

5.000 

50.000 

(8.782) 

(5.610) 

20.000 


(c) Output from the Symbolic Model : 


Range of u 


B 

z. 

0<u<9.624 

0.985 u 

0 

0.174 u 

9.624<u<14.860 


0 


14.860<u<29.654 

0.087 (u-29.654) + 14.057 

0 

0.996 (u-29.654) + 5.610 

29.654<u<79.654 

14.057 

0 

u - 9.654 


Figure 4.11 Instantiation of a 25 dia. 4R Side & Bottom Angle Cutter 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm : 


Is arc 
tangentia' 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

25.000 

0.000 

10.000 

5.000 

50.000 

(12.500; 

(2.204) 

20.000 


(c) Output from the Symbolic Model : 


Range of u 

Xc 

B 

2c 

0<u<12.693 

0.985 u 

0 

0 

- 

- 

- 

- 

12.693<u<30.557 

0.087 (u-30.557) + 14.057 

0 

0.996 (u-30.557) + 20.000' 

30.557<u^0.557 

14.057 

0 

u - 10.557 


Figure 4.12 Instantiation of a 25 dia. OR Side & Bottom Angle Cutter 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm ; 


Is arc 
tangentia 

d 

r 

a 

b 

hi 

e 

f 

h 

yes 

25.000 

0.000 

-lO.OOC 

5.000 

50.000 

(12.500 

1 (-2.204; 

20.000 


(c) Output from the Symbolic Model : 


Range of u 

Xc 

B 

Zc 

'0<u<12.693 

0.985 u 

0 

-0.174 u 

- 

- 

- 

- 

r2.693<u<34.982 

0.087 (u-34.982) + 14.443 

0 

0.996 (a-34.982) + 20.000 

34.982<u<84.982 

14.443 

0 

u - 14.982 


Figure 4.14 Instantiation of a 25 dia. OR Side & Bottom Angle Cutter 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm : 


Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

25.000 

4.000 

10.000 

-5.000 

50.000 

(8.196) 

(5.507) 

20.000 


(c) Output from the Symbolic Model : 


Range of u 

Xc 

B 


0<u<9.028 

0.985 u 

0 

0.17411 

9.028<u<14.962 

4 cos [0.250 (u-9.028)-1.396] 

+ 8.196 

0 

4 sin [0.250 (u-9.028)-L396] 

+ 5.507 

■14.962<u<29.160 

- 0.087 (u-29.160) + 10.943 

0 

0.996 (u-29.160) + 20.000 

29.160<u<79.160 

10.943 

0 

u - 9.160 


Figure 4.15 Instantiation of a 25 dia. 4R Side & Bottom Angle Cutter 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm : 


Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

25.000 

0.000 

10.000 

-5.000 

50.000 

(12.500; 

(2.204) 

20.000 


- (c) Output from the Symbolic Model : 


Range of u 

Xc 

6 

Zc 

0<u<12.693 

0.985 u 

0 

0.174 u 

- 

- 

- 

- 

12.693<u<30.557 

-0.087 (u-30.557) + 10.943 

0 

0.996 (u-30.557) + 20.000 

30.557<u<80.557 

10.943 

0 

u - 10.557 


Figure 4.16 Instantiation of a 25 dia. OR Side & Bottom Angle Cutter 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm : 


Is arc 
tangential 

d 

r 

a 

b 

hi 

e 

f 

b 

yes 

25.000 

4.000 

-10.000 

-5.000 

50.000 

(8.061) 

(2.640) 

20.000 


(c) Output from the Symbolic Model : 


Range of u 

Xc 

B 

Zc 

0<u^.480 

0.985 u 

0 

-0.174 u 

■'7.480<u<14.810 


0 


14.810<u<31.886 

-0.087 (u-3 1.886) -t- 10.557 

0 

0.996 (u-3 1.886) + 20.000 

31.886<u<81.886 

10.557 

! 0 

u - 11.886 


Figure 4.17 Instantiation of a 25 dia. 4R Side & Bottom Angle Cutter 






















Figure 4.18 Instantiation of a 25 dia. OR Side & Bottom Angle Cutter 
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Figure 4.20 Instantiation of a 25dia. 4R Non-Fillet Type Side & Bottom Angle Cutter 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm ; 

ts arc , ~ 
. . d I r I a I b 

tangentia. 


yes 25.000 4.000 0.000 -5.000 50.000 (8.135) (4.000) 20.000 


(c) Output from the Symbolic Model : 


Range of u 


y, 

0<u<8.135 

u 

0 

■8.135<u<14.767 

4 cos [0.250 (u-8.135)-1.571] 

+ 8.135 

0 

14.767<u<30.478 

-0.087 (u-30.478) + 10.750 

0 

30.478<u<80.478 

10.750 

0 



0.996 (\i-30.478) + 20.000 


u-1 0.478 


Figure 4.21 Instantiation of a 25 dia. 4R Pick-feed Cutter 






















Chapter 4. Symbolic Definition of a Generic Cutter 


100 



































Chapter 4. Symbolic Definition of a Generic Cutter 


101 

































Chapter 4. Symbolic Definition of a Generic Cutter 


102 
































Chapter 4. Symbolic Definition of a Generic Cutter 


103 











Chapter 4. Symbolic Definition of a Generic Cutter 


104 



Figure 4.27 Instantiation of a Beam in a Laser Cutting Process 
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(a) Profile of the Cutter : 



(b) Input Data to the Symbolic Algorithm : 

ts arc . , ^ 1 I 

. .. ,d r a b h, e f b 


Is arc 
tangentia! 

d 

r 

no 

25.000 

-4.000 


(c) Output from the Symbolic Model : 


Range of u 

Xc 

0^u<8.500 

u 

8.500<u<14.783 

4 cos [-0.250 (u-8.500) + 1.571] + 
12.500 



14.783<u<80.783 



Figure 4.28 Instantiation of a 25 dia. 4R Concave Cutter 
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Figure 4.29 Instantiation of a Conical Cutter Used in 5 axis Machining 
































Chapter 5 


ILLUSTRATIVE EXAMPLES 


5.1 Introduction 

The conjugate geometry model described in Chapter 2 along with the model of a 
generic cutter explained in Chapter 4 can be effectively used for modeling a variety of 
manufacturing processes. A manufacturing engineer generally has a set of machine tools 
and a collection of cutting tools and he will be required to produce the given surface 
within the tolerance specifications at the minimum cost by appropriately choosing a 
machine tool and a cutting tool. These models will aid him in evaluating the various 
alternatives available to him without performing the actual cutting. In this chapter, the use 
of these models has been demonstrated by means of several examples in design and 
manufacture. 


5.2 Modeling of Gear Manufacturing Processes 

Gear is one of the most common machine elements. They are manufactured using 
forming processes such as casting and blanking. They are also produced by a variety of 
machining processes such as milling, shaping, planing, bobbing, wire EDM etc. The 
selection of the manufacturing process depends on the size and the shape of the gear, the 
quantity required and the accuracy dictated by the application. - 

The teeth of a gear can be cut either by a direct process or by a generating 
process. In a direct process, cutting is performed with form-cutters which must have the 
form of the space between the teeth. Figure 5-la shows the milling of a gear using a disk- 
type gear milling cutter and the milling of gear with a Diametral Pitch cutter is shown 
in Figure 5- lb. These two processes are direct processes. In the direct process, the size 
and the shape of the space between the teeth of a gear is different for different number 
of teeth even for the same pitch. Since it is uneconomical to have a separate cutter for 
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each combination of module and number of teeth, generally a set of cutters is used with 
each cutter catering to a range of module regardless of the number of teeth. Hence, the 
profile of the gear produced by such a cutter will be exact only for a particular value of 
the module and number of teeth and for all other values it will be approximate. Moreover, 
the direct method is very slow due to indexed motion. 


The generating processes use screw-type cutters with profiles identical to the basic 
rack corresponding to the pitch and module of the gear. This process is more accurate as 
it generates the correct involute profile regardless of the number of teeth. Furthermore, 
the cutting in most of the generating processes is continuous and hence more productive. 
Gear planing vdth a rack cutter, gear shaping and gear bobbing are some of the gear 
generating processes used for the mass production of gears. In this section, two gear 
generating processes, viz., gear planing with a rack cutter and gear bobbing will be 
illustrated using the proposed conjugate geometry model. 


5.2,1 Modeling of a Gear Planing Process with- a Rack Cutter 

Planing by means of rack type cutters is one of the oldest methods of gear 
production that is still extensively used for large gears. This is a true generating process 
since it utilizes the principle that an involute curve can be formed by a straight line 
generator when a circle is made to roll without slip relative to the generator. It is, 
perhaps, the simplest practical method of demonstrating the application of the basic rack 
to gear tooth generation. 

In this process of planing a spur gear, the blank is made to rotate about its axis 
at a suitable speed and the cutter is moved longitudinally at the same linear speed as the 
blank. The teeth of the rack are straight and are relieved to form die cutting edges. The 
cutting action is obtained by reciprocating the cutter parallel to the axis of the gear blank. 
During the working stroke the rotation of the blank and the longitudinal translation of the 
cutter proceed continuously while on the return stroke the cutter is moved out of contact 
with the work to prevent rubbing. After a tangential advance of the cutter for a distance 
slightly exceeding one pitch, rotation of the blank is stopped, the cutter moves back the 
amount it has advanced and the process is repeated until the blank has completed one 
revolution. The gear is formed by sinking the cutter deeper into the blank until the 
required depth and tooth thickness are obtained. 

The first step in modeling this gear planing process is to obtain the biparametric 
equation of the trapezoidal profile of the rack tooth shown in Figure 5-2. This equation 
can be obtained using the generic cutter definition for the values given in Table 5-1. 
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Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

3.000 

0.000 

0.000 

20.000 

0.000 

(1.500) 

(0.000) 

6.500 


Table 5-1. Input Data to Instantiate a Rack Cutter 


Since the coordinate frames i”, and in the definition of the cutter are apart by a 
distance hH along Y direction, the matrix for this case will be 



(5-1) 


Therefore, the biparametric definition of the rack cutter will be 



(5-2) 


where 


for the interval -1!^- — ^ ^ u <-— 
2 cos a 2 



z 


1 


= V 
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for the interval « <2!L 

2 2 



Zj = V 


and for the interval — < w < ^ + 


X, = 


= 


w 


+ « - 


M’ 


U - 


W' 


cos a - 


2 

sin a 
h 


cos a' 


Z, = V 


The axis of the gear blank is the Zj axis which coincides with the Zq axis of the 
global frame of reference. The blank rotates about its axis with a constant angular velocity 
of 0). The pitch line of the rack cutter is the axis and it is at a distance of r from the 
Xq axis. The rack translates with a velocity of -/to along theyT, axis. After cutting a tooth 
in this forward stroke, it lifts to disengage from the blank and retracts exactly the same 
distance. The twelve parameters describing the relative motion between the gear blank and 
the rack cutter in the forward stroke are given in Table 5-2. 


Cutter 

Blank 

a, = -r €0 / 

a, = 0 

b, =r 

o 

li 

c, =0 

o 

il 

0, =0 

0=0 

2 
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o 

li 

<t>2 =0 

o 

II 

> 

II 

e 


Table 5-2. Parameters Describing the Relative Motion in a Gear Planing Process 


The tables 5-1 and 5-2 fully describe the surface of the cutter X, and the 
parameters of motion. These two inputs are used by the symbolic algorithm to determine 
the surface Xj which is the surface of the gear produced in this process. The surface of 
a gear tooth is output as algebraic equations which are given below. There will be three 
equations corresponding to the three segments of the rack tooth, viz., the left flank, the 
tip flank and the right flank. 



(5-3) 


where 

due to the left flank of the tooth 


— 

+ 

U + 

sin a - r (at 

cos cot + 


( 

M + 


2 1 

w 

- + 

2j 

U + 

sin a -rat 

sin cot + 


2j 

< \ 
w 

W + 


2 

2 , 



1 

2j 




\ 

root 


sin a - — cos a 
2 


cos Qc + 
cos a + fr 



sin cot 
cos cot 


due to the tip flank of die tooth 



2 . 

2 . 


sin cot 
cos cot 


V 




Figure 5-3. Surface of the Gear Obtained in a Gear Flanini 
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due to the right flank of the tooth 


^2 = 


z, = V 


^ 

2 [ 2 } 

( w'" 

— + h/ - — 

2 [ 2 } 

H 


Sin a - r(ot 


sin a -rat 


cos (at + w 


sm 0)/ + M 


and M = — + 
2 


+ r © t sin a + — cos a 


cos a + 
cos a + |r 



sin ©/ 
cos ©r 


The gear surface described by these equations is plotted for w = 3, h = 6.5, r = 25, a = 
20°. It is shown in Figure 5-3. 


5.2.2 Modeling of a Gear hobbing Process 

Robbing is regarded as the most accurate as well as the most productive process 
of gear machining. It is a generating process in which the hob teefli trace out the involute 
curves. 


A typical gear hob and the schematic diagram of the gear hobbing process are 
shown in Figure 5-4a and 5-4b. The hob has basically the same shape as a screw with one 
or more threads and each thread is cut by a number of gashes to form the cutting edges. 
Both the gear blank and the hob rotate about their axes, the blank slowly and the hob 
faster. Since the hob has the shape of a screw, its threads appear to move in the direction 
of its axis so that they simulate the movement of a rack cutter. In addition, since the hob 
is rotating, the cutting faces move continuously in the tangential direction and this 
movement provides the cutting velocity. In order to cut each tooth across the entire face 
width of the gear, the hob is fed slowly in the direction parallel to the axis of the blank. 
The distance between the axes of the blank and the hob is set to get the required tooth 
depth on the blank. 

In order that the cutter faces should meet the gear blank at the same angle as a 
rack, the hob is tipped over by a small angle as shown in Figure 5-4b. This angle is 
known as the swivel angle. The swivel angle is adjusted so that it is equal to the lead 
angle a of the hob. The tangential velocities of the hob and the blank at the pitch 
diameter are same for spur gears. The helical gear is obtained by superimposing a 
compensatory motion to this basic motion by means of a differential gear arrangement. 
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Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

3.000 

0.000 

0.000 

20.000 

0.000 

(1.500) 

(0.000) 

6.500 


Table 5-3. Input Data to Instantiate a Gear Hob 

The modeling of the gear hob using the generic cutter is illustrated in Chapter 4. 
The tooth profile of the hob can be obtained by substituting the values given in Table 5-3. 
The tooth profile of the hob tooth is given in the plane. The surface of the hob is 

obtained by performing a sweep operation which consists of the following three 
transformations : 


a) a shift along Zj axis by a distance of -d-^H where rf, is the diameter of the hob, 

b) a twist about Z, axis by an angle of -o where a is the lead angle of the hob 
and c) a helical sweep about X, axis described by pitch p and the parameter v. 

Therefore, the sweep matrix [ ] for this case will be 



1 0 
0 1 
0 0 

0 0 


0 

0 

1 

2 


COSO 

sin o 



-sin a 

COSO 

0 

0 


0 

0 

1 

0 


1 

0 

0 

/v 

27t 


0 

cos V 
-sin V 

0 


0 0 
sin V 0 
cos V 0 

0 1 



cos o -sin o cos v 
sin o cos o cos v 
0 -sin V 

/ V ^1 • 

sin V 

27t 2 


-sino sinv 0 

COSO sinv 0 

cos V 0 

1 

COS V 1 

2 


(5-4) 


Hence, the biparametric surface of the hob obtained from the definition of generic cutter 
is 
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(5-5) 


where 


„ , . 1 w h . . w 

for the interval - — ^ u < ~ — 

2 cos a 2 


= 




/ \ 

2 I, 2 ] 


Sin a 


w I H' 1 . 

■ — + M + sin a 

2] 


2 

w 


( 

M + 

I ^ 


Sin a 


COSO + 


Iv 


2 n 

sin a COS v + 
sino sinv - 


j 

(J + u + — 
2 


\d, + |w + 




cos a 


cos a 


sin V 


cosv 


w w 
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The axis of the gear blank is the axis which coincides with the Zo axis of the 
global frame of reference. The blank rotates about its axis with a Yrn 

of o. Initially, the plane X, - 7^ coincides with the Xo - Yo plane and as die bl^k rotat s 
it rises ivith a linear velocity off. The axis of the hob is the J, a^as md it is a istanc 
of (d, + d^ - 2/i)/2 where dj is the outer diameter of the hob, is the outer 
the blank and h is the tooth depth of the gear. The hob is tilted about 7, y 
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constant angle of o. The only movement that the hob has is the rotation about its own 
axis with an angular velocity of «co where n is the number of teeth that will be cut using 
a single-start hob. The twelve parameters describing the relative motion between the gear 
blank and the hob for machining a spur gear are given in Table 5-2. 


Cutter 

Blank 

a. =0 

! 

«2 = 0 

/ \ 
d, + d, 

b= ' ^ - h 

'[ 2 ) 

=0 

c^ = ^ 

c, =ft 

0, = neat 

02 = 0 

(j), = a 

O 

II 

II 

o 

Xj/j = to t 


Table 5-4. Parameters Describing Relative Motion of Gear Hobbing Process 

The symbolic algorithm uses these data describing tiie hob surface and the relative 
motion to determine the surface of the gear generated. The equations describing this gear 
surface are given below ; 

= [^2 Tj ^2 ( 5 -^) 


where 


Due to the left flank of the tooth, 
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where the condition of contact reduces to the following two constrmnts ; 

|1 + sin (n w r + v)] = 0 and 

n (d2 - h) 


V = 


/ cosa 


due to the tip flank of the tooth 


u sin^ a [l + sin (n (o r + v)] 
d, 


X, = - 
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where the condition of contact reduces to the following two constraint 

ri + sin (n w r + v)] = 0 and 
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■ sin © / 
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Due to the right flank of the tooth. 
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where the condition of contact reduces to the following two constraints ; 


[1 + sin (n w t + v)] = 0 

V = - 


and 

It (dj - A) 0) t 
/ cosa 


Here again, we will get exactly the same equations for the surface of the gear as in the 
previous example of gear planing since the profile of the hob tooth and the rack tooth are 
same. The gear surface described by these equations is plotted for iV = 3, /i = 6.5, r = 25, 
a = 20°. It is shown in Figure 5-3. 


5.3 Simulation of Alternative Machining Strategies 

A manufacturing engineer generally has a set of machine tools and a collection of 
cutting tools and he will be required to produce the given surface within the tolerance 
specifications at the minimum cost by appropriately choosing a machine tool and a cutter. 
Often, there will be more than one way to produce a component. Sometimes, the engineer 
has to explore alternate ways to manufacture a component because of the unavailability 
of the machine tool and/or the cutting tool (for which the component is earmarked) either 
due to a break-down or due to a busy schedule. These alternatives will have to be 
evaluated under the constraints of accuracy, time and cost in order to select the most 
optimal one. The selection of the most suitable alternative will be greatly simplified and 
speeded up if one can simulate the process using computer-based models such as the 
conjugate geometry model proposed in this thesis and perform evaluation of these 
alternative machining strategies with the aid of computers. In this section, the simulation 
of the alternative machining strategies for machining a helicoidal surface and a hyperbolic 
paraboloidal surface is discussed. 


5J.1 Alternative Machining Strategies to Produce a Helicoidal Surface 

A helicoidal surface can be realized if the relative motion between die blank and 
the cutter consists of a rotation about the axis of the blank and a translation along the 
same axis. The manufacture of a helicoidal surface whose slot has a hemi-spherical cross- 
section across its helix is considered in this example. Three c^es of realizing this 
helicoidal surface using the machine tools shown in Figures 3-1, 3-2 and 3-3 are 
discussed below. In each of these cases, the description of the expressions for the twelve 
degrees of freedom and the helicoidal surface produced are given. 
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Case (i) : In this case, a disc type of cutting tool similar to the one shown in Figure 4-24 
is used on the 4 axis vertical milling machine shown in Figure 3-1. This disc type of 
cutter can be instantiated using the generic cutter model using the values given in Table 
5-5. 


Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

00 

5.000 

0.000 

-90.000 

50.000 

20.000 

5.000 

10.000 


Table 5-5. Input Data to Instantiate a Disc Type Cutter 
The resulting expressions for the cutter surface Zj are given below : 


'p = [x, p, z, 1] (5-7) 


where 


for the interval — - r ^ s - + r(it - 1) and 0 ^ v ^ 2 it 

2 2 


x^{u, v) = 

3’i(«5 v) = 


d 

. ^ 2^^ d 

r sm + - - r 

r 2 

' d ^ 

. “ 2 '’^ d 

f Sin + _ - /• 

T 2 


cos V 


sin V 


v) - r 


1 ~ cos 


d ^ 
2 


and d is the diameter of the cutter and r is its comer radius. 


The machining scheme is shown in Figure 5-5a. The blank is set on the rotaiy 
table such that their axes are collinear. The tool is held on the vertical spindle of the 
machine such that the axis of the disc is parallel to the table axis and offset from it along 
X direction so as to get the required depth of cut. The helical formative motion is 
obtained by rotating the rotary table and moving it along Z direction. Both these 
movements are coordinated to get the required lead of the helix. 
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Cutter 

Blank 

d+d 

a= - r 

’ 2 

a, =0 

ft 

o 




c, = 

c,-ll 
^ 2% 

e, =0 

e, =0 

II 

o 

I 

4)2 =0 

i! 

O 

V 2 = / 


Table 5-6. Parameters Describing Relative Motion in a Helical Milling Process - Case (i) 


From the above, the geometry of the machined surface Xj can be expressed as a 
biparametric equation in u and t as follows ; 

>] 


(5-8) 
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where 


for the interval f"'’ ^ ^ | and d ^ t ^ t. 
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sm V cos t 
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Z^iu, 0 = H 


1 - cos 


d ^ 
u--+r 
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11 

2 % 


sm V = 


-Kid^+d-lr) tan 


d 

u--+r 

2 


The plot of this surface is given in Figure 5-5b. 


Case (ii) : In this case, a ball-nose end milling cutter similar to the one shown in Figure 
4-6 is used on the 5 axis horizontal milling machine shown in Figure 3-2. This cutter can 
be instantiated using the generic cutter model using the values given in Table 5-7. 



Figure 5-6. MachinWg of a Helicoidal Surface - Case (il). 
(a) Schematic Diagram of the Process and 
fb) The Helicoidal Surface Obtained in This Process. 
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Is arc 
tangential 

d 

r 

a 

b 

h, 

e 

f 

h 

yes 

10.000 

5.000 

0.000 

0.000 

50.000 

(0.000) 

(5.000) 

10.000 


Table 5-7. Input Data to Instantiate a Ball-Nose End Milling Cutter 


The resulting expressions for the cutter surface I, are given below : 



T, 


'] 


where 


for the interval 0 ^ « 
v) = 
3’i(«. v) = 


^ and 0 
2 

r sin — cos V 
r 

. u . 
r sm — sm V 

r 


^ V ^ 2 


Zi(«, v) 



(5-9) 


The machining scheme is shown in Figure 5 -6a. The blank is set on the horizontal 
rotary table such that their axes are collinear. The cutter is held on the horizontal spindle 
of the machine such that the axis of the cutter is perpendicular to die table axis and the 
tip of the cutter is offset from the table axis along Y direction so as to get the required 
depth of cut. The helical formative motion is obtained by rotating the rotary table and 
moving it along X direction. Both these movements are coordinated to get the required 
lead of the helix. 


Cutter 

Blank 

a, =0 

pi 

fl, = 

^ 2n 
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2 

^ = 0 

o 

Ii 

S =0 

e, = -1 
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o 

11 


<1)2 = 0 

V|/, = 0 

Vj - 0 

1 


Table 5-8. Parameters Describing Relative Motion of Helical Milling Process Case (ii) 

From the above, the geometry of the machined surface y u.. u* • j 
, . . . c u can be obtained as a 

biparametnc equation in u and t as roiiows ; 


> = [i, y, 1] 


(5-10) 
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The plot of this surface is given in Figure 5-6b. 
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Case (iii) : In this case also, the same ball-nose end milling cutter used in Case (ii) is 
used. However, the machine tool used here has a rotary table which has its axis vertical 
as shown in Figure 3-3. The cutter is described by Table 5-9 and the following equation 
which is exactly same as before. 


Is arc 
tangential 

d 

r 

a 

b 

h. 

e 

f 

h 

yes 

10.000 

5.000 

0.000 

0.000 

50.000 

(0.000) 

(5.000) 

10.000 


Table 5-9. Input Data to Instantiate a Ball-Nlose End Milling Cutter 


The resulting expressions for the cutter surface E, are given below ; 


= [-. 

where 

for the interval 0 s h ^ 
x^iu, v) = r 
3'i(«, V) = r 
v) = r 


y, 2, 1] (5-11) 


and 0 s v ^ 2 rt 

2 

. u 

sm — cos V 
r 


. u . 

sm — sm V 

T 



The machining scheme is shown in Figure 5-7a. The blank is set on the vertical 
rotary table such that their axes are collinear. The cutter is held on the horizontal spindle 
of the machine such that the axis of the cutter is perpendicular to the table axis and the 
tip of the cutter is offset from the table axis along Y direction so as to get the required 
depth of cut. The helical formative motion is obtained by rotating the rotary table and 
moving it along Z direction. Both these movements are coordinated to get the required 
lead of the helix. 
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Figure 5-7. Machinirig of a Helicoldal Surface’- Case OH)- 
(a) Schematic Diagram of the Process and 
fbl The Helicoldal Surface Obtained in This Process. 
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Table 5-10. Parameters Describing Relative Motion of Heycal Milling Process - Case Oii) 

From the above, the geometry of the machined surface S, can be expressed as a 
biparametric equation in u and t as follows : 

-P = [^2 ^2 ^2 1 ] ( 5 - 12 ) 


where 


for the interval 0 s « s and 0 s t s t, 
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and 


It 

tan V = ^ 

P 

The plot of this surface is given in Figure 5-7b. 

In case (ii) and case (iii), the helicoidal surfaces obtained are identical and they 
are exactly same as that of the desired surface which is a function of only the nose radius 
of the cutter. However, the surface obtained in case (i) is influenced by both the diameter 
d and comer radius r and hence it is an approximation to the desired surface. It can be 
seen that as the diameter d of the disc cutter approached twice the value of its fillet radius 
r, (i.e., d tends to 2r), the error between the desired and produced surface decreases. 


53.2 Alternative Machining Strategies to Produce a Hyperbolic Paraboloidal 

Surface 

When a block of material is moving along a straight line with constant velocity, 
if an end mill type cutter is tilted with constant angular speed in a plane perpendicular 
to the movement of the blank, die resulting surface is called a hyperbolic paraboloid. Two 
alternative methods of realizing this surface is discussed in this section, using the 5 axis 
machines shown in Figures 3-4 and 3-5. These are described in Case (i) and Case (ii). 

Case (i) : In this case, the cutter used is an end milling cutter. Only the side of the end 
mill is in contact with the material. The cylindrical surface of this cutter can be expressed 
in parametric form as 


'P = [^1 >’i 1] 


(5-13) 


where 

for the interval Q ^ u ^ hj + r and 0 ^ v ^ 2 % 
jCj(u, v) = r cos V 
yj(tt, v) = r sin V 

^ Z^Oi, v) = u - hj- r 


Figure 3-4 shows the machine tool on which this surface is machined. It is a 5 axis 
machine with the tool head which can be tilted about the two principal directions X and 




Figure 5-8. Modeling the Machining of a Hyperbhiic Paraboloidal Surface - Case (i). 
(a) Schematic Diagram of the Process and 
(b) The Hyperbolic Paraboloidal Surface Obtained In This Process. : 
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Y. The scheme of machining is shown in Figure 5-8a for this case. The blank is mounted 
on the table and the table moves along the X axis with unit velocity. At the same time, 
the cutter that is mounted on the spindle head swivels at an angular velocity of ti/k about 
the X axis. The twelve degrees of freedom for this process is described in Table 5-11. 


Cutter 

Blank 

a, =0 

a, = -t 

b, =0 

b, =0 

li 

O 

11 ■ 


e, =0 

o 

il 

<1>2 =0 

¥, = 0 

o 

It 


Table 5-11. Parameters Describing Relative Motion of a Milling Process - Case © 

From the above, the geometry of the machined surface can be expressed as a 
biparametric equation in u and t as follows ; 
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where 


and 


for the interval Q ^ u ^ hj + r and ^ t ^ 
x^{u, t) = r cos V + t 

y^iu, r) = r sin V cos - (w - - r) sin 

K K 

z,(tt, t) = r sin V sin — + (« - - r) cos — + hf 

K ■' K ■' 


tan V = 


{U - - f) TZ 

This hyperbolic paraboloidal surface is shown in Figure 5-8b. 


Case (ii) : In this case, the cutter used is same as that of Case (i). Here also, only the side 
of the end mill is in contact with the material. Hence the equation of the cutter surface 
remains the same. However, the machine tool used here is a different one. Figure 3-5 
shows the machine tool on wdiich the same hyperbolic paraboloidal surface is machined. 
It is a 5 axis machine with die tool head \^dlich can be tilted about the two axes called 
A and B axes. The kinematic arrangement of this machine has been explained in detail 
in Chapter 3 by means of Figures 3-5 through 3-8. The scheme of machining is shown 
in Figure 5-9a for this case. The blank is mounted on the table. Since the table is fixed, 
the linear movement of the blank is replaced by the movement of the tool head in the 
opposite direction. At the same time, the cutter that is mounted on the spindle head 
swivels at an angular velocity of ti/k about the X axis. The twelve degrees of freedom for 
this process is described in Table 5-12. 


Cutter 

Blank 

a, =/ 

a, =0 

A, =0 

o 

II 

It 

c, =0 
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a _ Ttt 

o 

II 

<x> 

o 

II 

(j), =0 

o 

11 

> 

o 

II 


Table 5-12. Parameters Describing Relative Motion of a Milling Process - Case (ii) 


Since the rotary axes of the machine are not parallel to the principal directions, in order 
to get a tilt of the cutter about the X direction, the following relation between the a and 
P values (which are the joint coordinates pertaining to the rotary axes A and B of the 
machine tool) should be maintained. 


P = cos 


a = tan 


-1 


1 - cos 


Tt t 


1 - 


i 


sin^fi j 

cos 5 (1 - cos P) ' 
sin p ^ 


(5-15) 


From the above, the geometry of the machined surface can be expressed as a 
biparametric equation in « and t as follows : 

= [x, 1] (5-16) 


where 


for the interval 0 s « s + r and ^ t <. 
x^iu, t) = r cos V + t 
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Figure 5-9. 'Mbdeling the Machining of a Hyperbolic Paraboloidal Surfade r Case (ii)- 
(a) Schematic Diagram of the Process and 
(b): The Hyperbolic Paraboloidal Surface Obtained In This Process. 
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and 


tan V 


jC 

(u - - r) Ti 


This hyperbolic paraboloidal surface is shown in Figure 5-9b. 

Since the machine tool used in Case (i) has one of the rotary axes about the X 
direction, the hyperbolic paraboloidal surface can be machined exactly with a single NC 
block using this machine. However, several cutter path steps are required for the profiler 
used in Case (ii) in order to get the same surface to the required accuracy. In Case (ii), 
for each intermediate value of t, a,, a and p should be calculated wliereas in Case (i), the 
values of ^nd Gj are to be calculated only for the end values of t. Even though the 
programming effort is more for the machine tool shown in Figure 3-5, its tool head is 
very compact and the working space of the cutter is axisymmetrical about the cutter axis 
whereas the tool head is more bulkier and the range of tool reach is not axisymmetrical 
for the machine shown in Figure 3-4. 


5.4 Design and Manufacture of a Cylindrical Cam and Follower Mechanism 

A cam mechanism is a three link mechanism with the cam and the follower in 
higher-pair contact. Since the contacting portions of the cam and follower surfaces are 
conjugate to each other, one can use the conjugate geometry formulation proposed in this 
thesis to analyze these mechanisms. As the proposed algorithm is symbolic, one can 
simulate a variety of two- and three-dimensional cam and follower mechanisms simply 
by changing the input values. In fact, the idea of attaching moving coordinate frames to 
the conjugate surfaces which is the basis of the proposed conjugate geometry formulation 
was conceived initially for the analysis of these mechanisms. In this example, a 
cylindrical cam with a cylindrical translating follower is considered. 

For the purpose of modeling the cam and follower mechanism using the generic 
model of the machine tool, the follower and the cam are viewed respectively as an end 
milling cutter and a cylindrical blank respectively. Therefore, the follower's biparametric 
representation can be given as 

'p = [^, z, 1] (5-17) 
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where 


for file interval 0 ^ u ^ 2% and 0 ^ v ^ Vj 
Xi(u, v) = r cos u 
yi(«, v) = r sin « 

Zi(tt, v) = V 


The axis of the cam and the Zq axis coincide and the cam is rotating about its axis 
with an angular velocity of xrnity. The axis of the follower (Z, axis) is parallel to and 
offset from the Tq 3xis by a distance of -a. It has the freedom to translate parallel to the 
Zq axis. The translation required in this case is assumed to be sinusoidal with an 
amplitude of s/2. The motions of these two elements of the mechanism can be described 
using the values given in the Table 5-13. 


Follower 

Cam 

aj = -a 

S = 0 

o 

ii 


s • . 

c, = 0 

c=~ sin t 
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II 

CD 
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Ii 

O 

2 


o 

II 

% -0 
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II 

111, =/ 


Table 5-13. Parameters Describing Relative Motion of a Cylindrical Cam & Follower 
Mechanism 

From the above data, the surface of die cam profile as determined by die conjugate 
geometry model is given below ; 
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(5-18) 


where 

for the interval 0 s v ^ Vj and 0 ^ t ^ 2% 
t) = r cosM cosr - v sinf - a cost 
y^{u, t) = -r cosu smt - v cost + a sint 

t) = rsmB + — sint 
2 

and 


tan u 


2v 

5 cost 


This surface is shown in Figure 5-11. 

It is interesting to note that the manufacture of this cam profile can be modeled 
using the machine tool shown in Figure 3-2 or Figure 3-3 using exactly the same set of 
inputs where the follower will be rqilaced by an end milling cutter. 




Chapter 6 


DESIGN OF EXTRUDER SCREW SURFACES 


6.1 Introduction 

So far, we have discussed how conjugate geometry problems can be handled using 
symbolic, computational and graphical approaches. There are several types of free-form 
surfaces whose manufacturability needs to be evaluated before their design can be 
considered acceptable. One such case is that of design and manufacture of helicoidal 
surfaces (Saybold, 1992; Booy, 1978; Janssen, 1978; Battle et al., 1984). Flutes of twist 
drills, extruder screw surfaces and cylindrical cam surfaces are some examples of 
helicoidal surfaces. In the present chapter, it will be shown how helicoidal surfaces are 
designed keeping in mind the manufacturability constraint using the conjugate geometry 
model proposed in this thesis. In the next chapter, the specific problem of designing an 
end milling cutter and a side milling cutter for machining a helicoidal surface has been 
solved. The discussions in these chapters are focussed on the extruder screw surfaces 
which are used extensively in the polymer industry. 


6.1.1 Definition of a Helicoidal Surface 

Let 5 (“) be a 2-D profile in X - Y plane which is the generatrix. When this 
generatrix is given a helical sweep, the surface generated by the edges of the generatrix 
profile will be a helicoidal surface. When the generatrix has rotated by an angle v about 
Z-axis from its initial position, at that instant, it would have moved by a distance pv/ln 
along Z-axis where p is the pitch of the helix. Therefore, the helicoidal surface generated 
by the generatrix 8 (m) can be expressed in homogeneous coordinate system as; 
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P («, v) = giu) 


cosv sm V 
-sinv cos V 


0 

0 


0 

0 


0 

0 

1 

£V 

2n 


0 

0 

0 


(6-1) 


Figure 6-1 shows a helicoidal surface obtained by sweeping a straight line generatrix 

g = [ u 0 0 1 ] 25 ^ u i 15 ( 6 - 2 ) 


and an elliptical generatrix 

5=[acos« !>sm« 


1 ] 


(6-3) 


for a = 1 00, b = 50 and p = 200. 


6.1,2 Description of Extruders 

Twin-screw extruders are used in the polymer industry for extrusion and pumping 
of mass. In these equipments, one screw wipes its mate and vice-versa. It is this wiping 
action that makes the twin-screw' equipment attractive for handling many polymers since 
it eliminates dead spots where polymer can collect, stagnate and return degraded to the 
mainstream. 

The most common tv.in-screw extruder comprises a housing or a barrel, and two 
parallel shafts equipped with identical or opposite-hand screw or paddle elements (Figure 
6-2). The identical screws rotate in the same direction at the same speed and the opposite- 
hand screws rotate at the same speed but in opposite directions. However, twin-screw 
equipments with unequal number of lobes and hence unequal speed ratio and/or unequal 
diameter ratio are not uncommon. The cross-section of a typical screw is shown in Figure 
6-3. Each lobe of the screw consists of a tip segment, a root segment and two flank 
segments connecting these two segments. 

Based on the geometry of the pair of screws, these twin-screw equipments may 
be classified into three categories. They are; 
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Section 
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Figure 6-3. Cross-Section of a Typical Extruder Screw 


i) Twin-screw equipments with identical screws, 

ii) Twin-screw equipments with screws of unequal velocities but same diameters 
and 

iii) Twin-screw' equipments with screws of unequal velocities and unequal diameters. 

In the first category, the shape and the size of both the screws are identical and hence 
they rotate at the same speed. In the second category', the number of lobes, velocity, etc. 
may be different for both the screws w'hereas in the third category, the screw diameters 
also will be different. In general, the cross-sections of the screws are composed of circular 
arcs. However, in the case of the third category, only one screw can be designed to have 
the cross-section composed of circular arcs and the flank curves of the other one will be 
non-circular. 

Another major classification is based on an important geometric attribute of the 
pair of screws, viz., the clearance between the screws. They are; 


i) Fully-wiped screws. 
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ii) Screws with constant clearance and 
iii) Screws with variable clearance. 

The fully-wiped screws have practically no clearance between them. These are used for 
the extrusion of highly homogeneous mass. Depending on the flow properties and 
homogeneity of the mass, one has to introduce the required clearance. In some screws, 
the clearance is designed to be variable so as to achieve the desired design specifications 
such as mass flow per revolution. 

There are other less important classifications. The screws may be (i) Single-lobed 
screws, (ii) Bi-lobed screws, (iii) Tri-lobed screws etc. based on the number of lobes and 
(i) Co-rotating and (ii) Contra-rotating based on the direction of rotation. 

The cross-section of a twin-screw has a unique shape for a given diameter, 
centerline distance and number of parallel channels because of the requirement that one 
screw must, in any position, wipe its mate. This unique shape is determined by means of 
a well-known kinematic principle. Kinematic analysis yields equations that relate the 
major design parameters and in gr^hical form, serve as the basis for design and selection 
of screw dimensions. Volumes and surface area can be calculated once the major 
parameters have been selected. 


6 ^ Design of Screw I of the Extruder 

One can think of a screw with a given cross-section as being obtained by 
uniformly twisting a bar with that cross-section. Thus the screw is also an infinite set of 
infinitesimally thin discs, each rotated slightly with respect to its neighbor. Since all pairs 
of discs are fully wiped, we have to only study one pair to determine the cross-sectional 
screw cross section. 

Let screw I rotate with an angular velocity of Oj about O, and screw II rotate with 
an angular velocity of o, about O. (Figure 6.4). Let the center distance between O, and 
O 2 be Tq. Let PQ be the tip segment of disc I. We want to know the shape of disc II that 
is wiped by the tip PQ which means that we must construct the path or trajectory of PQ 
with respect to disc II. Constructing such a trajectory becomes simpler when we add a 
motion to the system of discs I and II which brings disc II to rest. Therefore the disc II 
becomes standstill when we superimpose a rotation -©2 about Oq. As a result of this, disc 
I slows down to (©, - ©,) and also goes round Oq with an angular velocity of -© 2 . 
Therefore, the trajectory of is a circle of radius Tq with its center at Gq. Let r, and r 2 
be the tip radii and «, and n 2 be the number of lobes of screws I and II respectively. Let 
k be the center-distance ratio of any screw which is the ratio of its tip radius and the 
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Figure 6-4. Kinematic Diagram of the Twin-Screw Mechanism 
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center distance of the pair of screws and m be the velocity ratio of the pair of screws. 
Therefore, 

k, = ^ (6-4) 

^0 


*2 = 


(6-5) 


m 



rh 


( 6 - 6 ) 


6.2.1 Determination of the Flank Curve 


The line 0^0-^ rotates with an angular velocity of - 03 . Let P is any generic point 
on the tip segment of disc I and C>, be the position occupied by Of at any instant. Also 
the line rotates with an angular velocity of (©] - Oj). Since it is a radial line of the 
tip segment, it has a length of r,. Hence 


'll 

e. 


(0, - CO, 

-1 i = 1 - m 


(6-7) 


CO, 


Now the flank curve of the disc II which is the trajectory' of the point P can be expressed 
in terms of the parameter 0 as 


(6-8) 


Xp = Tq cos 62 - Tj cos Y2 = ^0 cos 02 - Tj cos [(1 - JH ) 62] 
jp = Tq sin 02 - Tj sinYo = Tq sin02 - Tj sin [( 1 -m) ©2] 

In polar coordinates and in terms of the center-distance ratios, this can be wTitten as 

(6-9) 


= ^0 [1 + ^ cos(7n02)] 


To determine the tip angle, we must know angle cp where the trajectory of P intersects the 
tip circle of disc I. At this point, r = r 2 and 0=6'. Hence, 
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or 


2 

h = '■o V ■ 2*1 cos(i«02)]^ 


( 6 - 10 ) 


02 = -^ COS"*^ 

m 


1 + 


2 *1 


( 6 - 11 ) 


Next, we must determine angle Vj/j. This is obtained by considering the triangle 
OoP,'0,'. Therefore, 


= COS 


-1 


2 2 2 
To + 


1 + *2^ 

2 To Tj 


2*2 ■ J 


(6-12) 


Therefore the flank angle (pj of disc II is 
< 1>2 = ^2 + 


= COS 


-1 


1 + *2^ - 


2k, 


+ — cos 
m 


\ + - k^ 


2 *, 


(6-13) 


6.2.2 Relation between the Tip and the Root Angles 

a-, and a. are the tip angles and Pj and are the root angles of discs I and II 
respectively. Since the tip segment of disc I wipes the root segment of disc 11, during the 
period the point of contact between the tip segment of disc 11 and the root segment of disc 
I moves from one end of the tip segment of disc II to the other, the disc I would rotate 
by (1 - m) P;. Therefore, 

«i = P2 ■ (1 " "») h 

= m^. 


(6-14) 
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Similarly, 


Pj = aj - (1 - m) 

= OT a. 


(6-15) 


Equations (6-14) and (6-15) imply that 

Pi «2 

Thus the ratio of the root to tip angles of one disc is the inverse of its mate. 


(6-16) 


The disc I has «i lobes hence it has 2ni flanks and similarly the disc II has 2«2 flanks. 
Hence, 


«i[«i Pi 24 )i] = 2 it 

Ti^aj + P 2 + 2 4 ) 2 ] = 2 7c 


(6-17) 


“2 P2 = 


27t 


ttj + Pi 


271 



- 2 m cos”^ 


2 cos 


(6-18) 


From equation (6-16) and (6-18), we can arrive at a, and p, for the given parameters. 
The following relation exists between the tip and root radii and the center distance; 
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6.2.3 Algorithm for the Design of Disc I 

The parameters that describe disc I are obtained from the following steps: 

1. Calculate a 2 and Pj from Equations (6-16) and (6-18). The ratio of the angles 

a,/p, or p 2 /a,, n,, /I 2 . ^2 ^0 wi^^ be known from the design specification. 

2. Calculate a, and P, from Equation (6-18). 

3. Calculate <{), from Equation (6-13). 

4. Define V , the position vector of a point P on the disc I for all the four arcs of 
each lobe. There will be 4n, arc segments that form the disc I. 

=1^1 yi 0 1 ] ( 6 - 20 ) 


where, 

=Pi * 

3^1 = qi + r. sina,. 

% ^ ^ «i2 

1 £ i s 4nj 


6.3 Design of Screw II of the Extruder 

6.3.1 Finding the Conjugate Segments of Disc II Corresponding to Each Segment 
of Disc I 

The profile of disc I can be obtained from the algorithm given in the preceding 
section. The twelve parameters describing the motion of the disc I and disc II are also 
known. Since the two discs are fully-wiped, the profiles of the discs are conjugate to each 
other. Hence the profile of disc II can be determined using the conjugate geometry model 
with twelve degrees of freedom proposed in Chapter 2. 

The steps involved in finding the profile of disc n are as follows; 

Define !p ; This definition is given Equation (6-19). 


1 . 
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2. Define the 12 parameters of motion: After superimposing a rotation of angular 
velocity Oj about Zq as described in Section 6.2, disc II is stationary and disc I 
rotates about Z, axis with an angular velocity of ©, - ©j and also goes round Og 
with an angular velocity of ©j. The twelve parameters of motion for this case are 
given in Table 6-1. 


Disc 1 

Disc li 

a, = Tq cost 

a, = 0 

sin t 

II 

o 

c, =0 

• II 

o 

0, =0 

e, =0 

o 

It 

-0 

II 

1 

V T 0 


3. 


Table 6-1. Parameters Describing Relative Motion between the Two Screws. 


The transformation matrix [ ] for transforming the position vector V and 

the normal vector n p j from the cutter coordinate system to the global 
coordinate system Sf^ is; 


ms Twi-Tt f 

V'*- - 



-sin (m-l)t 
0 


sill (w-1) t 


cos (w-l)t 
0 


recast Tq sin t 


0 0 ■ 
0 0 
1 0 
0 1 


( 6 - 21 ) 


Matrix j is^the time-differential of the matri5|jM 1 which converts the 
position vector P in the cutter coordinate system o, into velocity vector 
V PI in the global coordinate system Sq. 
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-(w-l)sm (w-l)cos 

-(m-l)cos (w-l)r -(in-l)siii (w-l)r 
0 0 
-r^mt recast 


0 O' 
0 0 
0 0 
0 0 


( 6 - 22 ) 


4. The transformation matrices ^ j and|^ ^ j for transforming the position 
vector from the blank coordinate system Oj the global coordinate system 
and vice-versa are respectively; 

1 0 0 
0 1 0 
0 0 1 
0 0 0 



0 

0 

0 

1 


(6-23) 



10 0 0 
0 10 0 
0 0 10 
0 0 0 1 


(6-24) 


Matrix j is the time-differential of the matrix jwhich converts the position 
vector P in the blank coordinate system into velocity vector ^ in the global 
coordinate s>’stem S^.. 



0 0 0 0 
0 0 0 0 
0 0 0 0 
0 0 0 0 


(6-25) 
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2 

The position vector P is 


?,.cos - ^,-sm (m-l)r + r.cos [«,. + (7n-l)r]l 


+ TqCOS t 


p.sin + g,C 06 (w-l)f + r-sin [«,- + (m-l)r] 


+ Tq sm r 


(6-26) 


The normal vector pj of the cutter surface is obtained as follows; 
=[-r,.sintt. Tjcos 0 Oj 

Vv =[ 0 0 10] 


(6-27) 


^ r cos « 


r sm u 


(6-28) 


7. The unit normal vector p i of the cutter surface is obtained as follows; 

I '^,,1 I = 


’"p.i = [cos « 


0 0 ] 


(6-29) 


8. The normal vector ” ^ ^ of the cutter surface in the global coordinate system is; 

= [cos [«,. + (/n-l)t] sin [u^ + 0 Oj (6-30) 


9. The relative velocity vector p.n in the global coordinate system can be now 
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or 


sin (m-2)t - — ' 


u, = tan 


-1 


cos (m-2)t + 


(m-l)Pi 


(6-35) 


1 1 . The profile Ej of screw 11 can be now represented as a parametric equation in 
either or t as given below: 


X2 = PjCos (m~l)t - 9,. sin 

+ r.cos [u. + (m-l)t] + r^cos f 
>>2 = ;j,-sin + g^cos 

+ r-sin [Uf + (m-l)t] + am t 


(6-36) 


where 


t = 



m-2 


or 


/ (w-1)^. 

sin (m-2)t - — 


u, = tan"* 


cos (m-2)t + 


^0 j 


and 


My S M. S My 

Is i s 4nj 


63.2 Finding Conjugate Segments of Disc II Corresponding to the Cusps of Disc I 
When the profile of disc I is a closed curve and has C'-continuity throughout, then 
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the above eleven steps will yield the profile of disc II completely. However, since the disc 
I of the extruder is generally composed of circular segments, it does not satisfy this 
condition and hence there will be cusps. Cusp is a point that is common to both the 
preceding and the following curves which have different normals at that point. Hence, at 
the cuspoidal point, the normal abruptly changes. The procedure described in the previous 
section maps each point of profile of disc I and hence it considers each cusp for the two 
end values of its normal only and not for the values in-between. Hence the presence of 
a cusp is observed as a discontinuity in the Ej profile obtained at Step 11. Mapping of 
each cusp will have to be done for each value of its normal. In order to do this, one needs 
to detect the presence of the cusps and then perform their mapping. Figure 6-5 shows the 
cross-section of Screw I and the cross-section of Screw II produced using the conjugate 
geometry formulation. In this figure, the solid lines are conjugate lines to the segments 
of the profile I and the dotted lines were generated by the cuspoidal points of this profile. 


63.2.1 Detection of Cusps 


In the definition of disc I as given in Equation (6-19), u, which is the absolute 
angle of any point on an arc with respect to its center (p^, g,). Let u be a parameter wdiich 
is the absolute angle of any point P on the profile of disc I with respect to O, (Figure 6- 
6a). The relation between u and is given by: 




( 6 - 37 ) 


The plot of u versus u, is shown in Figure 6-6b for the profile of disc I shown in 
Figure 6-6a. The data pertaining to the eight segments of disc I of this profile is given in 
Table 6-2. 


/ 

A 

Q, 

r, 

A, 

Ai 

1 

-99.804 

-6.262 

150.000 


0.78540 

2 

0.000 

0.000 

100.000 

1.50811 

1.63349 

3 

99.804 

-6.262 

150.000 

2.35619 

j.v foyv 

Hi 

0.000 

0.000 

50.000 

3.07890 

3.20428 

1 5 

99.804 

6.262 

150.000 

3.20428 

3.92699 
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6 

0.000 

0.000 

100.000 

4.64970 

4.77508 

7 

-99.804 

6.262 

150.000 

5.49779 

6.22049 

8 

0.000 

0.000 

50.000 

6.22049 

6.34587 


Table 6-2. Biparametric Definition of the Profile Shown in Figure 6-6a. 


It can be observed that there are four cusps and the graph shows the corresponding four 
discontinuities. This property can be used to collect the cusps as given below; 

Let w^i and be the starting and ending values of u, for the /th segment 
of the profile. If m ,2 ^ then the end point of the fth segment is a 
cuspoidal point. 


63.2.2 Mapping of Cusps 

Each cu^ Cj identified earlier can be mapped onto the Sj coordinate system using 
the following proc^ure; 

1 . At the cusp Cj which is Ae end point of segmwit / and starting point of segment 

(i+l), Ae normal vector changes from Aat of Ae rA segment to Aat of (f+l)A 
segment. The interval in t, say, and tj 2 , in vAich Ais change takes place is 
obtained by substiAting Ae values of (u,^, p„ q, and r, ) and j, and 

in Equation (6-34). 

2. The cusp can be assumed to be an arc of zero raAus. Hence, use Equation (6-35) 
for mapping Ae cusp for Ae interval to tj, wiA {p^ q^ as Ae cuspoidal point 
in 5,, r, = 0 and u, = 0 

6.4 Illustrative Examples 

As mentioned earlier in Section 6.1.2, Ae pair of extruder screws are classified as 
symmetrical and unsymmetrical. The unsymmetrical screw pairs are furAer classified as 
equal diameter screw pairs and unequal diameter screw pairs. In Ais sections, several 
examples of symmetrical screw pairs and unsymmetrical screw pairs of boA types are 
considered to illustrate Ae application of the conjugate geometry model for the design of 
Aese screws. 
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6.4.1 Design of a Twin-Screw Extruder of Symmetrical Screws 
Since the screws of a symmetrical screw pair are identical. 


r = r = r 

= a 

p, = p, - p 


(6-J8) 


As they rotate with equal speed. 


= n, - n 
m = 1 


(6-39) 


Equations (6-14), (6-15), (6-38) and (6-39) can be reduced as 

a, = = a, = 3, = a (say) 

From the above. Equation (6-18) can be simplified as 


/ ^ 


a = - - 2 cos 
n 




(6-40) 


(6-41) 


where a > 0. 

Therefore, the design procedure for this class of screws can be stated as follows: 

From the given design specifications, determine the number of lobes n, the tip 
radius r, the tip angle 23 and the center distance using Equations (6-41) and the 
following two equations which are applicable to all kinds of screw pairs : 

/-+/-,= (6-42) 

Examples of five types of symmetrical screw pairs are explained in this section. 


6.4.1. 1 Design of a Twin-Screw Extruder of Mono-Lobed Symmetrical Screws 

r 


Given : 


The tip radius 
root radius 


100 mm. 
25 mm. 
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Figure 6-7. A Mono-Lobed Screw I. 


The value of for this case as obtained from Equation (6-42) will be 125 mm. 
Therefore, the values of the tip and root angles a = p are calculated as 77.3648° from 
Equation (6-41). Table 6-3 describes the parametric equation of the profile of screw I as 
given in Equation (6-20). The plot of this profile is shown in Figure 6-7 and the generated 
the cross-sectional profile of screw 11 using the algorithm given in Section 6.3 is shown 
in Figure 6-8. 


















Figure 6-8. Generation of the Profile of Screw II for a Mono-Lobed Screw I 
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3 

62.500 

78.063 

125.000 

3.14159 

4.03726 

n 

0.000 

0.000 

25.000 

4.03726 

5.38752 


Table 6-3. Biparametric Definition of the Profile Shown in Figure 6-7. 


6.4.1.2 Design of a Twin-Screw Extruder of Bi-Lobed Symmetrical Screws 



Figure 6-9. A Bi-Lobed Screw I. 


Given : The tip radius 

root radius 
Number of lobes 


r =100 mm. 

= 50 mm. 

n =2. 


The value of for this case as obtained from Equation (6-42) wiE be 150 mm. 
Therefore, the values of the tip and root angles a = 3 are calculated as 7.18076° from 
Equation (6-41). Table 6-4 describes the parametric equation of the profile of screw I as 








Figure 6-10. Generation of the Profile of Screw II for a Bi-Lobed Screw I 
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given in Equation (6-20). The plot of this profile is shown in Figure 6-9 and the generated 
the cross-sectional profile of screw II using the algorithm given in Section 6.3 is shown 
in Figure 6-10. 


/ 

A 

9, 

6 

Ar 

As 

1 

-99.804 

-6.262 

150.000 

0.06269 

0.78540 

2 

0.000 

0.000 

100.000 

1.50811 ■ 

1.63349 

3 

99.804 

-6.262 

150.000 

2.35619 

3.07890 

B 

0.000 

0.000 

50.000 

3.07890 

3.20428 

5 

99.804 

6.262 

150.000 

3.20428 

3.92699 


0.000 

0.000 

100.000 

4.64970 

4.77508 

7 

-99.804 

6.262 

150.000 

5.49779 

6.22049 

g 

0.000 

0.000 

50.000 

6.22049 

6.34587 


Table 6-4. Biparametric Definition of the Profile Shown in F^ure 6-9. 


6.4.13 

Design of a Twin-Screw Extruder with Bi-Lobed Symmetrical Screws 


with Sharp Tips 



Given : 

The tip radius 

r 

= 100 mm. 


Tip angle 

a 

= 0° 


Number of lobes 

n 

= 2. 


The value of for this case as obtained from Equation (6-41) will be 141.414 
mm. Therefore, the values of the root diameter of the screws is calculated as 41.414 
mm. from Equation (6-42). Table 6-5 describes the parametric equation of the profile of 
screw I as given in Equation (6-20). The plot of this profile is shown in Figure 6-11 and 
the generated the cross-sectional profile of screw II using the algorithm given in Section 
6.3 is shown in Figure 6-12. 
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/ 

A 

<7/ 

ri 


U,2 

1 

-100.000 

0.000 

141.414 

0.00000 

0.78540 

2 

0.000 

0.000 

141.414 

2.35619 

3.14159 

3 

100.000 

0.000 

141.414 

3.14159 

3.92699 

B 

0.000 

0.000 

141.414 

5.49779 

6.28319 


Table 6-5. Biparametric Definition of the Profile Shown in Figure 6-11 



Fioure 6-11. A Bi-Lobed Screw I with a = 0. 

6-4.1.4 Design of a Twin-Screw Extruder with Tri-Lobed Symmetrical Screws 

The tip radius r 

Tip angle a 


Given : 


= 100 ram. 
= 20 ° 
































Figure 6-12. Generation of the Profile of Screw II for a Bi-Lobed Screw I 
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Figure 6-13. A Tii-Lobed Screw I. 


Number of lobes n =3. 

The value of Tq for this case as obtained from Equation (6-41) will be 187.939 
mm. Therefore, the values of the root diameter of the screws is calculated as 87.939 
mm. from Equation (6-42). Table 6-6 describes the parametric equation of the profile of 
screw I as given in Equation (6-20). The plot of this profile is shown in Figure 6-13 and 
the generated the cross-sectional profile of screw n using the algorithm given in Section 
6.3 is shown in Figure 6-14. 


B 

A 

9, 

r, 


U.2 

1 

-76.604 

-64.279 

187.939 

0.69813 

1.04720 

2 

0.000 

0.000 

100.000 

1.39626 

1.74533 

3 

76.604 

-64.279 

187.939 

2.09440 ■ 

2.44346 























Figure 6-14. Generation of the Profile of Screw II for a Tri-Lobed Screw I 
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n 

0.000 

0.000 

87.939 

2.44346 

2.79253 

5 

93.969 

-34.202 

187.939 

2.79253 

3.14159 

6 

0.000 

0.000 

100.000 

3.49066 

3.83972 

7 

17.365 

98.481 

187.939 

4.18879 

4.53786 

8 

0.000 

0.000 

87.939 

1 

4.53786 

4.88692 

9 

-17.365 

98.481 

187.939 

4.88692 

5.23599 

10 

0.000 

0.000 

100.000 

5.58505 

5.93412 

11 

-93.969 

-34.202 

187.939 

0.00000 

0.34907 

12 

0.000 

0.000 

87.939 

0.34907 

0.69813 


Table 6-6. Biparametric Definition of the Profile Shown in Figure 6-13. 


6.4.1^ Design of a Twin-Screw Extruder with Quadri-Lobed Symmetrical 

Screws 

Given ; The tip radius r =100 mm. 

Tip angle a =10° 

Number of lobes n =4. 

The value of for this case as obtained from Equation (6-41) will be 190.743 
mm. Therefore, the vedues of the root diameter of the screws is calculated as 90.743 
mm. from Equation (6-42). Table 6-7 describes the parametric equation of the profile of 
screw I as given in Equation (6-20). The plot of this profile is shown in Figure 6-15 and 
the generated the cross-sectional profile of screw II using the algorithm given in Section 
6.3 is shown in Figure 6-16. 
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Figure 6-15. A Quadri-Lobed Screw I. 


‘ P. Pi n iiii U-a 

1 -64,281 -76,603 190.743 0.87264 1.17808 

2 0.000 0.000 100.000 1.48352 1.65805 

3 64,281 -76.603 190.743 1.96349 2.26893 

4 0.000 0.000 90.743 2,26893 2.44344 

5 76.603 -64.281 190.743 2.44344 2.74887 

f OOQO 0.000 100.000 3.05431 3.22885 

2 76.603 64.281 190.743 3.53429 3.83972 

^ QQQQ 0000 90.743 3.83972 4.01423 
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Figure 6-16. Generation of the Profile of Screw II for a Quadri-Lobed Screw 
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9 

64.281 

76.603 

190.743 

4.01423 

4.31967 

10 

0.000 

0.000 

100.000 

4.62514 ■ 

4.79964 

11 

-64.281 

76.603 

190.743 

5.10508 

5.41052 

12 

0.000 

0.000 

90.743 

5.41052 

5.58503 

13 

-76.603 

64.281 

190.743 

5.58503 

5.89047 

14 

0.000 

0.000 

100.000 

6.19591 

3.22884 

15 

-76.603 

-64.281 

190.743 

0.39269 

0.69813 

16 

0.000 

0.000 

90.743 

0.69813 

0.87264 


Table 6-7. Biparametric Definition of the Profile Shown in Figure 6-15. 

From the above five examples, it can be observed that the angular offset between 
the two screws of a symmetrical screw pair is 

0 if n is odd and 

it/n if n is even. 


6-4.^ Design of a Twin-Screw Extruder of Unsymmetrical Screws of Equal 
Diameters 


Since the diameters of the screws are equal. 


K = /• = /- 

'2 ' 


( 6 - 43 ) 


As the shapes of these screws will be different, they rotate with imequal speed. 


m - 


ffi, ^ 


(6-44) 


Therefore, Equations (6-14), (6-15), (6-43) and (6-44) yield 
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3 , 

2 


{ \ 


n 

n. 


- (m +1) cos"’ 



(6-45) 


where a, > 0 and P, > 0. 

Therefore, the design procedure for this class of screws can be stated as follows: 

From the given design specifications, determine the number of lobes n,, the tip 
radius r, the tip and root angles ttj, P, and the center distance using Equation (6-42) 
and Equation (6-45). Example of a screw pair of this type is illustrated in this section. 


6,4.2.1 

Design of a Twin-Screw Extruder of Unsymmefrical Screws of Equal 


Diameters - Bi-Lobed Screw I 


Given : 

The tip radius 

r 

= 100 mm. 


root radius 

rr 

= 55 mm. 


Number of lobes 


= 2. 


Number of lobes 

«2 

= 1. 

= Pi 


The value of for this case as obtained from Equation (6-42) will be 155.000 
mm. and m = 0.5. From Equation (6-45), it is calculated that ttj + pj = 62.415° and hence 
ttj = Pj = 31.208°. Table 6-8 describes the parametric equation of the profile of screw I 
in the form given in Equation (6-19). The plot of this profile is shown in Figure 6-17 and 
the generated the cross-sectional profile of screw n using the algorithm given in Section 
6.3 is shown in Figure 6-18. 


B 

A 

A 

r, 

A, 


B 

-99.757 

6.965 

155.000 

6.23824 

0.61438 

2 

0.000 

0.000 

100.000 

1.29846 

1.84313 

3 

99.757 

6.965 

155.000 

2.52722 

3.18654 

B 

99.757 

-6.965 

155.000 

3.09664 

3.75597 

5 

0.000 

0.000 

100.000 

1 

4.44005 

4.98473 
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6 

-99.757 

-6.965 

155.000 

5.66881 

0.04495 


Table 6-8. Biparametric Definition of the Profile Shown in Figure 6-17. 



Figure 6-17. A Bi-Lobed Unsymmetrical Screw I. 


6.43 Design of a Twin-Screw Extruder of Unsymmetrical Screws of Unequal 

Diameters and Unequal Speeds 

These are more general case of screw pairs whose kinematic behavior is controlled 
by the Equations derived in Sections 6.2.1 and 6.2.2. The design procedure for this can 
be summarized as below; 

From the given design specifications, determine the number of lobes the tip 
radii r^, rj, the tip and root angles a,, aj, ^2 carter distance using 

Equations (6-14), (6-15, (6-16), (6-18) and (6-41). This has been explained with the help 












Figure 6-18. Generation of the Profile of Mono-Lobed Screw II for a Bl-Lobed Screw I 
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of three examples. 


6.4.3.1 

Design of a Twin-Screw 

Extruder of Unsymmetrical Screws of 


Unequal Diameters and Unequal Speeds - Mono-Lobed Screw I and 


Velocity Ratio m = 2 



Given : 

The tip radius of Screw I 


= 100 mm. 


The tip radius of Screw I 


= 66.667 mm. 


Center Distance 

^0 

= 200 mm. 


Number of lobes 

«I 

= 1. 


Number of lobes 

«2 

= 2. 


ot] P, 



Figure 6-19. A Mono-Lobed Unsyrtwnetrical Screw I for m * 2. 


The tip and root radii of screw II and r,; for this case are respectively 133.333 
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Figure 6-20. Generation of the Profile of Bi-Lobed Screw H for a Mono-Lobed Screw I 
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mm. and 100.000 mm. as obtained from Equation (6-19). Hence, m = 2 and A, = 0.5 and 
Aj = 0.66667. From Equation (6-45), it is calculated that ttj + p, = 181.780” and hence 

а, = P, = 90.895”. Table 6-9 describes the parametric equation of the profile of screw I 
in the form given in Equation (6-20). The plot of this profile is shown in Figure 6-19 and 
the generated the cross-sectional profile of screw II using the algorithm given in Section 

б. 3 is shown in Figure 6-20. 


B 

A 

A 

n 

■ 

U.2 

1 

0.000 

0.000 

100.000 

0.77759 

2.36401 

2 

127.690 

49.697 

200.000 

3.03912 

3.64491 

3 

0.000 

0.000 

66.667 

3.91918 

5.50560 

D 

-127.000 

49.697 

200.000 

5.77986 

6.38566 


Table 6-9. Biparametric Definition of the Profile Shown in Figure 6-19. 


6.43.2 


Given : 


Design of a Twin-Screw Extruder of Unsymmetrical Screws of 
Unequal Diameters and Unequal Speeds - Mono-Lobed Screw I and 
Velocity Ratio m = 3 


The tip radius of Screw I 
The tip radius of Screw I 
Center Distance 
Number of lobes 
Number of lobes 
Number of lobes 
Number of lobes 


=100 mm. 

= 26.667 mm. 
Tq • = 266.667 mm. 
n, = 1. 

/ij 3. 

a, = 16.757”. 

p, = 16.75r. 


The given data fully describe the specifications for screw I and the relative 
motions. Table 6-10 describes the parametric equation of the profile of screw I in the 
form given in Equation (6-20). The plot of this profile is shown in Figure 6-21 and die 
generated the cross-sectional profile of screw 11 using the algorithm given in Section 6.3 
is shown in Figure 6-22. 
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F^ure 6-21 . A Mono-Lobed Unsymmetrical Screw I for m = 3. 


/ 

P, 

% 


Pi 


1 

0.000 

0.000 

100.000 

1.42456 

1.71703 

2 

20.343 

38.262 

70.000 

2.09298 

4.22348 

3 

0.000 

0.000 

26.667 


5.20130 

4 

hbh 



5.20130 

7.33180 


Table 6-10. Biparametric Definition of the Profile Shown in Figure 6-21. 






























•e 6-22. 
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Design of a Twin-Screw Extruder of Unsymmetrical Screws of 
Unequal Diameters and Unequal Speeds - Mono-Lobed Screw I and 
Velocity Ratio m = 4 



Figure 6-23. A Mono-Lobed Unsymmetrical Screw I for m = 4. 


Given : The tip radius of Screw I 

The tip radius of Screw I r , , 
Center Distance Tq 

Number of lobes n, 

Number of lobes 
Number of lobes a, 

Number of lobes P, 


= 100 mm. 

= 20.833 mm. 
= 208.333 mm 
= 1 . 

= 4. 

= 0.34172°. 

= 0.34172°. 


Tne given data imiy descriDe me ^lecincahons for screw i and me re’iMve 
motions. Table 6-11 describes the parametric equation of die profile of screw I in the 
form given in Equation (6-20). The plot of this profile is shown in Figure 6-23 and the 




■e 6-2' 
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generated the cross-sectional profile of screw II using the algorithm given in Section 6.3 
is shown in Figure 6-24. 


/ 

P, 

Qi 

n 



1 

0.000 

0.000 

100.000 

1.56781 

1.57378 

2 

12.287 

39.290 

62.000 

1.77520 

4.40927 

3 

0.000 

0.000 

20.833 

4.40927 

5.01557 

4 

-12.287 

39.290 

62.000 

5.01557 

7.64958 


Table 6-11. Biparametric Definition of the Profile Shown in Figure 6-23. 






















Chapter 7 


MANUFACTURE OF EXTRUDER SCREW SURFACES 


7.1 Introduction 

The design of cutters for the manufacture of helicoidal surfaces such as the ones 
used in the extruder screws is an important area of research. These surfaces are 
manufactured mostly by milling processes using form cutters of the end mill type, the side 
mill type or the side-and-face mill type (i.e., disc type). The problem of design of cutters 
for machining helicoidal surfaces is not yet satisfactorily solved (Saybold, 1992). In this 
chapter, a methodology is proposed to solve this problem. Using tiiis method, one can 
determine die two-dimensional profile of the cutter from the given cross-sectional profile 
of a helicoidal surface. This two-dimensional profile is similar to the profile of the generic 
cutter described in Chapter 4. The proposed methodology can also be used to determine 
the shape of the cross-sectional profile and hence the helicoidal surface that will be 
obtained using a given cutter. The problems of the design of both the end mill type and 
disc type cutters have been solved. The algorithm proposed for the design of the end mill 
type cutters is also valid for the design of form turning tools if it is required to realize the 
helicoidal surface by a turning process. 


7.2 Methodology 

Let ^8 (w) be the cross-sectional profile of the helicoidal surface which can be 
considered to be a set of C'-continuous arc segments as given below : 

= :v, 0 1] 


(7.1) 
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where. 


^2 =/?,. + r.cosUf 
>2 = Qi + sinwj 

tty i tt,. i ttj2 

1 i i i 4n 

Let Sq (Oq- Xq, Yq, Zq) be a fixed coordinate frame and let $2 (Oj- Xj, 72. ^ 2 ) be a 
coordinate frame to which the profile S (“) is initially attached in the X 2 - Y 2 plane as 
shown in Figure 7-1. If S (^) is considered to be a generatrix of the helicoidal surface 
then by applying the following helical sweep to the generatrix, we can obtain the 
helicoidal surface P (tt, v) .Note that p is the pitch of the helix. 


V (tt, V) = ^g{u) 


cos V sin V 0 

-sinv cos V 0 

0 0 1 

0 0 ^ 

2 It 


0 

0 

Q 

1 


(7-2) 


On the other hand, we can consider a plane that cuts the helicoidal surface normal 

to Ae helix (Figure 7-2). If Ais plane is moved along Ae axis of Ae helix, at an 
appropriate position, Ae cross-section of Ae helicoidal surface captured on Ais plane will 
be Ae cutter profile. The same helicoidal surface as Ae one given by Equation (7-2) will 
be obtained if Ais profile of Ae cutter, say ^8 (tt) , is given exactly Ae same helical 
sweep. 


°P (tt, v) = ‘^^(tt) 


cosv sm V 0 0 

-^v cos V 0 0 

0 0 10 

pv . 


0 


0 


2 tt 


(7-3) 


From Ais, it is evident Aat every point of Ae profile ^8 00 can be mapped onto Ae 
profile ^8 (tt) and vice-versa. Therefore, if a helicoidal surface is known, one can obtain 
the surface of Ae cutter Aat can be used to machine it and one can also determine Ae 
geometry of a helicoidal surface Aat will be produced by a giv^ cutter. 
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7.2.1 Definition of the 52and Coordinate Frames 

As shown in Figures 7-1 and 7-2, Sq (Oq- Xq, Yq, Zq) is a fixed frame of reference 
attached to the helicoidal surface. The cross-sectional profile 
^8 0*) initially coincides with the Xq - Tq plane. It is given a helical sweep described 
by the sweep parameter v and the pitch of the helix p in order to obtain the helicoidal 
surface. Let a be the lead angle of the helix. If r is the tip radius of the cross-sectional 
profile ^8 (**), then 


a = tan 


-1 P 


2 Tcr 


(7-4) 


Consider another coordinate frame S, (O^- X^ Y„ Z,) which is obtained by 
applying the following two transformations to the system : 


i) A rotation about Yq axis by an angle of ^ ^ 


w * 

ii) A translation along Zq axis by a distance of — — 

2:c 


where 


Therefore any point 
following relation : 


°P = 



V* 

7t 

n 


the i 

space 

can 

be mapped 


■ -sino 

0 

-cos a 0 


0 

1 

0 0 


cos o 

0 

-sin 0 0 


0 

0 

pv j 

2it 


(7-5) 


Sq space using the 


(7-6) 


The initial shape of the blank on which die helicoidal surface is machined will be 
a cylindrical block of radius r which is the tip radius of the cross-sectional profile 

^8 («) . The shape of the cutter profile will be decided by the root and flank arcs of the 
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cross-sectional profile. When v — v , the mid-point of the root arc of the cross-sectional 
profile will lie in the Y, axis. Hence, the profile of the cutter will be fully obtained in the 
plane by finding the intersection curve of this plane with the helicoidal surface. 
It is enough to consider the range of v between 

V* - — and V ’ + 
n n 

7.2.2 Determination of the Profile of the Cutter (**) Lying on a Plane Parallel 
to and at a Distance of w from the Plane from the given cross-sectional 
profile \ («) 

Ths >"0 Zq coordinates of the helicoidal surface in terms of the cross- 
sectional profile ^8 («) can be written from Equation (7-2) as 

= JCj cos v - sin v 
= Xj sin V + ^2 cos V 

= 

2 % 


•'0 

3^0 


From Equation (7-6), one can write the equation of the plane parallel to and at a distance 
of w from - 7^ as 


Xg = -x^ sin a + w cos a 

3’o =7, (7-8) 

pv* 

z - £ — - X cos a - w sin a 
® 271 " 

The profile of the cutter which is the intersection curve between the helicoidal surface and 
the above plane will be obtained by solving the Equations (7-7) and (7-8) for x^ and y^ 
as follows : 

From the expressions for x,, and y^ of the above two equations, one can solve for 
x^ and y^ as 
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X = 


w 


cos V - sm v 


tan a sm 0 

= Xj sin V + >>2 cos V 


(7-9) 


Substituting these values of x^ and 3;^ in the Zq expressions, we get the following 
equation in v. 


p tan o 

y, sin V - cos v + v 

^ " 2 7t 




p tano 
27C 



0 (7-10) 


The solution of the above equation has to be found numerically to find the value 
of V for a given set of values. Once v is thus known, the corresponding point of 
the cutter profile (x^, j,) can be calculated. It is observed that the convergence of 
the solution is assured if v’ is given as the starting value for v in the iterative 
process. 


7.2.3 Determination of the Cross-sectional Profile of the Helicoidal Surface E (“) 
from the given profile of the Cutter 00 at tv = 0 

By comparing the expressions of Equations (7-8) and (7-9), we can find the value 
of V as 


V 


= V 


2n 

P 


X cosa 

C 


(7-11) 


Again comparing the expressions for Xq and yo of Equations (7-8) and (7-9), we can 
obtain the cross-sectional profile of the screw 00 as 


X, = y sin V - x sin 0 cos v 
= y^ cos V + x^ sin 0 sin v 


(7-12) 


V = V* 


2k 


X cos 


p 


a 


where 
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7.3 Illustrative Examples for the Design of End Milling Cutters 

In this section, the methodology described in Section 7.2.2 will be used to 
determine the cutter profile at plane (i.e at w = 0). An end milling cutter to 

machine the helicoidal surface will be a surface of revolution of this profile about the f, 
axis. The method of obtaining the cutter profile for this case will be illustrated by means 
of examples for the machining of extruder screws. 


7.3.1 Design of an End Milling Cutter for the Machining of a Mono-Lobed 
Extruder Screw 

The cross-sectional profile of a mono-lobed extruder screw is similar to the one 
shown in Figure 6-7. It has a tip radius of 1 00 mm. and a root radius of 66.667 ram. As 
mentioned earlier, only the root and flank arcs influence the design and the tip arc is 
realized in the preparation of the blank. The dimensions of one root arc and the two flank 
arcs adjacent to it are given in Table 7-1. 


B 

A 

A 

ri 

A, 

U,2 


113.855 

69.389 

200.000 

3.22844 

3.68893 

Bl 

0.000 

0.000 

66.667 

3.68893 

5.73585 

1 

-113.855 

69.389 

200.000 

5.73585 

6.19634 


Table 7-1 . Biparametric Definition of the Root and Adjacent Flank Arcs of 
a Mono-Lobed Extruder Screw. 


Using the method described above, the profile of the cutter was calculated for 
pitch p = lS mm .,/1 = 150 mm. and p = 225 mm. and the plots of these three profiles are 

shown in Figure 7-3. The coordinates of the cutter profile is listed in Table 7-2. 





X. 


Xf 

V. 

X- 

V. 

-25.105 

99.956 

-49.179 

99.346 

-71.330 

97.066 





























Chapter 7. Manufacture of Extruder Screw Surfaces 


199 


-24.174 

94.623 

-47.249 

94.029 

-68.296 

91.824 

-23.169 

89.549 

-45.177 

88.975 

-65.067 

86.864 

-22.083 

84.792 

-42.951 

84.244 

-61.633 

82.248 

-20.907 

80.417 

-40.563 

79.902 

-57.988 

78.044 

-19.639 

76.502 

-38.006 

76.027 

-54.135 

74.330 

-18.274 

73.129 

-35.283 

72.701 

-50.088 

71.186 

-16.819 

70.383 

-32.408 

70.008 

-45.874 

68.692 

-15.289 

68.346 

-29.410 

68.028 

-41.536 

66.920 

-13.700 

67.086 

-26.326 

66.827 

. -37.130 

65.928 

-12.082 

66.651 

-23.208 

66.449 

-32.721 

65.748 

-12.046 

66.651 

-23.140 

66.450 

-32.625 

65.753 

-9.937 

66.657 

-18.512 

66.530 

-26.107 

66.083 

-7.227 

66.661 

-13.884 

66.589 

-19.584 

66.339 

: -4.818 

66.664 

-9.257 

66.632 

-13.058 

66.521 

-2.410 

66.666 

-4.628 

66.658 

-6.530 

66.630 

0.001 

66.667 

0.000 

66.667 

0.000 

66.667 

2.410 

66.666 

4.628 

66.658 

imiii^y 

66.630 

4.818 

66.664 

9.257 

66.632 

13.058 

66.521 

7.227 

66.661 

13.884 

66.589 

19.584 

66.339 

9.637 

66.657 

18.512 

66.530 

26.107 

66.083 

12.046 

66.651 

23.140 

66.450 

32.625 

65.753 

12.082 

66.651 

23.208 

66.449 

32.721 


13.700 

67.086 

26.326 

66.830 

37.130 

65.928 

15.289 

68.346 

29.410 

68.028 

41.536 

66.920 

16.819 

70.383 

32.408 

70.008 

45.874 

68.692 

18.274 

73.129 

35.283 

72.701 

50.088 

71.186 
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19.639 

76.502 

38.006 

76.027 

54.135 

74.330 

20.909 

80.417 

40.563 

79.902 

57.988 

78.044 

22.082 

84.792 

42.951 

84.244 

61.633 

82.248 

23.169 

89.549 

45.177 

88.975 

65.067 

86.864 

24.174 

94.623 

47.249 

94.029 

68.296 

91.824 

25.105 

99.956 

49.179 

99.346 

71.330 

97.066 


Table 7-2 Coordinates of the Profile of the Cutter to Machine the Mono-Lobed Extruder Screw 

Since there is only one flute on this screw, the full helicoidal surface will be obtained 
in one steps. 

7.3.2 Design of an End Milling Cutter for the Machining of a Bi-Lobed Extruder 
Screw 

The cross-sectional profile of a bi-lobed extruder screw is similar to the one shown in 
Figure 6-9. It has a tip radius of 100 mm. and a root radius of 75 mm. As mentioned 
earlier, only the root and flank arcs influence the design and the tip arc is realized in the 
pr^aration of die blank. The dimensions of one root arc and the two flank arcs adjacent 
to it are given in Table 7-3. 


/ 

A 

<7, 

6 

U,i 


7 

-92.388 

38.268 

175.000 

5.3851-3 

5.89049 

8 

0.000 

0.000 

75.000 

5.89049 

6.67588 

Hi 

-92.388 

-38.268 

175.000 

6.67588 

718124 


Table 7-3. Biparametric Definition of the Root and Adjacwit Flank Arcs of 
a Bi-Lobed Extruder Screw. 

Using the method described above, the profile of the cutter was calculate! for pitch p=15 
mm.,/?=150 mm. and /?=225 mm. and the plots of these three profiles are ^own in Figure 
7-4. The coordinates of the cutter profile is listed in Table 7-4. 
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225 mm 


- 17.115 

99.979 

- 36.530 

99.639 

- 58.194 

93.799 

- 16.154 

95.809 

- 34.547 

95.492 

- 54.843 

90.026 

- 15.132 

91.888 

- 32.430 

91.596 

- 51.221 

86.568 

- 14.048 

88.258 

- 30.169 

87.995 

- 47.307 

83.470 

- 12.898 

84.966 

- 27.756 

84.735 

- 43.089 

80.780 

- 11.682 

82.064 

- 25.188 

81.867 

- 38.562 

78.540 

- 10.402 

79.600 

- 22.467 

79.438 

- 33.745 

76.792 

- 9.063 

77.623 

- 19.605 

77.497 

- 28.673 

75.566 

- 7.676 

76.176 

- 16.623 

76.084 

- 23.406 

74.883 

- 6.254 

75.294 

- 13.552 

75.231 

- 18.021 

74.754 

- 4.813 

74.998 

- 10.432 

74.961 

- 18.022 

74.754 

- 4.812 

74.998 

- 10.431 

74.961 

mm 

74.842 

- 3.849 

74.999 

- 8.345 

74.975 

- 10.812 

74.911 

- 2.886 

74.999 

- 6.259 

74.986 

- 7.208 

74.961 

- 1.925 

75.000 

- 4.172 

74.994 

- 3.604 

74.990 

- 0.963 

75.000 

- 2.086 

74.998 

0.000 

75.000 

0.000 

75.000 

0.000 

75.000 ' 


74.990 

0.963 

75.000 

2.086 

74.998 

7.207 

74.961 

1.925 

75.000 

4.172 

74.994 

10.811 

74.911 

2.886 

74.999 

6.258 

74.986 

14.416 

74.842 

3.849 

74.999 

8.345 

74.975 

18.022 

74.854 

4.812 

74.998 

10.431 

74.961 

18.022 

74.754 

4.813 

74.998 

10.431 

74.961 

23.406 

74.883 

6.254 

75.294 

13.553 

75.231 

28.674 

75.566 
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7.676 

76.176 

16.624 

76.084 

33.745 

76.792 

9.063 

77.623 

19.605 

77.497 

38.563 

78.540 

10.402 

79.600 

22.467 

79.438 

43.088 

80.780 

11.681 

82.064 

25.187 

81.867 

47.308 

83.470 

12.898 

84.966 

27.756 

84.735 

51.221 

86.568 

14.048 

88.258 

30.169 

87.995 

54.842 

90.026 

15.133 

91.888 

32.430 

91.596 

58.193 

93.799 

16.154 

95.809 

34.547 

95.492 



17.114 

99.979 

36.530 

99.639 




Table 7-4 Coordinates of the Profile of the Cutter to Machine the Bi-Lobed Extruder Screw 


Since there are two flutes for this screw, the full helicoidal surface will be 
obtained in two steps. The blank will be indexed by it radian while machining the second 
flute. 


733 Design of an End Milling Cutter for the Machining of a.Tri-Lobed Extruder 
Screw 

The cross-sectional profile of a bi-lobed extruder screw is similar to the one shown 
in Figure 6-11. It has a tip radius of 100 mm. and a root radius of 87.939 mm. As 
mentioned earlier, only the root and flank arcs influence the design and the tip arc is 
reali 2 Bd in the preparation of Ae blank. The dimensions of one root arc and Ae two flank 
arcs adjacent to it are given in Table 7-5. 


i 

— 

A 

A 

n 

A, 


7 

-93.969 

-34.202 

187.939 

0.00000 

0.34907 

8 

0.000 

0.000 

87.939 

0.34907 

0.69813 

1 

-76.604 

-64.279 

187.939 

0.69813 

1.04720 


Table 7-5. BiparameAc Defiratkm of Ae Root mtd Adjacent Flank Arcs of 
a Tri-Lobed Extruder Screw. 
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Using the method described above, the profile of the cutter was calculated for 
pitch p=15 mm., ;7=150 mm. and p=225 mm. and the plots of these three profiles are 
shovwi in Figure 7-5. The coordinates of the cutter profile is listed in Table 7-6. 






V, 

Xf 

V. 

X, 

Vf 

-10.643 

99.993 

-22.714 

99.862 

-38.091 

99.173 

-9.882 

97.835 

-21.112 

97.719 

-35.471 

97.108 

-9.096 

95.855 

-19.453 

95.755 

-32.744 

95.224 

-8.287 

94.068 

-17.739 

93.983 

-29.910 

93.532 

-7.456 

92.487 

-15.974 

92.417 

-26.975 

92.043 

-6.603 

91.125 

-14.159 

91.069 


90.770 

-5.734 

89.993 

-12.301 

89.950 

-20.826 

89.721 

-4.848 

89.100 

-10.406 

89.070 

-17.634 

88.904 

-3.949 

88.457 

-8.481 

88.436 

-14.383 

88.325 

-3.044 

88.069 

-6.537 

88.056 

-11.089 

87.990 

-2.132 

87.939 

-4.581 

87.933 

-7.771 

87.900 

-2.132 

87.939 

-4.581 

87.933 

-7.772 

87.900 

-1.705 

87.939 

-3.665 

87.935 

-6.217 

87.914 

-1.280 

87.939 

-2.749 

87.937 

-4.663 

87.925 

-0.853 

87.939 

-1.833 

87.939 

-3.109 

87.933 

-0.426 

87.939 

-0.916 

87.939 

-1.554 

87.937 

0.001 

87.939 

0.000 

87.939 

0.000 

87.939 

0.426 

87.939 

0.916 

87.939 

1.554 

87.937 

0.853 

87.939 

1.832 

87.938 

3.109 

j 

87.933 

1.280 

87.939 

2.748 

87.937 

4.663 

87.925 

1.706 

87.939 

3.664 

87.935 

6.218 

87.914 
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2.133 

87.939 

4.580 

87.933 

7.772 

8' 

2.133 

87.939 

4.580 

87.933 

7.772 

87 

3.043 

88.069 

6.536 

88.056 

11.089 

87 

3.950 

88.457 

8.482 

88.436 

■B 

88 

4.848 

89.100 

10.405 

89.069 

17.634 

88 

5.734 

89.992 

12.301 

89.950 

20.826 

89 

6.604 

91.125 

14.159 

91.069 

23.943 

90 

7.456 

92.487 

15.973 

92.417 

26.975 

92 

8.287 

94.068 

17.739 

93.983 

29.910 

93 

9.096 

95.855 

19.453 

95.755 

32.744 

95 

9.882 

97.834 

21.112 

97.719 

35.472 

97 

10.642 

99.993 

22.714 

99.861 

38.091 

99 


Table 7-6 Coordinates of the Profile of the Cutter to Machine the Tri-Lobed Extruder Sc 

Since there are three flutes for this screw, the full helicoidal surface w 
obtained in three steps. The blank will be indexed by 27t/3 radian every time 
machining the first and second flutes. 

73.4 Design of an End Milling Coffer for ffae Machining of a Quadri-1 

Exfruder Screw 

The cross-sectional profile of a bi-lobed extruder screw is similar to the one s 
in Figure 6-13. It has a tip radius of 100 mm. and a root radius of 90.743 mn 
mentioned earlier, only the root and flank arcs influence die desifft and the tip i 
realized in the preparation of the blank. The dimensions of one root arc and Ae two 
arcs adjacent to it are given in Table 7-7. 


/ 

P, 

9, 

r, 



15 

-76.604 

-64.279 

190.743 

0.39269 

0.69813 
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16 

0.000 

0.000 

90.743 

0.69813 

0.87264 

1 

-64.281 

-76.603 

190.743 

0.87264 

1.17808 


Table 7-7. Biparametric Definition of the Root and Adjacent Flank Arcs of 
a Quadri-Lobed Screw. 

Using the method described above, the profile of the cutter was calculated f 
pitch p=15 mm., /j= 150 mm. and p=115 mm. and the plots of these three profiles a 
shown in Figure 7-6. The coordinates of the cutter profile is listed in Table 7-8. 


p = 75 mm. 


V = 225 mm. 

Xf 

V, 

X, 

V, 

X^ 



-8.275 

99.996 

-16.211 

99.930 

-23.534 

99.68 

-7.606 

98.315 

-14.891 

98.259 

-21.596 

98.05 

-6.923 

96.783 

-13.543 

96.735 

-19.622 

96.56 

-6.223 

95.407 

-12.167 

95.368 

-1.7.614 

95.22 

-5.509 

94.195 

-10.768 

94.165 

-15.574 

94.05 

-4.783 

93.156 

-9.345 

93.133 

-13.508 

93.04 

-4.048 


-7.903 

92.279 

-11.417 

92.21 

-3.302 

91.620 

-6.445 

91.609 

-9.307 

91.56 

-2.550 

91.134 

-4.976 

91.127 

-7.182 

91.10 

-1.793 

90.841 

-3.498 

90.838 

-5.048 

90.82 

-1.032 

90.743 

-2.016 

90.742 

-2.909 

90.73 

-1.033 

90.743 

-2.015 

90.742 

-2.908 

90.73 

-0.826 

90.743 

-1.612 

90.742 

-2.327 

90.74 

-0.620 

90.743 

-1.209 

90.743 

^ -1.745 1 

90.74 

-0.413 

90.743 

-0.806 

90.743 

-1.163^ 

90.74 

-0.207 

90.743 

-0.403 

90.743 

■Ell 

90.74 
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-0.001 

90.743 

0.000 

90.743 

0.000 

90.743 

0.207 

90.743 

0.403 

90.743 

0.581 

90.743 

0.413 

90.743 

0.805 

90.743 

1.163 

90.742 

0.620 

90.743 

1.209 

90.743 

1.744 

90.741 

0.826 

90.743 

1.612 

90.742 

2.326 

90.740 

1.032 

90.743 

2.015 

90.742 

|||[||^^ 

90.738 

1.032 

90.743 

2.014 

90.742 

2.908 

90.737 

1.792 

90.840 

3.497 

90.837 

5.047 

90.825 

2.549 

91.133 

4.975 

91.127 

7.181 

91.102 

3.301 

91.620 

6.445 

91.608 

9.306 

91.567 

4.047 

92.295 

7.903 

92.278 

11.416 

92.216 

4.783 

93.155 

9.345 

93.132 

13.507 

93.046 

5.509 

94.195 

10.767 

94.164 

15.574 

94.051 

6.223 

95.406 

12.167 

95.367 

17.613 

95.224 

6.922 

96.782 

13.542 

96.735 

19.621 

96.559 

7.606 

98.315 

‘ 14.891 

98.258 

21.595 

98.049 

8.275 

99.995 

16.211 

99.929 

23.534 

99.685 


Table 7-8 Coordinates of die Profile of the Cutter to Machine the Quadri-l.obed Extruder Screw 


Since there are four flutes for Ais screw, Ae full helicoidal surface will be 
obtained in four steps. The blank will be indexed by Jc/2 raAan every time while 
machining the second, Aird and flutes. 


7.4 Determination of the Surface of Ae Cutter from Ae Profile of the 
Cutter "giu) 

The procedure to obtain Ae two-dimensional profile of Ae cutter ^g (*) from Ae 
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given cross-sectional profile («) of the helicoidal surface was explained earlier and 

this profile was determined for four extruder screws in the preceding section. This is the 
basic shape required to obtain various types of cutters in order to manufacture the 
helicoidal surface. One can obtain a turning tool, an end milling cutter, a side milling 
cutter or a disc type cutter by suitably manipulating this profile. In this section, the 
method of obtaining these different cutters and the respective manufacturing schemes are 
explained. In all these cases, the blank used will be a cylinder of radius equal to the tip 
radius of the cross-sectional profile. 


7.4.1 Machining of the Helicoidal Surface Using a Turning Process 

The simplest process to produce a helicoidal surface is turning since any turned 
surface is composed of a number of helices. The axis of the blank coincides with the Z 
axis (longitudinal axis) of the machine tool. The profile of the cutter is ground on the 
turning tool such that the plane of this profile is tilted by an angle of o with respect to 
the X - Z plane of the machine tool as can be seen in Figure 7-7. The translation of the 
cross slide (i.e., movement of the turning tool) is coordinated with the spindle rotation to 
get the required pitch p. The manufacturing scheme for this case is shown in Figure 7-7. 


7.4.2 Machining of the Helicoidal Surface Using a Milling Process with an End 
Milling Cutter 

Figure 7-8 shows the end milling cutter used for the machining of the helicoidal 
surface. The surface of the end mill is obtained by revolving the cutter profile about the 
F ^ axis. While machining, the axis of the end milling cutter shall be radial to the cylinder. 
The tip of the end mill will be at a distance equal to the root radius of the cross-sectional 
profile at the final depth of cut. 


7.4.3 Machining of the Helicoidal Surface Using a Milling Process with a Disc "Pype 
Cutter 


Many a time, milling with a side-and-face cutter is preferred to end-miltog 
to higher rate of metal removal. When a high surface finish is required, it may 
to grind the helicoidal surface after milling using a disc type of 
surfaces of the side-and-face cutter and the grinding wheel for machiMg^ Wioudai 
surface can be obtained as surfaces of revolution of die cutter jho e ® 


axis 
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Figure 7-9. Machining of the Extruder Screw Using a Mil ling process with a Disc Type Cutt ^ 
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Figure 7-10. Machining of the Extruder Screw Using a Milling process with a Cylindrical Cutter 
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parallel to the axis as shown in Figure 7-9. The axis of this disc cutter should be at an 
angle equal to the lead angle o with respect to the axis. 


7.4.4 Machining of the Helicoidal Surface Using a Milling Process with a Side 

Milling Cutter by Scan-Milling 

Scan-milling is a milling process by which the desired surface is realized by 
employing several passes of a cutter rather than a single pass using a form-cutter. TTiis 
method gives a surface which is an approximation to the desired surface since it produces 
scallops on the surface. This method is employed extensively in the machining of 
sculptured surfaces such as the dies of car bodies and aircraft components using end 
milling type cutters. By suitably choosing the distance between adjacent passes and by 
employing dressing and deburring, sufficiently accurate surfaces can be obtained. 

Since it is not possible to design a side mill which produces the helicoidal surfaces 
encountered in the extruder screws in single pass, these surfaces may be produced using 
cylindrical side mills or even by a wire-EDM process by a scanning process. In diis 
process, the axis of the wire or the side mill is at an angle ct with respect to the X axis 
and it lies in a plane parallel to the Xg - Zg plane as shown in 7-10. Generally, rough 
machining is preceded by scan-milling. A profile wdiich is offset from the cutter profile 
by an amount equal to the radius of the cutter is shown in the figure. Initially, the center 
f)f the cutter is at one end of this offset profile and at this condition, the blank moves 
from one end to the other helically. After one pass, the cutter is moved by a small amount 
along die center line profile and then the blank repeats the helical motion but in the 
opposite direction. This process is repeated till the cutter reaches the other end of the axis 
profile. 


7.5 Design of Side Milling Cutter for Machining the Helicoidal Surface in Single 

Pass 

The scan-milling is an approximate method and it also takes more time to 
machine. Hence it is preferable to machine the helicoidal surface in single pass. This is 
possible for all kinds of helicoidal surfaces when we employ a turning process or a 
milling process either with an end milling type cutter or a disc type cutter as describ«l 
in SectLons 7.4.1, 7.4.2 and 7.43, Tbete are some helicouial ffluf a ces fiat which, a side 
milling cutter exists using which the helicoidal surface can be produced exactiy ^ 
pass. In this section, the procedure to determine such a side milling cutter is Ascribed (i 
it exists) and the procedure is illustratwl by means of m example. 




Figure 7-11. Machining of the Extruder Screw Using a Milling process 

- Isometric View — 
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7.5.1 Procedure for the Design of a Side Milling Cutter 


In all the cases discussed earlier, the profile 0 *) of the cutter was obtained at 

the plane, i.e., at w = 0. This profile can be obtained for different values of w as 

well. The axis of the side milling cutter that is being designed will be parallel to the - 
plane and it will make an angle a with the X^ axis. The procedure to obtain the 
surface of this cutter is as follows : 


1. Obtain the cutter profile for various values of w at sufficiently 

close intervals. 

2. Choose an appropriate value for d where the axis of the cutter passes 

through a point (0, d, pv'Hit). From the axis of the cutter, find the 
minimum distance of each of the cutter profiles obtained in the last step. 
Draw circles each with center at (w cos o, d, pvl2% - w sin o) and the 
respective minimum distance as radius. These circles are drawn parallel to 
the plane. 

3. The cutter surface will be the surface that envelopes all these circles. 


7.5.2 Illustrative Example 


Consider the following generatrix ^8 («) for a helicoidal screw \\4iose pitch is 
60% and lead angle is %/4. 

>2 1] P-W) 


= 2 

72 = 4m 
Zj = 3 m 

-6 :£ M S 6 


where. 
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Let the axis of the side mill be at 30 units from the Jt; - Z, plane. 


The profile of the cutter (w) was calculated from w = -30 to w = 30 at an 

interval of 5 units. The resulting profiles are plotted and are shown in Figures 7-1 1 and 
7-12. The axis of the cutter is parallel to the axis and passes through the point (0, -30, 
0) with respect to the coordinate system. The shortest distance between the axis of the 
cutter and each of these cutter profiles are given in Table 7-9. 


w (=zj 

Min. dist.of the profile 
from cutter axis 

-30.0000 

24.9967 

-25.0000 

27.1884 

-20.0000 

28.9538 

-15.0000 

30.3042 

-10.0000 

31.2501 

-5.0000 

31.8131 

0.0000 

32.0000 

5.0000 

31.8131 

10.0000 

31.2501 

.15.0000 

30.3042 

20.0000 

28.9538 

25.0000 

27.1884 

30.0000 

24.9967 


Table 7-9 Coordinates of the Profile of the Side Milling Cutter for the 
Example Given in Section 7.5.2. 

From Table 7-9, a two-dimensional profile can be obtained in the plane at 

-30.000 if w is plotted along axis and the minimum distance along X^ axis. When 
a curve fitting was done, it was found that these points define a circle with a radius of 































Figure 7-12. Machining of the Extruder Screw Using a Milling process with a Side Miiiing Cutter 
-Front View 
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67.7857 units and center at (*35.7857, -30.0000, 0) in - Z, plane. The side mill is a 
surface of revolution obtained by revolving this profile about the axis of the cutter. This 
side milling cutter and the desired surface are depicted in Figure 7-13. 

The side milling cutter described above exists only for simple surfaces such as the 
one shown in this example. Since most of the extruder screw surfaces are more 
complicated, the side milling cutter to machine them in one pass does not exist. Hence, 
for realizing an extruder screw surface in single pass, one has to resort to one of the other 
three processes, viz., a turning process, a milling process using an end milling cutter or 
a milling process using a disc type cutter. 



Chapter 8 


CONCLUSIONS AND FUTURE EXTENSIONS 


8.1 Concluding Remarks 

Today, the speed of computers and their graphic rendering capabilities with the 
help of high resolution display devices have improved so much that we are able to see 
the objects more realistically on the screen. What is called Virtual Reality in the parlance 
of Computer Animation has become possible. These advancements have triggered research 
in the geometric modeling of motion simulation of robots and mechanisms such as the 
aircraft landing gear mechanism. Extensive research is also being done by several 
researchers in modeling the metal removal processes geometrically (Wang, Diysdale etc.). 

The geometric modeling of shape realization processes is an important aspect for 
a manufacturing engineer from the viewpoint of design for manufacturability. In order to 
be able to evaluate the feasibility of realizing a surface as defined by a designer, a 
manufacturing engineer is required to define the geometry of the tool, the kinematic 
structure of the machining process and the law of relative motion between the tool and 
the work piece. In many cases, the surface realized is not exact but is a close 
approximation of the surface defined by the designer. In such cases, it is always desirable 
to explore different feasible machining strategies. The model proposed in the present 
thesis provides a methodology for depicting the geometry of the machined surfaces for 
the alternative machining strategies. This has been accomplished using unified symbolic 
models and procedures. The proposed methodology is a useful predictive tool for selecting 
an appropriate machining strategy for realizing a specified surface. 


8.1.1 Conjugate Geometry Formulation 

The present thesis outlines the methodology of developing geometric models of 
several shape generating processes. The interaction of the geometric parameters of the 
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cutter surface as well as the kinematic parameters of the manufacturing process 
significantly affect the geometrical definition of the manufactured surface. It has been 
found that the concept of conjugate geometry enables us to distinguish between the 
mathematical definition, the display definition and the shape realization definition of a 
surface. The shape realization definition, as can be seen from the examples and the case 
studies, is different from the other two definitions. Generally, a design engineer is more 
familiar with the mathematical definition (implicit or explicit form) and the display 
definition (parametric form) of a shape. The sh^e realization definition (envelope form), 
it is felt, is useful in designing for manufacturability of a surface. This definition can be 
used for designing a machining process, the tool geometry, a feasible tool path and 
ensuring a smooth machined surface. 

The concept of conjugate geometry has been used to define analytically four 
methods of realizing either a generatrix or a directrix curve of a machined surface. This 
classification enables us to define about twelve different groups of shape realization 
processes. One can see that almost all conventional as well as modem machining 
processes can be characterized with the help of these twelve basic groups. The proposed 
conjuagate geometry formulation enables the modeling of all tiiese twelve groups of 
manufacturing processes in a unified manner. 

As compared to the early ^sterns such as the APT system wherein a single fixed 
coordinate frame is used, in the conjugate geometry formulation, three coordinate frames 
are used. One frame is fixed and the other two are moving coordinate frames each being 
attached to the blank and the cutter. The movements of the cutter and the blank are 
considered independent of each other. The ^nchronization of their movements is achieved 
by suitably describing the twelve parameters of motion in order to obtain the desired 
surface. This mediodology is particularly advantageous since on many machine tools it 
is not only the movement of the cutter but die combined movement of both die cutter and 
the blank that produces die surface on the blank. Employing the homogeneous 
transformation matrices for transforming the positicm, normal and velocity vectors from 
one frame to the odier makes die calculaticms simple and elegant. 


8.1.2 Generic Cutter 

The generic cutter definition is very elaborate and it covers a wide ^ectrum of 
cutting tools used in shape realization by metal removal. The nine parameters describing 
the generic cutter can also be used to store the information of cutters employed in a 
manufacturing shop in a data base precisely. 
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8.1.3 Symbolic Manipulation 

The use of symbolic manipulation technique for developing conjugate geometry 
model of the machine tool and the generic model of the cutter for the shape realization 
process has been found to be an effective modeling and computational tool. The ^mbolic 
algorithm proposed here is general in nature and can be used for modeling any shape 
generating process. Since the algorithm, is implemented symbolically, just by changing the 
nine parameters describing the cutter surface and the twelve motion parameters, one can 
obtain the algebraic expressions describing the surface produced on the blank. This is very 
helpful when a manufacturing engineer has to evaluate various alternatives to manufacture 
a surface as he or she is freed from deriving the expressions for each case. 


8.2 Scope for Future Work 

The theory developed in this thesis can be used in developing a user-friendly 
menu-driven software. The work makes some basic assumptions. It assumes that the 
cutter, work piece and the elements of the machine tool are rigid bodies. While cutting 
is taking place, it is assumed, the stresses generated are local and do not affect the 
neighborhood. If the elastic properties of different members of a machining process are 
taken into account, then the surface actually machined on die blank will be different from 
what is obtained using the present modeling. The proposed method for the cutter path 
simulation can be linked with an analysis package in order to amve at the optimum NC 
code. These extensions of this thesis work are explained in die following sections. 


8.2.1 Development of an Integrated Package Using the Conjugate Geometry 

Formulation 

In its present implementation, die ^mbolic algorithm proposed in this thesis 
accepts only one expression for each of the twelve motion parameters describing the 
motions of the cutter and the blank. However, in reality, the surfaces machined are 
composite and hence they cannot be expressed by one algebraic equation. For each 
segment of the machined surface, there will be a different set of values for the motion 
parameters. Therefore an integrated software is required to be developed using the 
conjugate geometry theory proposed in the present diesis. Such a software will use the 
proposed symbolic algorithms to determine the geometry of the composite machined 
surface. 

The conjugate surface obtained in this work is nothing but the swept volume of 
the cutter. In other words, the cutter is assumed to be fully immersed in the blank 
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material. However, actually it may no be so. Hence, in order to get the actual surface, we 
shall have to trim the portion of this surfaces by the surfaces defining the initial shape of 
the blank. 

Presently, the algorithm determines only the surfaces and it does not give any 
indication as to which side of the surface the solid volume is present. One can determine 
this by following some conventions in defining sweep parameter v for the cutter profile. 
For instance, if v is positive in the counterclockwise direction about the Z, axis, then the 
solid material lies on the left side as one moves along the positive direction of v. By 
using similar conventions while defining the initial shape of the blank, we can obtain the 
required information that indicates on which side of the blank surface the materia! is 
present. 


8.2.2 Modeling the Effect of Tolerances 

No dimension can be produced exactly and hence each dimension has a tolerance 
limit within which it is produced. Each of the eight parameters describing the generic 
cutter will have a tolerance. There will be clearance on the spindle bearing which causes 
spindle run-out. These two tolerances will influence the geometry of the cutter surface Ej. 
Each of the links forming the kinematic chain of the machine tool will also have a 
tolerance. Similarly, there will be clearances on the bearings between the links. All these 
tolerances will influence the definition of the twelve motion parameters. However, For the 
modeling of a manufacturing process in this work, it has been assumed that all the 
dimensions are deterministic. Hence the actual machined surface Ej will be slightly 
different from the one predicted by the conjugate geometry model here. In order to obtain 
more realistic results, the effect of the tolerances of the cutter surface and the kinematic 
links of the machine tools will have to be taken into consideration using the stochastic 
approach. 


8.2.3 Modeling of Mechanistic Properties and NC Code Optimization 


What is considered in this thesis is the solution to a geometric problem. The 
surface of the cutter is considered to be a continuous surface without flutes which create 
the cutting edges. The mechanics of cutting is not dealt with. However, under the 
influence of the cutting forces, the cutter, the blank and the elements of the machine tool 
structure wilT have deflections. Due to the combiiied effect ot these detlections, the shape 
of the surface produced will be different from what is predicted by the present symbolic 
algorithm. It is required to ensure that the deviation of the actual surface is well within 
the tolerance requirements. Under extreme conditions, the deflection may exceed the 
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Figure 8-1 . Schematic Diagram of the System for NC Program Optimization 
Using the proposed Conjugate geometry Formuiation 
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elastic limit and cause the breakage of the cutter or even damage to the machine tool. It 
is essential to predict the cutting stresses beyond which this damage might occur so that 
preventive measure such as reduction of cutting speed, feed rate and depth of cut could 
be done. 

In order to take the elasto-dynamic behavior of the manufacturing process in to 
account, a system similar to the one shown in the block diagram of Figure 8-1 is to be 
used. The primary requirement of such a system is a solid modeler preferably of the B-rep 
type. Initially, the conjugate geometry algorithm receives the description of the cutter and 
the system is fed with the initial shape of the blank. The analysis module is fed with data 
related to the material properties of the cutter and the blank, the number, shape and size 
of flutes and the allowable deviation of the actual surface from the desired surface. After 
the initialization, the first two NC blocks are read and they are converted into the twelve 
algebraic expressions as a function of time and fed to the twelve degrees of freedom 
model. Using the definition of the cutter and these twelve expressions, the twelve degrees 
of freedom model calculates the swept volume. 

The swept volume of the cutter generated and the current shape of the cutter are 
sent to two modules, one is a trimming algorithm which trims the swept volume in order 
to get the actual shape of die metal removed from the blank; the other module is a 
difference Boolean operator which calculates the shape of the blank after the metal is 
removed. The output of this module will be used to update the current shape of the blank 
at the end of the cycle. 

The most vital part of this system is the analysis module. The input to the analysis 
module is the sh^e of die metal removed by the current NC block. The analysis module 
will calculate the deflection of die cutter at eveiy instant of cutting. The module will 
break the current NC block into suitable number of NC blocks and fit optim^ feed rate 
values for each of diem so that die cutter deflection is within the allowable limit. Thus, 
it converts each NC block ■wdiich may be unoptimal to several NC blocks which are 
optimized. 

There are two approaches for the analysis of cutting forces for this situation. One 
approach is the lumped mass approach. In this method, the analysis is done only at few 
representative points. All die forces that are induced during the cutting and the resulting 
deflections are transformed to these points by incorporating suitable factors. After this, 
the static and dynamic analyses are performed. This approach is approximate but will be 
computationally economical. The second approach is the finite element analysis. In this 
approach, the matrix of the volume is broken into very small elements and each of these 
finite elements are analyzed. This approach will give the most accurate results but this 
will be computationally intensive. 
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Since the conjugate geometry algorithm generates the surface as a first step, it will 
be more appropriate if the integrated software proposed is based on a B-rep solid modeler 
incorporating the Boolean operators. 
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Appendix 


Definition of the Model of the Generic 
Machine Tool - Listing of MACSYMA Program 


/* FILE NAME : MCHN.MAC */ 

y************************^ 


kill (all) ; 

/* Declaring Dependencies of variables */ 

/* on Parameters ul,vl,u2,v2, and t */ 

/* */ 

depends ( [xl, yl, zl], [ul, vl] ); 

depends { ix2, y2, z2], [u2, v2] ); 

depends { [aal^ bbl, ccl, thetal, phil, sil], t ); 

depends ( [aa2, bb2, cc2, theta2, phi2, si2], t ); 


/* 

/* 

/* 

pi 

/* 


CUTTER SURFACE: 


Defining 
: matrix 
Position 


/ 

====== */ 

Point P on Cutter 
{ [xl, yl, zl, 1] ) ; 
Transformation matrix 


Surface in Ml system */ 
from Ml to MO */ 


r_thetal_xl 

: matrix { 

[ 1 , 

0 , 

0 , 

0 ] . 



[ 0 , 

cos(thetal) , 

sin (thetal) , 

0 ] , 



[ 0 , 

-sin{thetal) , 

cos (thetal) , 

0 ] , 



[ 0 , 

0 , 


1 ] ) ; 

r_phil_yl 

: matrix ( 

[cos (phil) 

. 0 , 

-sin(phil} , 

0 ] , 



[ 0 , 

1 , 

0 , 

0 ] , 



[sin (phil ) 

, 0 , 

cos (phil } , 

0 ] , 



[ 0 , 

0 , 

0, 

1 ] ) ; 

T_sil_zl 

: matrix { 

[cos (sil) , 

sin(sil) , 

0, 

0 ] , 



[-sin (sil) 

, cos (sil) , 

0. 

0 ] , 



[ 0 , 

0 , 

1 , 

0 ] , 



[ 0 , 

0 , 

0, 

1 ] ) ; 

tr_01_to__00 

: matrix { 

[ 1 , 

0 , 

0, 

0 ] , 



[ 0 , 

1 , 

0, 

0 ] , 



[ 0 , 

0 , 

1. 

0 ] , 



[aal. 

bbl. 

ccl. 

1 ] ) ; 


tlO 


r_thetal_xl . r_phil_yl . r__sil__zl . tr_01_to_00; 


/* Position Transformation matrix from MO to Ml */ 
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tOl : tlO^^-1; 

tOl : trigsimp (tOl) ; 

/* Transformation matrix for normal from Ml to MO: 


mlO 

: matrix { [tlO[l,l], tl0[l,2], 

tl0[l,3]. 

0], 



[tl0[2,l], tl0[2,2], 

tl0[2,3], 

0], 



[tl0[3,l]. tl0[3,2], 

tl0[3,3]. 

0], 



o 

o 

0, 

0]) ; 

/* Transformation matrix for normal from MO 

to Ml: */ 


mOl 

matrix { [t01[l,l], t01[l,2]. 

t01[l,3]. 

0], 



[t01[2,l], t01[2,2], 

t01[2,3]. 

0], 



[t01[3,l], t01[3,2]. 

t01[3,3]. 

0], 



[ 0, 0, 

0, 

0]); 

/* Velocity Tranformation Matrix from Ml to MO: */ 


tmpll 

'diff 

(tl0[l,l], t); 



tmpl2 

•diff 

{tl0[l,2]. t); 



tmplS 

'diff 

(tl0[l,3], t); 



tmpl4 

•diff 

(tl0[l,4], t); 



tmp21 

‘diff 

(tl0[2,l], t); 



tmp22 

‘diff 

(tl0[2,2], t); 



tmp23 

•diff 

(tl0[2,3], t); 



tmp24 

'diff 

{tl0[2,4], t); 



tmp31 

•diff 

(tl0[3,l]. t); 



tmp32 

•diff 

{tl0[3.2], t); 



tmp33 

•diff 

{tl0[3,3], t); 



tmp34 

•diff 

(tl0[3,4], t); 



tmp41 

•diff 

{tl0[4,l]. t); 



tmp42 

•diff 

(tl0[4,2], t); 



tmp43 

•diff 

(tl0[4,3], t); 



tmp44 

•diff 

(tl0[4,4], t); 



tlOv : 

matrix 

( [tmpll, tmpl2, tmpl3, 

tmpl4] , 




[tmp21, tmp22, tmp23. 

tmp24] , 




[tmp31, tmp32, tmp33. 

tmp34] , 




[tmp41, tmp42, tmp43. 

tmp44] ) ; 


/* Velocity of P in 

Cutter Surface w.r.t. MO: 

*/ 


vOpl : 

pi . tlOv; 



/* BLANK SURFACE: * 

/ 



/* =-= 

- * 




/* Defining Point P 

Dn 6lank Surface in M2 system */ 


p2 : 

matrix 

( [x2, y2, z2, 1] ) ; 



/* Transformation matrix from M2 to MO */ 



r theta2 x2 : 

matrix 

( [ 1, 0, 

0, 

0] , 



[ 0, cos{theta2) 

, sin{theta2), 

0], 



[ 0, -sin {theta2) 

, cos{theta2); 

0], 



[ 0, 0, 

0, 

1]); 

r_phi2_y2 : 

matrix 

{ [ cos{phi2), 0, 

-sin (phi2 ) , 

0}, 



[ 0, 1, 

0, 

0], 



[ sin(phi2), 0, 

cos {phi2 ) , 

0], 



[ 0, 0, 

0, 

1]) ; 

r si2 z2 : 

matrix 

( [ cos(si2), sin{si2); 

0, 

0], 



[-sin{si2) , cos (si2) , 

0, 

0], 
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tr_„02„to 00 


[ 

0, 

0, 

1, 

0], 

[ 

0, 

0, 

0, 

1]) 

[ 

1. 

0, 

0, 

0], 

[ 

0, 

1, 

0, 

0], 

[ 

0, 

0, 

1, ■ 

0], 

[ 

aa2. 

bb2. 

cc2 , 

1]) 


t20 


r_theta2_x2 . r_phi2_y2 . r_si2_z2 . tr__02_to_00; 


Transformation matrix from MO to M2 */ 
t02 : t20"^"^-l; 

t02 : trigsimp (t02); 


/* Transformation matrix for normal from M2 to MO: */ 


m2 0 

: matrix ( [t20[l,l]. 

t20[l,2]. 

t20[l,3]. 

0], 


[t20[2,l]. 

t20[2,2] , 

t20[2,3], 

0], 


[t20[3,l]. 

t20[3,2]. 

t20[3,3]. 

0], 


[ 0, 

0, 

0, 

0]) ; 

/* 

Transformation matrix for normal from MO 

to Ml: */ 


m02 

: matrix ( [t02[l,l], 

t02tl,2]. 

t02[l,3]. 

0], 


[t02[2,l], 

t02t2,2]. 

t02 [2,3] , 

0], 


[t02[3,l]. 

t02[3,2]. 

t02[3,3]. 

0], 


[ 0, 

0, 

0, 

0]); 


/* Transformation matrix for normal from MO to Ml: */ 
1 "*^ Velocity Tranformation Matrix from M2 to MO: */ 


tmpll 

,diff 

(t20[l,l]. 

t) ; 


tmpl2 

•diff 

(t20[l,2]. 

t) ; 


tmpl3 

•diff 

{t20[l,3], 

t); 


tmpl4 

•diff 

(t20[l,4]. 

t) ; 


tmp21 

•diff 

{t20[2,l]. 

t) ; 


tmp22 

•diff 

(t2D[2,2], 

t) ; 


tmp23 

•diff 

(t20[2,3]. 

t); 


tmp24 

■diff 

(t20[2,4]. 

t) ; 


tmp31 

•diff 

{t20[3,l]. 

t); 


tmp32 

‘diff 

{t20[3,2]. 

t) ; 


tmp33 

■diff 

{t20[3,3]. 

t); 


tmp34 

■diff 

{t20[3,4], 

t); 


tmp41 

■diff 

(t20I4,13, 

t) ; 


tmp42 

•diff 

<t20[4,2], 

t); 


tmp43 

■diff 

(t20[4,3]. 

t) ; 


tmp44 

■diff 

{t20[4,4]. 

t); 


t20v 

matrix 

( [tmpll, 

tmpl2. 

tmpl3, tmpl4] , 



[tmp21, 

tmp22, 

tmp23, tmp24] , 



[ tmp3 1 , 

tmp32, 

tmp33, tmp34] , 



[tmp41. 

tmp42 , 

tmp43, tmp44] ) ; 


/* Velocity of P in Blank Surface in MO: */ 
v0p2 : p2 . t20v; 


v0p2 : p2 . t20v; 

/* CUTTER SURFACE NORMAL: */ 

/* = = = = = = = = = = = */ 

/* Normal of Cutter Surface in Ml: */ 
dxl_dul : ’diff (xl, ul) ; 

dyl.dul : *diff (yl, ul); 

dzl^dul : 'diff (zl, ul); 
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' di f f { xl , vl ) ; 

•diff (yl, vl); 

'diff {zl, vl); 

matrix { [dyl_dul * dzl^dvl - dyl_dvl * dzl_dul, 

dxl_dvl * dzl_dul - dxl_dul * dzl„dvl, 

dxl_dul * dyl_dvl - dxl_dvl * dyl__dul, 

0 ]) ; 

Magnitude of Normal of Cutter Surface: */ 
magnitude_nlp : sqrt (nlp[l,l]^2 + nlp[l,2]'^2 + nip [1 , 3 ] ^^2) ; 

/* Unit Normal of Cutter Surface in M2:*/ 
unit_nlp : matrix ( [nlp[l , 1] /magnitude_nlp; 

nip [1,2] /magni tude_nlp , 
nip [ 1 / 3 ] /magni tude_nlp , 

0 ]); 


dxl_dvl 

dyl_dvl 

dzl__dvl 

nip 


/* Unit Normal of Cutter Surface in MO: */ 


nOpl : unit_nlp . mlO; 

/* BLANK SURFACE NORMAL: */ 

/* = = = = = = =: = = =: = = = = = = = = = = = */ 

/* Normal of Blank Surface in M2: */ 

dx2_du2 : 'diff (x2, u2); 

dy2_du2 : 'diff {y2, u2) ; 

dz2_du2 : 'diff (z2, u2) ; 

dx2_dv2 : 'diff {x2, v2); 

dy 2_dv2 : 'diff { y 2 , v2 ) ; 

dz2_dv2 : 'diff {z2, v2); 


n2p : matrix { [dy2_du2 * dz2_dv2 ~ dy2_dv2 * dz2_du2, 

dx2„dv2 * dz2_du2 - dx2_du2 * dz2_dv2, 

dx2_du2 * dy2_dv2 - dx2_dv2 * dy2_du2, 

0 ]); 


/* Magnitude of Normal of Blank Surface: */ 

magnitude_n2p : sqrt {n2p[l,l]^2 + n2p[l,2]'^2 + n2p[l,3]^2); 

/* Unit Normal of Blank Surface in M2:*/ 
unit_n2p : matrix { [n2p[l, 1] /magnitude_n2p, 

n2p[l, 2] /magnitude_n2p, 
n2p [1,3] /magni tude__n2p , 

0 ]); 


/* Unit Normal of Blank Surface in MO: */ 
n0p2 : unit_n2p . m20; 

/* CONDITION OF CUTTING: */ 

/* = = = = =: = =: = = =: = = = = = = = =: = = = */ 

/* Relative Velocity in MO: */ 
v0pl2 : v0p2 - vOpl; 

constraint! : v0pl2 . transpose (nOpl ) = 0; 
constraint! : trigsimp (constraint!) ; 

constraint2 : v0pl2 . transpose {n0p2 ) = 0; 
constraint2 : trigsimp {constraint2 ) ; 


save { "mchn. sav" , all); 


Appendix B 


Definition of the Model of the Generic 
Cutter - Listing of MACSYMA Program 


^★★★★★★★*************** j 

/* FILE NAME : CU.MAC */ 

kill (all) ; 

/* NOMENCLATURE */ 

/* ============ */ 


/* ull : value of ul at point P for the generic cutter */ 
/* ul2 : value of ul at point Q for the generic cutter */ 
/* ul3 ; value of ul at point S for the generic cutter */ 
/* ul4 : value of ul at point T for the generic cutter */ 

/* si ; value of ul at point P for the current cutter */ 
/* s2 : value of ul at point Q for the current cutter */ 
/* s3 : value of ul at point S for the current cutter */ 
/* s4 : value of ul at point T for the current cutter */ 


/* GENERIC CUTTER */ 

/* */ 

e<in_OcP : Zc = Xc * tan (a); 

eqn^PQ : (Xc - e) "2 + (Zc - f ) "*2 = r^2; 

®«3n_QS_when__b_is_0 : Xc = Zc * tan{b) - d/2 * tan (a) * tan{b) +d/2; 

eqn_QS_for_others : Zc = Xc * cot{b) + d/2 * (tan(a) - cot(b)); 

eqn_ST : Xc = h * tan{b) - d/2 * tan{a) * tan(b) + d/2; 

Point_Oc : [0,0,0]; 

XcA : e; 

ZcA : f; 

Point_A : [XcA, 0, ZcA] ; 

XcR : d/2; 

ZcR : d/2 * tan(a); 

Point^R : [XcR, 0, ZcR] ; 

dist_AR : sqrt { (XcA-XcR) ^^2 + (ZcA-ZcR)"^2 ); 

XcS : h * tan{b) - d/2 * (tan{a) * tan(b)) + d/2; 
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ZcS : h; 

Point_S : [XcS, 0, ZcS]; 

convert_a_b_to_radians {) := BLOCK( a : a * %pi / 180, 

b : b * %pi / 180) ; 

calculate_e_f {) := 

BLOCK (if Is_Fillet_Radius then 

(e : d/2 - r * {cos(a) - sin (b) ) /cos (a+b) , 
e : ev{e) , 

f : (e * sin{a) + r)/cos{a), 
f : ev(f) ) ); 

calculate_h () := 

BLOCK (if b = 0 and Is_Fillet_Radius then h : f, 

if b = 0 and Is_Fillet_radius = FALSE then h : ZcQ) ; 

Xc(ul) :=BLOCK { 

if (ul >= 0 and ul < si) 
then ( 'ul * cos (a) ) 
else if (ul >= si and ul < sl+s2) 
then 

if r = 0 then {'ul * cos (a)) 

else if r = INF then XcP + ( (XcQ-XcP) /s2 ) * ( ' ul-sl ) 
else (e + abs(r) * cos ( atan2 (sl*sin (a) -f , sl*cos(a)-e) 

+ Cul-sD/r)) 

else if (ul >= sl+s2 and ul < sl+s2ts3) 

then ( ( *ul-sl-‘s2“s3 ) *sin (b) + ev{XcS) ) 
else if (ul >= sl+s2+s3 and ul <= sl+s2+s3+s4) 
then '(ev(XcS))); 

Zc(ul) :=BLOCK ( 

if (ul >- 0 and ul < si) 
then ( *ul * sin (a)) 
else if (ul >= si and ul < sl+s2) 
then 

if r = 0 then ( 'ul * sin(a)) 

else if r = INF then ZcP + ( (ZcQ-ZcP) /s2 ) * ( ‘ ul-sl) 

else (f + abs(r) sin ( atan2 (sl*sin (a) -f , slices (a) -e) 

+ ( 'ul-sD /r) ) 

else if (ul >= sl+s2 and ul < sl+s2+s3) 

then ( { *ul--sl“'s2-s3) *cos (b) + ev(ZcS)) 
else if (ul >= sl+s2+s3 and ul <= sl+s2+s3+s4) 
then { { 'ul-sl-s2"s3) +ev{ZcS))); 

Xl{ul,vl) := Xc{ul) * (cosCvD); 

Yl(ul,vl) Xc(ul) * (sinCvD); 

Z1 (ul, vl) := Zc (ul) ; 


save ( "cu . sav" , all ) ; 
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/* FILE NAME : CUTTER. MAC */ 


nuirier : TRUE; 

convert_a_b_to_radians { ) ; 


/* 

CALCULATION 

OF E & F; */ 

/* 


=-==-===- */ 

calculate_e_f { 

) ; 

/* 

POINT P: */ 


/* 

==—==— */ 



checkup 

: ( (e + f*tan(a)) * (cos (a) "'2) ) ''2 

-(e-'2 + f''2 - r''2) * (cos(a)''2); 


checkup 

: ev{ checkup); 


if abs (checkup) < 0.000001 then check__P : 0 


else if asksign (checkup) = 'NEG 


then 

checkup : POINT_P_INVALID; 


XcP 

: (e + f*tan(a)) * (cos (a) ^2) - sqrt (check_P) ; 


ZcP 

: XcP * tan (a) ; 

/* 

POINT Q: */ 


/* 

==- */ 



check_Q 

: ((f + (e + (d/2) * {tan(a) *tan(b) -1) ) *tan{b) ) 

* (cos{b)"2) )'^2 

- {(f'"2 + (e + (d/2)*(tan(a)*tan(b)-l) )"‘2 - r‘^2) 

* (cos(b)"^2) ) ; 


check_Q 

: ev{check_Q); 


if abs{check_Q) < 0.000001 then check__Q ; 0 
else if asksign<check_Q) = 'NEG 

then check^Q : POINT_Q_INVALID; 

ZcQ : (f + (e + {d/2 ) * { tan (a) *tan (b) -1 ) ) * tan(b)) 

* (cos(b)'^2) + sqrt (check_Q) ; 

XcQ : ZcQ * tan(b) ~ (d/2) * ( tan (a) *tan (b) -1) ; 

/* CALCULATION OF H */ 

/* =============:==== 

calculate__h (}; 

/* CALCULATION OF SI, S2, S3 & S4: */ 

si : sqrt (XcP^2 + ZcP'^2) ; 

si : ev{sl) ; 

if abs(b) = 0 then s3 ; 0 

else s3 : sqrt { ( (h - d/2* ( tan (a) -cot (b) ) ) *tan (b) - XcQ) ^2 

+ {h - ZcQ)"‘2 ) ; 
s3 : ev(s3) ; 

if r = 0 then s2 : 0 

else s2 : r * ( atan2(ZcQ - f, XcQ - e) - atan2(ZcP - f,XcP - e) ) ; 
s2 : ev{s2); 

s4 : hi; 
numer : FALSE; 




